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USING THE METHOD OF MOMENTS TO CALCULATE 
THE SPACE-ENERGY DISTRIBUTION OF NEUTRON DENSITY 
FROM FLAT AND POINT SOURCES IN AN INFINITE MEDIUM 


A. P. Ptitsyn 


Translated from Atomnaya Energiya, Vol. 10, No. 2, pp. 117-126, February, 1961 
Original article submitted April 6, 1959 


A method is given for calculating the space-energy distribution of neutron densities from flat and point 
sources in an infinite medium. 


N 
The neutron density (x, E) is sought in the form E)= Sa; K - To a large 


i={ 


degree the form of the function K (x) is arbitrary; its setection is based on physical principles. From 
the 2N space moments of the function (x, E), 2N parameters aj, bj dre found, The neutron density 
distribution is found in hydrogen and water. The calculations for water are compared with experi- 
mental data. A comparison with the accurate solution of Wick [1] in the case of retardation of 
neutrons by hydrogen shows that from four moments the suggested method can, with sufficient ac- 
curacy, find the spatial distribution of neutrons at distances up to 20 free path lengths. 


In [1-4] in principle the problem was solved of finding the space-energy distribution of photons and neutrons in 


an infinite medium, In [1,2] a study was made of the asymptotic form of the solution for very large distances from 
the source. 


An analytical solution was found in [1] for the case of a constant free path length for strongly retarding neutrons. 


For cross sections depending on energy, in [3] a semianalytical method was developed to calculate the distribution at 
very large distances from the source, 


For comparatively small distances the method of polynomial expansions can be used for the calculation [4]. It 
is precisely these distances (~15-20 free path lengths) which are of interest to us. The idea of the polynomial method 
consists of expanding the required function into a series in terms of known polynomials Uj(x): 


p(x, E)=e (E) U; (Bx), (1) 


where the coefficients A;(E) are found from a system of integral equations. This method was developed in application 


to photons,and at distances of 15-20 path lengths it can find the photon density by means of four polynomials, which 
amounts to finding four space moments. 


In the case of weakly retarding media this method is also suitable for determining neutron density; however, 
with strong retardation, especially when the free path length decreases weakly with decrease in energy, the series in 
terms of polynomials converges poorly. Here the use of four terms of the series [1] at distances equal to 5-6 paths leads 
to negative values of neutron density, This is presumably due to the unsuitably of the polynomial expansion method 


for these problems, although a knowledge of the four space moments y p (E) is in itself sufficient to find the neutron 
density at the given distances for any values of E*. 


As is known, neutron density mainly depends exponentially on the distance x; however,the exponent is not al- 
ways a linear function. Thus, in the case ofa one source 6(E—Ep) the neutron density, where there is not a 
single collision, at all values of x is proportional to e~* r, Being retarded, these primary particles will give at any 
distances a neutron spectrum proportional to e-*/ . The accumulation of neutrons, as shown in [1,2], leads to a 


* The author of [5] showed that a value of the parameter 8 can be chosen in the polynomial expansion (1) for which 
the series in terms of polynomials converges well. 


109 


still slower reduction in density at large distances. On the other hand, at small distances x< In E)/ E, according to 
the theory of age,the density of neutrons retarded to an energy E(E< Eg), is proportional to ¢—**/‘*, 


An attempt could be made to find the density » (x, E) in the form of a linear combination of functions of the 
given type 
p(x, E) =a, ™ 4 ag (2) 


However, in this case to determine the coefficients aj, bj from the spatial moments of the function y, a system of non- 
linear equations is obtained of a type such that the solution cannot be obtained in a clear form. 


We will seek the function yp (x, E) in the form 
N 
p(x, E)= a, (F)K 2] (3) 
i=1 


(with regard to the form of the functions K (x) see section 3). The advantage of this form of approximation is that ap- 
proximations of any order can be obtained in the same way,and the main advantage is that the system of equations 
for the coefficients aj, b; permits an accurate analytical solution. As will be shown, two terms of the sum (3) describe 
with sufficient accuracy the neutron density for any energy E at distances of x to 20 free path lengths. 


1. Integral Equations for Moments 


To obtain equations which areesatisfied by the spatial moments of neutron density in the case of a flat mono- 
energy source we will proceed, as usual, from the kinetic equation 


a(u)p’ — da, du'ds; (u’) \ dQ’G; (uo, u— u’)p’ (2, vu’, p’)= 5 (x) d(u) S (), (1.1) 
0 


where u= In E)/E; Ep is the energy of the source neutrons;, is the cosine of the angle between the axis x and the 
direction of movement of the neutron; {9 is the cosine of the scatter angle of the neutron; »*(x. u, #)=vn(x,u,H); 
v is the speed of the neutron; n(x,u, «)dx du dQ is the number of neutrons in an element dx du dQ; 


— 


S(#) is the angular distribution of neutrons leaving the source; o (u) and Og; (u) are the macroscopic cross sections 
(total and elastic scatter at nuclei with a mass Mj,respectively); aj = (Mj + 1)*/ 4Mj. 


For simplicity the scatter is assumed to be elastic and isotropic in a system of a center of masses. Extension to 
the case of inelastic and nonisotropic scatter is elementary. Nonisotropicity is allowed for by introducing after the 
integral sign the factor w; [(Mj+ 1)*/2M4)x e~(u-u 1y where (cos $) is the indicatrix of scatter 


1 
at nuclei of the type i, standardized 7— \ @; (cos $)dQ=1. Inelastic scatter is allowed for by introducing the inte- 


gral of nonelastic collisions 
u 
0; (w')W;(u', u) po (2, u’)du’, 
0 
inel, 
where 0 ‘. is the cross section of inelastic scatter; Wj(u’, u) is the probability of inelastic scatter from u’ to u, 


standardized so that W,(u’, u)du=1, 


To eliminate peculiarities in the right side of equation (1.1) we will separate from the function p*(x, u, M) the 
part which describes those neutrons of the source which have not undergone a single collision with nuclei, i.e. we 
will represent ~* in the form 


a 
4 u 
ae 1 
G; (Hg, u) == In ( 
oe 
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The function v(x, #) apparently satisfies equation (1.1) without 4n integral of collisions 


tr 1.3 
The solution of this equation is 
t 
(x, =P (ux), (1.4) 
where 9 (y) is a known discontinuous function 
1, y>0, 

8(y)= | 0, y<0. (1.5) 


We will multiply equation (1.3) by Py (#) [P, (#) Legendre polynomials) and integrate for all angles 2. We 
obtain 


tr tr 
OY, 
21-+1 [ 1) +1 | +0 (0) = S). (1.6) 


Here y and are the coefficients of expansion of the functions and in terms of Legendre polynomials 


5, = (w). 


t co 
Le gett. \ 


(1.7) 
where 8 is a certain constant. Multiplying equation (1.6) by (6 "+ Vat) x" and integrating with respect to x from - 
to + ©, we find for the coefficients mz, the following systems of algebraic equations: 


Mig = 
min = | 1) + |, (1.8) 
nmi, 2, ... 
where 


The system of equations (1.8) is recurrent. Each coefficient with an index n is expressed only through a coef- 
ficient with an index n-1. The system is solved in the following order: from a determination of the coefficients (m;, = 
=y Sz ) it follows that 


my = [+ o |. 


Hence it is possible to find mz,,etc. 


The coefficients mp can be obtained directly by substituting in the integrals (1.6) and (1.7) the function » t 
from (1.4). 


We will deal with the system of equations (1.8) in detail since this will be useful at a later stage. 


Assuming similar to (1.7) that 


tr 

| 
| 
| 
nti 
Min (x, ude. (1.10) 
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Substituting ¥'(x, u,u) (x, x (x, 6(u) in equation (1.1), multiplying it by (B"*2/n!) a"P, (u) » and 
integrating with respect to dQ and dx, we obtain, considering the determination of (1.10) and (1.7), the system of 
integral equations which we require to determine the functions Myn (u)*: 


u 


(u) M,, (u)— 2a, \ 6 + e—(u-u’) P, (u,)M,, (u’) du’ = x 
1 0 


x + 1) Miss, n-1+1M Mine (0) P, (ue)6 (q; — u); 


M,, (u) = 0, 


Hi = =0, qi = In 
(1.12) 
The system of integral equations (1.11) for M,, and also the system (1.8), is recurrent and is solved in exactly 


the same order. 


The initial values of the functions M,, can readily be found. In fact, for u= 0 the integrals in the left side of 


(1.11) revert to zero, the Lengendre polynomials in the right side revert to unity (4 j(0)= 1) and for Mj, (0) we obtain - 
a system of algebraic equations 


(0) M,,, (0) = (1+ 1) Migs. (0) + min (0). 


Since the coefficients mjp satisfy the system of equations (1.8) we find the solution of the system (1.13) [6]: 


(0) 
(0) =(n+ 1) m,,, 


(1.14) 


For the direct solution of equations (1.11) the coefficients mj, must be found in a clear form; in their turn these 
depend on the coefficients $7 of angular distribution of the source S (i). 


Since the kinetic equation (1.1) and hence the system of equations (1.11) are linear, the solution of the latter 
for an arbitrary source 


(1.15) 
S*(u)=S*(—p), = —S-(—p), 


can be found as the sum of solutions corresponding tosymmetrical St (u) and antisymmetrical § (4) sources. Therefore 
*In the case of nonisotropic scatter instead of Pj(u) in equation (1,11) it is essential to substitute w (cos v) P7(u). Al- 


u 
lowing for inelastic scatter leads to an additional integral term (u’) W; andtoa 
i 


corresponding free term d10MipZo wel (0) W; (0, u). 
1 
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. 
(1.11) 

where 
u—u’ u—u’ 
M;+41 2 M;—1 

‘3 8) =-— ; 

(1.13) 

where 


we will later assume that the source S(u) is either an even or an odd function of 1. Without limiting the generaliza- 


tion it is apparently sufficient to find the coefficients min and my, for the symmetrical 


= [8 — +9 10) (1.16) 
and the antisymmetrical 


= | — Ho) — 8 (4+ My) (1.17) 
40 
5 -sources, respectively. 


As A eimunigiay above, the coefficients mJ can be found not only ffom the system (1.8) but also directly, substi- 
tuting yt (x, MH) in (1.7). Taking as S(u) (1.16) we find 


n+ 


(the value of y is determined by formula (1.9), Hence 


O for odd (3+n), 
Min = 


y"*!P. (Uo) for even (8 +n), 
Similarly, using (1.17) we find 


(Mo) WO for odd (b+ n), 


Min = 
™ 0 for even (2 +n). (1.20) 


We will consider one particularly important case of a symmetrical source, where the angular distribution & () 
is independent of the direction (isotropic source); 


(1.21) 


For the coefficients mj, according to (1.21) and (1.8) we have [5] 


2nd 


0 for odd(n-2)orn<l, 


i 
a n 2!'n! ! 
2 (4) dy ay ) if ) is a non- (1.23) 
!(n-4-1-4-4)! negative even number. 


An important equality follows from (1.23) and system (1.11) as well as the obvious relationship Mjp (u)= 0, if (2 +n) 
is an odd number. 


M,,(u)=0, if - (1.24) 


2. Point Source 
The “moments” MP in the case of a point isotropic source 


(1.19) 
where 
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ned t 
ve (u) = \ (r, u)4nr?dr (2.1) (2.1) 
0 
for even values of 2 + ncan be found knowing the moments Mjp for a plane isotropic source. In fact, with regard to 
the known connections between the solutions for plane and point sources [2], 


u) = 2n \ drrwp Pf (r,u) P, (+) (2.2) 
Ix} 
We have [5] 
M,,, (u) = 2n ( dxx™ ( drrypP* (r, u) (+) =%,,M (2.3) 
|x| 


In particular, for moments pt the neutron density 
= (2k +1) Mo, on k=0, 1, 2, 


3. Method of Approximating a Function With Respect to Its Moments 
Let there be 2N space moments 


x"g (x) dx 


of a certain function g(x), given for the range, (0, ©), We wish to approximate the function g(x) by the sum 


N 


g (x) ~ 2 a,K (b,x), (3.1) 


= 


where the function K (x) is such that its moment 


K,= \ x"K (a) dx 


exist. In [3] a method is given which can be used to express the coefficients aj, bj by moments gy. This method is 
essentially as follows. We will write the sum of (3.1), somewhat changing the designations, in the form 


N 
(x)= (ni/s) K (2/84). (3.2) 
Multiplying (3.2) by x" and integrating from 0 to ©, we obtain a system of nonlinear equations 


N 
= > (3.3) 


or, dividing (3.3) by Ky and introducing the designations fn =g,/Kp: we obtain 


+ NoSg+ ... + NySy = fii 


Let f (s) be a function such that the right sides ca of the system (3.4) are its moments, We then obtain this system of 
equations if in the N-point quadratic formula of the Gauss type with weight f (s) 


(3.4) 


oo N 


i=1 
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we substitute instead of H (s) successively all powers of s from zero to 2N-1, for which the equality in the approximate 
expression (3.5) is fulfilled accurately, As is known, sj in expression (3.5) is the zero of the polynomial hy(s), ortho- 


gonal with the weight f (s) to all polynomials of lower powers, and 1; are the so-called Christoffel numbers, deter- 
mined by the formula 


dhy(s) 1-10... hy (s) 


ds 8; 
i 


Therefore, the solution of the system of equations (3.4) is given by formula (3.6) and by the roots of the equation 


(3.7) 
hy (s) = 0. 

If the function f (s) is positive, then all roots of equation (3.7) are real and different. In the general case among 
the roots sj there can be some which are related in a complex way. Let sq and sq, 4 be two such roots. Then the two 
complex related terms in formula (3.2) will give 


(Na/Sa) K + (MNa+1/Sa41) K (2/sa41) = 2Re [(Na/Sa) K (2/sq)]. (3.8) 


Among the roots of equation (3.7) there can also be multiple roots. We will consider this case in greater detail. 


Let a certain root sy have a multiplicity m. To be precise we will assume that these are the roots 


§,=5,=...=s, = 6"). 


Formula (3.2) in this case then becomes 


™ N 


g (a) = Na) K (Ba) + (ni/si)K (2/s;), (3.9) 
i=m-- 
m 
and the Christoffel numbers  . according to (3,6) revert to infinity. It is true that their sum ( > Na) is finite and ac- 
a=1 


cording to the first equation of system (3.4) is equal to 


m 


however, the function 


N 
BCK (Bx) + (n;/s;) (2/s;) 


in the general case will not be a solution of the problem. Generally speaking, 2N moments of this function do not 
coincide with 2N given moments g,, of the function g(x). This is not a contradiction of system (3.4) since formally 
this system is “satisfied” by finite additions arising on the addition of infinite terms 


m 


Nasa = +- arbitrary constant. 


The solution of the problem in this case will be the function 


2m N 
(0) =K Br) (mi/si) K (3.11) 
n= 


In particular, if all N roots of equation (3.7) are equal, then according to (3.11) for g(x) we have 


2N 
§ = K (Br) a 


= N 

Na= fo n=C, (3.10) 
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and this is a polynomial expansion [1]. 


Returning now to the problem of finding the neutron density it can be stated that at the distances we are con- 
cerned with the neutron density behaves mainly as K (6 x) (for example as eB X), then the solution can be found in 


the form of a polynomial expansion (3.12). If the density at small and large distances is very different (different 
exponents) then the solution must be sought in a general form (3.2). 


We will write in a clear form all formulas for the case when the four moments gp, g), 82, gg Of the function 
g(x) are known, The coefficients of the polynomial 


ha (s) = s?—as +b, 
of the orthogonal with weight f (s) to all polynomials of a lower power, i.e. to 1 and s, are equal to 


fafo—fehi — FB 


The roots of equation h,(s)=0 are 


= 5 [a — (a2 — 


(3.14) 


5 [a+ (a? — 4b)"]. 


co 
a (6) dem 


(3.15) 
= fo— 


We will use the obtained formulas to approximate the neutron density ¥(x, u) using the moment ¥o,(u). If the 
angular distribution of the source S (4) is symmetrical, then the neutron density p(x, u) is an even function x. If S(K) 
is antisymmetrical, then the function %(x,u) is odd. It is therefore sufficient to approximate it only for positive values 
of x, i.e. for the range 0= x< ©; toward negative values of x the function p(x, u) should be continued symmetrically 
or antisymmetrically, For an even value of ~(x) the function K (x)should therefore be taken in the form 


(3.16) 
and for an odd value, in the form 
(3.17) 


For even functions from the system of integral equation ( 1.11), different-from-zero even moments are obtained, 
and for odd functions, odd moments are obtained. Therefore, in the range (0,%) for the required function wp» (x) we 
will only know whether the moments are even or uneven, So that the above-described method of approximation of 

(3.1) can be used in this case, the even and the odd functions must be sought in the form 


» p(x) even ; 


N 


dd 
v(x) (x/Vs;) ad (n i/S; Vs; ye odd, 


Then for 7 ;, s; all the formulas derived above are valid if fp is taken to mean totitn and y fine 1 The 
functions Mj, (u) according to (1,10) were determined by the formula 
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| 
Hence 
| 
2s, - 


M 


in = \ (x, u) dx, 


which in view of the symmetry (or antisymmetry) ~7 (x) can be rewritten in the form 


co 


nti 
M,(u) = 2 \ op, (x, u) da = Pin (4), (3.20) 


where ~zp represents the nth space moment of the function 7 (x) in the range (0,). From (3.20) 


! 
Vin = er My. (3.21) 


dd 
According to (3,16) and (3.17) for 4 1 Wwe have 


Ki, = (2n)!; 
K on = (2n + 2)! (3.23) 


Using (3.21)—(3.23), for moments fn fron (3,4), determining the values sj and nj, we obtain 


Mo on 


3.24 

Po, ant Mo, 2n 

K » Wo (x) odd. (3.25) 


2n--1 


We will write the corresponding formulas forthe neutron density from a point isotropic source. According to 
(2.1) and (2.5) we have 


pt anti 
Mo,2n = u) dr =F an (u) 


2 1)! 


where 


u) =4ar%p,(7, uv). 


Hence, for moments f,, we obtain 


Pon 2 1 
in = Mo, 2ns (3.27) 
Kon 
and for the density po (r, u) 
N 
P(r, u) = (7, u) = (n/V's,) (3.28) 
i=1 


If the roots s, and s, are connected in a complex way then, for example, to approximate the even function p(x) 
with respect to its four even moments we obtain, according to (3.18) 


1/4 


sin sin cos — 


| 

| 
or 
| 

(3.26) 
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-=, 4b—a? ; 

1 1 w 
A= Ta ( —ysin 
1 1 w, 
B= (4 fosin $+yeos 5); 


The values of a and b in (3,30) are given by formulas (3.13). 


(3.30) 


4. Neutron Density in Hydrogen and Water 

Hydrogen, We will find the solution of the system of equations (1.11) for a plane isotropic source in hydrogen, 
assuming the free path length of the neutron to be constant (which is the case in the range of energies 10*-1 ev), 
(u)=0,(u)=8 = 1. We will introduce instead of (¥) the functions Nin (u)=e4M7, (u). According to (1.11) the 
functions Nj, (u) are satisfied by the equations 


u—u’ 


Nin (u) — P, (e Nin (u ) du = 41) Nias, n-1 +IN nil + (e—*/2): (4.1) 


) 


N,, -1(u) =0. 


To find four moments of neutron density it is necessary to solve the following ten integral equations of system (4.1): 


1) Noo— \ No au’ = 1, u 


2) Ny—Ve Ny, du’ Nog (u) + 


7) (e [Se-("— 4") — 3] Ngg du’ = 


3) Noo \ No du’ =Ny(u)+3 5 8) Nu—> \ — 1] Nog du’ = 
0 


4) Noo — \ [Se-(—u") — 1] N,, du’ = == (3N 33 + + ge — 1); 
0 uo uw’ 
== Ny, (u) + (Be -- 1); ) Nis \ e N,, du 


u 


u—u’ = 


2 Nygdu’ = (2N Ng) + 
1 10) Nos 


1 


By differentiation with respect to u all these equations become differential equations and are solved in an 
elementary way. In fact, from the first equation by differentiation with respect to u we obtain 


| 
F 
3 = No» + e7u/2 (5e-¥ 3); 
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Noo (4) — Nog = 0, Noo (0) = 1. 


No (2) 


From the second and third equations we find 
N,, (u) = 4 e¥/?; 
We will consider the fourth equation: 
u 
Nyy (u)— J — 1] Nag da’ = = Nyy (u) + (Be~™ — 1) = 
0 


Differentiating twice with respect to u, we obtain 


” 1 2. 
5 (u) +9" (uw) = — 


(0) = (0) = 2 ’ 


(0) = 9 (0) + 9' (0) = 


Nyx (u) = — 2 4 sin (w/V’2) 
22 (U) = Fe ri + 75 sin u/ ° (4.5) 
Omitting such simple calculations, we will write the final solution for the moments Mo, 9k(u): 


M = 1; 
Mos (u) + 


M (u) = 
— 
as “(cos 


In [1] for u= 10 the accurate space distribution of density (x) was found. We will calculate it for u= 10 from 
equations (4.6) and compare it with the Wick solution. 


From formulas (4.6) we find 


M,, (40) = 4, Ma (10) = 65,578, 
(10) = 10,667, = 278,44. 


Hence 
(4.2) 
(4.3) 
(4.4) 
Hence 
V3 | (4.6) 
+ [ 1 — 
3 
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TABLE 1. Comparison of the Results for Determining Neutron Density by Different Methods 


Ig Wo (x) — 1,140 | —1,046 | —1,065 | —1,229 | —1,346 | —1,773 | —2,102 | —3,052 | —4,254 
from Eq. (4.7) 


lz (x) -~1,046 | —1,057 | —1,089 | —1,240 | —1,340 | —1,741 | —2,100 —3,070 | —4,250 
‘from ick's solu- 


tion 


The relationships (3.13) and (3.24) give 
a = 8,7289, b= 27,558, 4b—a?= 34,038 > 0; 


i.e. the roots s; and sy are connected in a complex way. According to (3.30) 


5 =5,8342, A=0,03598, y= —0,54010, w= —0,5892, B= —0,25725. 
Substituting the values in (3.29), we obtain the formula for the neutron density 


Wo (x) = 1.07196 cos (0,1267Bx) + 0,51450 sin (0,1267B2)}. (4.7) 
In Table 1 the solution of (4.7) is compared with the Wick solution. 


Water. The solution of the system of equations (1.11) for a plane isotropic source of neutrons with energy Ep = 
= 2.5 Mev in water was found numerically, Formulas (3.27) and (3.28) should be used to transfer to a point source. 


In experiments carried out by I. S. Pogrevov (1957) the distribution of density was measured by a chamber from 
u™®, surrounded by a cadmium cover of thickness 0.5 mm, To obtain the distribution of the number of reactions R(r) 
in the chamber we find the moments of the function R(r). According to (2.5) we have 
15,7 


RE =(2n+ 1) Rok! =(2n-+1) | Mo,2n(u) o7°* (u) du. (4.8) 


Hence, according to (3.27) and (3.28) 


Here, R(r) is the number of reactions in the chamber, per one atom of u™* and one neutron of the source. We will 
introduce the designation 


B (r) = er4nr*R (r), (4.10) 


1 


=A(0)=5 cm, 


In Table 2 the results of calculation using formulas (4.8)-(4,10) are compared with experimental data. 


Unfortunately, due to the insufficient source strength the measurements could not be taken for distances r 
greater than 36 cm, 


In conclusion the author would like to thank Yu. A. Romanov and L. P. Feoktistov for their valuable discussions, 
L. S. Bychenkova and M. N. Kopylova for performing the numerical calculations. 


| 

4 

a 

4 2n+1 pPl 

and 

a where 
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TABLE 2. A Comparison of Calculated and Experimental Results 


1,72 | 2,22 


| 3,22 | 4,22 


10°4B 
from equation 


(r) 
(4.10) 


94,4 131 


210 294 


1024 B (r), 
mental data 


experi- 


G. Wick, Phys. Rev. 75, 738 (1949). 
U. Fano, J. Res. Nat. Bur. Standards 51, 95 (1953), 
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THE CREATION OF A MAGNETIC FIELD WITH 
AN AZIMUTHAL VARIATION 


R. A. Meshcherov and E. S. Mironov 


Translated from Atomnaya Energiya, Vol. 10, No. 2, pp. 127-130, February, 1961 
Original article submitted June 22, 1960 


A method is given for calculating the form of pole surfaces for creating a magnetic field with an azi- 
muthal variation of given depth and with a certain law of change of the mean field intensity along 
the radius. Experiments are described on the modeling of a magnetic field of a 14 -meter cyclotron 
with azimuthal variation of the magnetic field. 


Calculating the Shape of Pole Surfaces 
To calculate the shape of poles from a given distribution of intensity of magnetic fields in the middle plane 
H, = F(r,¢) a method can be used which is based on the assumption that pole surfaces are magnetic equipotentials. 


As for axially symmetrical fields [1], the scalar magnetic potential in a cylindrical system of coordinates can be re- 
presented in the form of a series 


(r, 2) 2 Vor_1 (7, 9) (1) 


For the coefficients of these series the following recurrent relationship holds, obtained from the Laplace equation: 


Obviously 
oP 


Using formulas (1)-(3) it is possible to calculate the shape of pole surfaces for an arbitrary distribution Hz(r,¢), which 
can be given both analytically and graphically. In particular, this method can calculate poles for creating a field 
H,, = Ho(1 + Ar’ +Br cos ng), first suggested by Thomas [2]. However, the poles of the magnet in this case have a very 
complex profile, and their preparation involves considerable technological difficulties. 


The variation of the magnetic field can be achieved much more simply (by means of plane sector covers [3]). 


There are no accurate methods for calculating the magnetic field of 
these poles. The problem was therefore solved on the assumption 


+, 


er that the magnetic field intensity between the middle of the sector 
~ 20 LL and the middle of the recess depends on ¢ the same as in the case of 
Peale Se --y  Plane-parallel poles, the shape of which is shown in Fig. 1. This as- 
S q sumption obviously holds only under the condition 
TT. Ra >d. R( = a ) >d 

Fig. 1. Shape of plane-parallel poles used to where R is the radius; a is the angle of the sector; n is the number 
calculate the distribution of magnetic field of sectors, 

intensity. 


Figure 2 gives regions on complex planes which are converted 
into one another by means of the formula 


Vv. <n Din d44— 


= (4) 


: 
: 

| 

d 
w=—In 

x 
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where x=d/D. 
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Fig. 4. Dependence of the magnetic field intensity, 
averaged with respect to ¢, on the radius, 


Fig. 2. Regions on complex planes w and z. 


Hyy(u,0)/H max | Hz/Hz max 


Fig. 5. Shape of sector covers 
and boundary of equivalent 
magnetic field (dotted lines). 


From the theory of the complex potential it follows that the 
magnetic field between stepped poles (see Fig. 1) is equal to 


H (w) =H,,+iH,= — tHmax 


-2 -1 0 1 a d e x? 
Fig. 3. A comparison of the calculated and experi- where Hmax is the field intensity on moving an infinite distance 
mental data. The continuous curve represents the to the left. 


calculated data; the points show the results of the 
measurements at different radii (in millimeters): 
@—50; x—70; O—90; A—110; O- 130; @-145, 


Using formulas (4) and (5) a calculation was made of the 
distribution of intensity H, =f (u) in the middle plane between 
the stepped poles of infinite length. The results of the calcula- 
tion are given in the form of a curve in Fig. 3. 


In many cases, especially when studying the movement of ions in an idealized magnetic field of the sector 
type [4], Hy =f (u) is best represented in the form of an equivalent stepped field (see Fig. 2). A characteristic value 


of this field is the parameter A—the difference in coordinates of the steps of the pole surface and the equivalent 
magnetic field: 


Hy jv 
A,-H 
Hmax 5 
| Hy (u,0) 2 
min 50 100 200 
3 
-5 
y @ f 
x : 
| 
40 
| 
\ 
0.8 
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0 co 


(6) 


Preliminary Experiments 

The experiments were carried out with an electromagnet with pole pieces of diameter 370 mm and 90 mm 
between them, Between the cylindrical poles there were two steel discs of diameter 370 mm and thickness 16 mm 
with sectors fastened to them at an angle of 52.5° and thickness 5 mm. The gap between the sectors was 40 mm, con- 
sequently d= 20 mm, D= 25 mm and x=0,8. The measurements were carried out with a magnetic field intensity in 
the center Hy = 6000 oe. The points of Fig. 3 show the values of the function H,(R,¢YHz max Obtained from experi- 
mental data, Hzma x being the greatest of all values of magnetic field intensity near the middle of the sector. As 
should be expected, the best coincidence is observed at large distances from the center. At R= 50 mm the experi- 
mental data differ considerably from the calculated data. 


In the range of radii 50<R< 150 mm it can be assumed that 


d 
D ~ H:sec’ 


where Hz sec and H, , are the values of the magnetic field intensity near the center of the sector and near the center 
of the recess respectively for R= const. Consequently, in this range the boundaries of the equivalent magnetic field 
created by the sector covers will be straight lines not passing through the origin of the coordinates. The magentic 
field intensity averaged with respect to ¢ therefore decreases smoothly with increase in R from 50 to 145 mm (curve 
1 of Fig. 4). The boundaries of the equivalent magnetic field are shown in Fig. 5 by dotted lines and the boundaries 
of the sector coveres are shown by continuous or bold lines, 


Using formula (6) it is possible to calculate the shape of the sector covers in the range of radii 50<R< 145 mm 
for any relationship of H,(R). For example, to create H,() =const the sector covers were reduced parallel to their 
edges by 4 mm on each side since for 4 mm on each side since for x=0.8 and d= 20 mm Aw 4mm, Measurements 
of the magnetic field with this shape of covers showed that in the range of radii 50<R< 145 mm F(R) const (see 
curve 2 of Fig. 4). However, for small values of R, when the described method of calculation is unsuitable, the mean 
value of the field intensity, as before, differs considerably from H,, for large values of R. To even out the field the 
apexes of the sectors were cut off and discs of 40 mm diameter were placed in the center, The thickness of the discs, 
calculated from the condition Hy =H,,, equalled 


6,=D— = 2,5 mm. 


d 
m 


To make the field smooth a two-stage transition was made from the sectors to the central disc (the continuous 
bolder lines of Fig. 6). 


Figure 4 (curve 3) shows the distribution of the mean 
value of magnetic field intensity obtained with this 
center. As can be seen from this diagram, H,, is constant 
with a good approximation. 


The shape of the sectors was then changed so that 
(broken lines of Fig. 6) the mean value of the magnetic 
field intensity increased with increase in R. The increase 
in the mean value of the field intensity AH, = f (R) due 
to change in the shape of the sectors is shown in Fig. 7. 
The difference between the experimental and calculated 
data (~0.1%) is explained by the inaccuracy in preparing 
Fig. 6. Shape of covers providing a given relationship the sectors and errors in the measurements. 
of H,(R) radii O< R< 145 mm. For H,(R)=const the 
boundaries of the sectors are shown by continuous bolder 
lines; for H,(R)—a monotonously increasing function 
—the boundaries of the sectors are shown by broken lines. 


Modeling the Magnetic Field of a 1}- 
Meter Cyclotron 

The magnetic field of a 1} -meter cyclotron was 
modeled on the mentioned electromagnet, the cylindri- 
cal pole pieces being replaced by tapered pieces of 300 
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mm diameter, Between the poles there were 300 mm diameter steel discs which served as a model of the covers in 
the accelerating chamber of the cyclotron, In order to create a magnetic field with a maximum depth of variation 
of approximately + 15% (allowing for the design of the chamber) we chode d=16 and D=22 mm, Three sectors were 
fastened on each disc with an angle of 60° and thickness D-d to 6 mm. Based on a value A for x= 0.73, each sector 
was reduced parallel to its edges by 4.4 mm on each side. The apexes of the sectors were cut off and in the center 
there was a disc of 42 mm diameter and thickness 3.5 mm, The measurements were carried out at Hy = 14,500 oe, 
therefore during the measurements the thickness of the central disc (6)= 3.2 mm) was made more accurate. 


2 


3 1 
0 75 


R, cm 
9 50 100 150 200 Fig. 9. The relationship H,(R) for the 1} -meter 
R, mm cyclotron (Hg = 14 000 oe), 


Fig. 7. Increase in mean value of field intensity 


AH, = f (R) due to a change in shape of sectors At > 110 mm there is a sharp drop in H,(R), caused by 


edge effects; therefore to retain the necessary shape H,(R) in 
the range of radii 110<R< 140 mm steel correcting elements 
were fitted. 


Since the distribution of the mean value of the field 
intensity, obtained as a result of the measurements, was in 

good agreement with the expected distribution, from the model- 
ing data lids were made for the chamber, sector covers and 
correcting elements for the 14 -meter cyclotron. Figure 8 shows 
the lid of the 14 -meter cyclotron chamber with azimuthal 
variation of the magnetic field. After final correction of the 
magnetic field, which was made by screws with an enlarged 
head, the required relationship of H,(R) was obtained (Fig. 9). 
Measurements of the field intensity on the cyclotron magnet 
and on the model gave such close results that no adjustments 
were needed to the sectors, the central discs and the peripheral 
correcting elements. 


Fig. 8. View of lid of 15 -meter cyclotron chamber. 
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THE THERMOELASTIC STRESSES IN THE WALLS OF 
A REACTOR HOUSING WITH INTERNAL SOURCES OF HEAT 
IN NONSTATIONARY STATES 


B. I. Maksimenko, K. N. Nikitin, and L. I. 


Translated from Atomnaya Energiya, Vol. 10, No. 2, pp. 131-137, pinta: 1961 
Original article submitted February 18, 1960 


Bashkirov 


In nonstationary processes, the thermoelastic stresses arising at various critical points and parts of 
a reactor, working over a wide range of change in heat load, can exceed the stresses arising under 
stationary conditions. Therefore, to ensure reliable operation with changing loads, the rate at which 
these processes takes place should be limited, The geometric form of a number of the critical points 
may be reduced to the two usual configurations of plane and cylindrical walls. By solving the non- 
stationary heat conduction problem for planeand cylindrical reactor walls with internal heat sources, 


comparative results have been obtained for the thermoelastic stresses arising in walls of this sort dur- 
ing transitional processes. 


In order to calculate the thermoelastic stresses in the walls of a reactor operating under nonstationary condi- 
tions, it is necessary to solve the nonstationary heat conduction problem under the following conditions; 1) the internal 
heat sources are uniformly distributed throughout the material of the wall; 2) the coefficient of heat conduction of 

the material is independent of temperature; 3) the magnitude of the thermoelastic stress does not exceed the elastic 
limit of the material, and the form of the wall remains constant; 4) the temperature distribution is one-dimensional. 


If there are internal heat sources, the maximum stresses, in the majority of cases, are found to be the tangential 
stresses at the surface of various parts. Therefore all the considerations which follow will be directed toward finding 
these stresses. 
Plane Wall 


The law governing the distribution of temperatures in a plane wall with internal heat sources (Fig. 1) is given 
by the well-known heat conduction equation [1]: 


()) 


We are solving this equation with the boundary conditions 


ot 


dy 


and the initial condition 


t = tos 


Here h= a@/) is the ratio of the coefficients of heat production and thermal conductivity; ty = ty + cT (°C) is the temper- 


ature of the coolant; tg (°C) is the initial temperature of the medium surrounding the wall, and c is the rate of change 
of the coolant temperature. 


We also assume that the temperature of the coolant varies linearly, while the strength of the internal heat sources 
remains constant with time (q, = const). In this case the expression for the temperature at any point of the plane wall 


takes the form: 
© 2sin By» cos ( Bux ) at 9 


With the same boundary conditions, but with the initial condition 


= 
Ot dy CpY 
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which corresponds to a stationary working state of wale qo and at a coolant temperature 


ty =t,+ct, where t, > ty, 


the solution of Eq. 1 will have the form 
wd? 2sin By cos (Bay ) 


where w = c-(qvkepy ) is the rate of change of temperature at any see of the wall, in the stationary state, and in the 
presence of internal heat sources. 


The tangential thermoelastic : stresses on the surface of a plane wall, 


2 of infinite length and thickness 26 , bathed by coolant on both sides, may 
a ae be expressed by the equation 
E 
o (t—#). (4) 
The expression for the mean temperature of the wall will then have 
re the following form: 


1) for the initial condition ¢|,—9 = ty 


=t,+ct— )+ 


Fig. 1. Diagram of a plane wall with 


3? 2sin® B, 5 
internal heat sources. Bn) exp ( — Bi 


t=t,-+er— 3 ) +) a exp ( — Br +) (6) 


In these expressions, 8 p is the nth root of the characteristic equation 


Depending upon the initial conditions, the tangential stress will be equal to 
a,E wd? 


(8) 


2 sin? B, at, 


2 sin By cos ( ) 


at 
(Bn-F Sin Bn Bn) exp( — Ba gr )- (10) 


| 

= 

= ‘ 

where 
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It follows from Eqs. (7) and (8) that: 


1) the thermoelastic stresses in the plane wall are direc- 


tly proportional to the difference w = C-(qy/cpy ) and the square 
Jum tr = of the half thickness of the wall, and are inversely proportional 


to the coefficient of thermal conductivity of the material: 


Cm | 
— — — — — — | 2) the maximum absolute value of the stress occurs at 
ae | | | the surface of the wall (at y=5); 

| | 
Gv, £2 | 3) with increasing time, the difference decreases, 
q to " l | and, accordingly, the absolute value of the stress at the surface 

| 

| te | tm-1| tm 1, 
0 


For sufficiently long times, and in particular for T = 
Fig. 2. Step-wise change in the coolant tempera- =[2+ (5/hd)](5 2/q ), the difference (®-%) becomes so small that 
ture, or in the strength of the internal heat sources. it may be neglected. In this case, theso-called stationary state 
of the second kind supervenes, in which the maximum value 
of the tangential stresses, for either of the two initial conditions considered, is given by the formula 


(11) 


This formula makes it possible to determine the thermoelastic stresses in a plane wall, with linear time varia- 
tion of coolant temperature, and with internal heat sources of constant strength. However, for practical purposes, if 
the thermoelastic stresses are being investigaged in the transientstate it. is necessary to consider the change in cool- 
anttemperature under the condition that the strength of the internal sources changes too. This can be done, using the 
heat conduction equation (1), for a step-wise change in the coolant temperature, or in the strength of the internal 
heat sources (Fig. 2). Leaving out the intermediate steps, in a way similar to the solution of Eq. (1) for a plane wall 


with linear coolant temperature change and internal heat sources of constant strength, we can finally write for the 
case of step-wise heating or cooling. 


where the value of ©, is given by Eq. (10); 


i y 
2 sin By cos (Br ) [ ] 
vont (Bn-t+ sin By cos Bn) 


2 sing, cos ( But ) 


cos Py, exp [ — Bi | and so on. 


and m is the number of steps or intervals in the process, 


The mean temperature across a section of the wall is given by the expression: 


1m 8? 2 > = 
t= ty (S475) [Wi Pr + — Wy) + (Wy OD (13) 


where the value of 4 is given by Eq, (9), and ads 


2 sin? a(t—t)]. 
Be BR (Basin Bn cos By) [ - ] 


tr 
a 
; 
‘ 1—p 2a 
Wmd 
= i 
; 
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2 sin? 8, 


Substituting Eqs. (12) and (13) in Eq. (4), we obtain the expression for the thermoelastic stresses 


As in the preceding cases, for T— Zz ™%> (2 _ a = all the differences (8- -®) become so small with 


i=1 
respect to the first term in the series that they can be neglected. 


As a result, the stresses may be calculated from the expression 


where, for the last step, Wry = Cy Cpy ). 


Examination of Eq. (15) yields the following conclusion: on going over to each succeeding step, the stresses 
increase (or decrease) very rapidly, and at the end of the step, for sufficiently large values of T;, T, etc, they re- 
main constant with time. 


From the general solution (14) it is not hard to find out how the maximum values of the stress vary at each step, 
for special cases involving various combinations of number of steps, change in strength of the internal sources, and rate 
of change of coolant temperature. 


Cylindrical Wall 


In cylindrical coordinates, the equation of heat conduction for a cylindrical wall of infinite length with internal 
heat sources takes the form: 


at at 1 at _ 
)+& ; (16) 


Cpy 


We shall solve Eq. (16) for two cases of heat removal. 


Case I. Internal heat removal with linear variation in coolant temperature and insulated external surface (Fig. 
3). 


The boundary conditions are 


= —h(t,—12); 


the initial condition is 


Fig. 4. External heat removal. 


Fig. 3. Internal heat removal. 


_ @E w 62 15) 
ot 
t | = Uo. 
aN 3 
1 
4 
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Here tr =ty +cT . With these two sets of conditions, the solution of Eq. (16) takes the form: 


2 
wr} 


t=t, +ct— 


[ t+ 1) — Ing | + =a o( Q; hry; k), (17) 


where k=r,/t,; P =1/ty, and rp, ty, and r are respectively the external and internal radii of the hollow cylinder and 
the variable radius of a point in between. 


The mean temperature over the thickness of the wall is given by the expression: 


Substituting the values of the ay gers: from Eq. (17) into the ey ane we obtain; 
3 1 _1 


(19) 
Subtracting Eq. (17) from Eq. (19) finally gives the ii of the temperature drops; 


1) on the external surface of the cylinder (r=r, and P=k) 


— 


2) on the internal surface of the cylinder (r=r, and p= 1) 


The thermoelastic stresses on the external and internal surfaces of the cylinder with internal heat transfer are 
equal respectively to: 


ext aE wr? 1 3 1 2k4 > int 


As in the case of the plane wall, the thermoelastic stresses produced by step-wise change in the strength of 
the internal sources, or coolant temperature may be represented in the form: 


On “4a (k, Q 


— + — (® 


m—1 
where m is the number of steps, or intervals, and 


{(k, (1-0) )- (1+—-)In k-+ Ing + k? — Q?. 


Comparing Eq. (24) with Eq. (14), we can draw the following conclusions: 


1) the changes in the thermoelastic stresses with multistep changes in the strength of the internal heat sources 
or coolant temperature occurring in the hollow cylinder are analogous. to those which occur in the plane wall; 


2) for sufficiently long time intervals, all the differences (®-) become so small, that they do not have to be 


A d 4a 
24 { 4 
\ +r ee (18) 
TL 
% (24) a 
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taken into consideration, and then the maximum values of tangential stress on the external and internal surfaces at the 
end of each step are given by the equations; 
aE wWyri 1 


int 27 3 
( zh 2 In k ). (26) 


Case II. External heat removal with linear variation in coolant temperature and insulated internal surface (Fig. 


4). 


The boundary conditions are: 


the initial condition is: 
= ly: 


Omitting the intermediate steps, we shall write the final equations for the thermoelastic stresses; 


1) on the external surface of the cylinder 


2) on the internal surface of the cylinder 


int FE wr? 1 k?Ink = int 


Experimental studies on the heating of a hollow cylinder, carried out in the Heating Laboratory of the Scientific- 
Investigational Tube Institute [2], have shown that a hollow cylinder can be considered as a rolled- up sheet of suitable 
wall thickness. This simplifies the calculation of the distribution of temperature and thermal stress in a hollow cylinder, 


From Eqs. (22) amd (23) for internal heat removal, and Eqs. (27) and (28) for external heat removal, the maxi- 
mum temperature drops entering into the expression for the thermoelastic stresses (4) are equal to: 


Aram, , here jhere my (k?+1)— 


If we consider the hollow cylinder as a rolled-up sheet with wall thickness 
6=r,—r,=r,(k—1), 


then, as follows from the solution (15), the maximum temperature drops entering into the expression for the thermo- 
elastic stresses (4) are equal to: 


1) on the surface bathed by coolant 


wr? 
a 


At, =n, , here ny = 


2) on the insulated surface 


at | = h(t, —t); 
ext aE wri 3 4 1 Ink ext]; 
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At, =n, here ny = (k —1)*. 


The quantities m and n may be given the name “form coefficients." Comparing the expressions for the tempera - 
ture differences in a hollow cylinder and a sheet, the conclusion can be drawn, that if the hollow cylinder is replaced 
by a rolled-up sheet, the thickness of the latter should be taken for calculational purposes to be 


=. (29) 


The values of the form coefficients m and n, and the ratio m/n, for different values of k=r,/r,, are given in 
the table. 


Values of the Form Coefficients and Their Ratios 


Values of k=1,/ry 


1,15 1,25 


m, 0,0016832 0 ,0037355 0 ,00667 0,0104 0,01495 0,04174 
0 ,003476 0,008061 0,01461 0,0234 0,03373 0,1015 
0,003165 0,007013 0,01218 0 ,018662 0,026365 0,0684 

0,001683 0003735 0 ,00667 0,0104 0,01495 0,04171 
—ny 0,00333 0,00751 0,01333 0,0208 0,03 0,0833 


0 ,00167 0,00376 0, 00666 0,0104 0,015 0,0416 


CONCLUSIONS 


1, The equations derived for thermoelastic stresses can be used for making calculations on heat generating 


elements, the reactor housing, and other points of interest during transitional states, as well as on starting and cooling 
off the reactor. 


2. The thermoelastic stresses at the insulated surface of a hollow cylinder, within the practical range of ratios 
of radii (up tos k= 1.5), may be calculated from the formulas for a plane wall. 


3. The thermoelastic stresses at the surface of a cylinder, bathed by a coolant, may also be calculated from 
the formulas for a plane sheet, but with the addition of a calculational wall thickness, given by Eq. (29). 
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THE REACTION BETWEEN SOLID UO, AND MnO, IN 
A SULFURIC ACID SOLUTION 


E. A. Kanevskii and V. A. Pchelkin 


Translated from Atomnaya Energiya, Vol. 10, No. 2, pp. 138-142, February, 1961 
Original article submitted July 5, 1960 


In connection with the extensive use of pyrolusite in the sulfuric acid leaching of uranium from ores, 
a study was made of the reaction between UO, and MnO, in a sulfuric acid solution,and possible mechan- 
isms were discussed. The experimental data show that the reaction apparently occurs at the points of 
contact of the hydrated surface layers of UO, and MnO,; the rate of the process is limited by steric 
hindrances, On this basis the role of iron ions was also considered in the reaction between a solid oxi- 


dizing agent (pyrolusite) and primary minerals of uranium in the acid leaching of uranium-containing 
ores, 


The nature of the reaction between solid UO, and MnO, in sulfuric acid solution is of considerable interest. It 
is known that at the present time pyrolusite is used extensively in the sulfuric acid leaching of uranium from ores [1-3] 
and it undoubtedly plays the part of an oxidizing agent. Processes occurring in the ore pulp containing a solid phase of 
complex composition (ore) a second solid phase (pyrolusite) and a number of ions and compounds in the liquid phase 
are very complex, For this reason the first stage in the investigation should be the comparatively simpler process—re- 
action between UO, and Mn0O,j in sulfuric acid solution. As yet this reaction has not been specially investigated, al- 
though there are some apparently correct ideas on the leaching of uranium from ores using MnO, (pyrolusite) in papers 
dealing with the technology of uranium leaching (for example see [3,4]). 


Starting Materials and Experimental Method 

Uranium dioxide was obtained from the mixed oxides which had first been subjected to peroxide purification 
by reducing with hydrogen at 900°C. The content of tetravalent uranium in the mixed oxide was 98% with respect to 
the sum of the tetravalent and hexavalent uranium, The MnO, used in the work was of the “pure” grade and was not 
subjected to further purification. The sulfuric acid was "chemically pure". The concentration of manganese in the 


solution after the completion of the process was determined colorimetrically [5]. The method for analyzing uranium 
was described in [6]. 


Possible Mechanisms 
As far as we know, there are no references in the literature to the stoichiometry of the reaction between UO, and 
MnOz, but this reaction is usually written in the form of the following equations: 
UO, -+ MnO, + 2H,SO, —> UO,SO, + MnSO, + 2H,0 
UO, + MnO, + 4H* —> UO2* + Mn** + 2H,0. 
Our experimental check of the correctness of this equation confirmed that the ratio of the number of UO, moles 


to thenumber of MnO, moles taking part in the reaction is equal to unity, For this reason in later work when consider- 
ing possible mechanisms we used this equation. 


or 


It might be assumed that in an acid medium at first one of the oxides is dissolved and then there is a hetero- 
geneous oxidation-reduction process, In this case the reaction of UO, and MnO, in acid medium is expressed by one 
of the following schemes: 


+ 4H fag) —> Ulagy+ 2H,0(1); (a) 


—> + + 2H,0c1 
+ MnO2,,, —> (b) 


MnOz, 4H (aa) —> y (a) 
UO2g)-+MnOx 5) + 4H faa) —> (11) | + —> Mnfig)-+ 


(I) 


; 


It can also be asumed that at first there is solution of UO, and MnO, and then homogeneous reaction of the tetra- 
valent and hexavalent uranium in the solution: 


4H fag) Ula) + 2H,O,1); (a) 
(111) + 4H (aq) 2H,0, 1) ‘b) 
4- + 2H,0q + + 4H (c) 


“(aq) 
It should be emphasized that in the suggested schemes the separate stages are not shown, but only the successive 
reactions, The possibility is not ruled out therefore that, for example, one of the stages of the reaction (II, b) or (III, c) 


is the formation of trivalent manganese, However, this problem will not be discussed here. 


A common feature of the schemes (1) and (II) is the fact that in not one of the successive reactions do the UO, 
and MnO, take part simultaneously. Because of this the rate of the total process should be in a linear relationship to 
the size of the UO, surface, if the (I, a) reaction is limiting. The size of the MnO, surface should not have any effect 
on the rate of the process. On the other hand, if reaction (II, a) limits the process, then its rate should be in a linear 
relationship only to the size of the MnO, surface. 


To check the correctness of these conclusions a determination was made of the dependence of the rate of reac- 
tion for a constant amount of UO, on the excess of MnO, with respect to the stoichiometric calculation (figure, MnO,/ 
/UO,) and then the dependence of the rate of the process on the excess of UO, (see figure, UO,/MnO,)* . An appropriate 
correction was introduced because the “pure” grade of manganese dioxide contains MnO impurity. 


However, the experimental data (see figure) show that the rate of the processes is affected by the size of both 
the UO, and MnO, surfaces. It should also be emphasized that the dependence of the rate of the process on the size of 
the UO, or MnO, surfaces is not linear, A similar shape of the curve presumably indicates that the UO, and MnO, play 
the same part in the reaction, i.e. they take part simultaneously in it, 


This proves that in the schemes considered above,not one of the 
30 -—""| stages can be considered limiting in any case, Furthermore, in the light 
U02 /Mn02 of obtained experimental data the possibility cannot be excluded that 
25 a the process could take place according to the following scheme, as- 
€ 4 suming “solid phase® reaction of the oxides: 
20 
: 15 — + + 
| (IV) U03,,) + 2H (a) > 1); (b) 
+ + 
10 Z Mn,5) -+ 2Hfaq) > 
: Yo The extensive use of MnO, to oxidize a number of substances 
S in neutral organic oxidizing agents at room temperature [7] apparently 
supports this assumption. 


5 10 15 20 25 


It is obvious that in aqueous solutions the reaction between UO, 
Ratio of components, M 


and MnOp, i.e.,stage (IV, a), is not a solid phase, reactionin the sense 
Rate of solution of UO, and MnOy,. in which it is understood in [8-10], In a number of other papers it is 
assumed that the so-called solid phase reactions in most cases proceed 
with the decisive participation of gases or liquids or the same and others simultaneously [11]. When considering stage 
(IV, a) it should be remembered that MnO, in aqueous solution is undoubtedly hydrated [7],and on the surface of UO, 
in acid solutions a hydroxide forms [12]. Consequently, the solid phase process, shown schematically in the form of 
equation (IV, a),in practice is much more complex, It is important that according to this scheme for the reaction 


between UO, and MnO, in acid solution, not containing foreign ions, it is essential to have contact between the solid 
phases, 


It is well known that UOg and MnO, in acid solutions dissolve rapidly; therefore all kinetic features of the total 
reaction, if the scheme (IV) is correct, are bound up with the reaction of UO, and MnO,; more accurately with processes 
occurring at the points of contact of the hydrated layers of these components. 

The observed dependence of the rate of the process on the size of the surface of UO, or MnO, will not be under- 


*The method for determining the rate of reaction was described in detail by the authors in an article submitted to the 
journal "Neorganicheskaya Khimiya." 
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standable if the reaction proceeded according to schemes (1), (II) or (III) on the whole surface of UO, and MnOg, since 
the rate of solution of solids in liquids is in a linear relationship to the size of the solid phase surface (for example see 
[13,14]). If the reaction occurs not on the whole surface of the UO, and MnO, particles, but only ina layer formed at 
the point of contact, it is obvious that the dependence of the rate of the process on the size of the surface of the react- 
ing components should not be linear, 


If the rate of the process is limited by the reaction (IV, a) then it should be expected that the concentration of 
hydrogen ions, not taking direct part in this stage of the process, should not affect the degree of reaction of UO, and 
MnO, in acid medium over a fixed length of time of the experiment. For a comparatively wide region of change in 
concentrations of sulfuric acid with an equimolar ratio of UO, and MnO, and 100 ml solution,this is confirmed by experi- 
mental data. The absence of an effect of sulfuric acid concentration on the degree of reaction of UO, and MnO, is 
difficult to reconcile with schemes (1), (II) or (III), since each of them includes the reaction of solution of oxides with 


the participation of hydrogen ions, these reactions depending to a large extent on the concentration of hydrogen ions 
[15,16]. 


The given experimental data agree most completely with scheme (IV). However, on their basis the possibility 
cannotbe excluded of the process occurring according to other schemes if the rate of the process can be effected in 
the same way by UO, and MnO, and the hydrogen ion concentration has no effect. 


Some Features of the Process 

We will consider at first the effect of preliminary grinding of UO, and MnO,. Samples of UO, (1.0 g) and MnO, 
(0.37 g) were ground for an hour in a mechanical SMB mortar, After this for 4 hr at 20°C in 50 ml of 0.5 N H,SO,, 
oxides were dissolved, separately or ground together and, in parallel, oxides were dissolved which had not been subjec- 
ted to additional grinding. The grain size of the initial UO, and MnO, was 0.074 mm, their degree of reaction was 
determined from the concentration of uranium in the solution after the end of the experiment. The data on the reac- 
tion of UO, and MnO, are as follows: 


Conditions for preparing speci- 
mens 

UO, (initial) 0.3 
UO, + MnO, (initial) 12 
UO,+ MnO, (ground separately) 45 
UO, + MnO, (ground together) 83 


These data show that the grinding of UO, and MnOg especially when they are ground together, has a very favor- 
able effect on the reaction of these substances in sulfuric acid solution, which confirms the correctness of scheme (IV). 
However, an x-ray study showed that the solid-phase reaction during the grinding of dry oxide for several hours only 
occurs to a very small extent —less than 0.1%, Consequently, the role of grinding is mainly to facilitate the formation 
of pairs of UO, and MnO, particles in close contact, 


Degree of reaction, % 


Other assumptions can also be put forward for the mechanism of reaction of UO, and MnO, in sulfuric acid med- 
ium. For example, we will asume that the reaction occurs at the moment of fast impact of free particles of these di- 
oxides and then these particles diverge. In this case the completeness of the reaction should depend on the number of 


impacts of particles and, consequently, should increase with decrease in the ratio S:L. The effect of this factor should 
therefore be studied. 


Experiments were carried out at room temperature in glass test tubes fastened on a disc rotating at a speed of 
60 rpm. The weight of UO, was 0.5 g, that of MnO, was 0.185 g. The oxides were ground together for an hour in a 
mortar. The H,SO, concentration was 5 N, time of the experiments ,30 min. When calculating the ratio S:L the weight 
of the solid was taken as the total weight of UO, and MnO». The degree of reaction of UO2 and MnO, in sulfuric acid 
solution in this case was determined by the concentration of uranium. The following data were obtained for the effect 
of the ratio S:L on the reaction of UO, and MnO,: 


S:L Degree of reaction 
1:1 22.6 
1:2 21.0 
1;3 21.5 
1:4 21.5 
1:5 23.8 
1;10 21.2 
1:15 21.2 
1:20 21.0 


The lack of dependence of the degree of reaction of components on the ratio S:L indicates that the reaction 
between UO, and MnO, in sulfuric acid solutions occurs not as a result of rapid reaction of the particles in brief con- 


tact with one another but probably occurs in the boundary layer which is strongly hydrated, forming at points of close 
and constant contact of UO», and MnO, particles. 


The relative strength of these pairs of particles is also confirmed by results of experiments in which a study was 
made of the effect of the intensity of mixing. On changing the speed of the four-blade mixer from 50 to 500 rpm the 


degree of reaction of UO, and MnOy, in 1 N or 5 N H,SO, remains practically unchanged both in the stoichiometric 
ratio of these oxides and with a large excess of MnOg, 


From the combination of all experimental data it can therefore be assumed that the reaction of UO, and MnO, 
in sulfuric acid solution proceeds according to the scheme (IV) in the absence of foreign ions, At first there is reaction 
in the surface hydrated layers of the UO, and MnO, particles in close contact; then the reaction products are rapidly 
dissolved,reaction again occurs on the fresh surface of the particles, etc. However, it must not be assumed that the 
process cannot occur according to schemes (1), (II) or (III) if at separate stages of successive reactions,unstable, rapidly 
decomposed substances are formed. In this case, as when the process occurs according to scheme (IV) if there is not 
sufficient contact between the reacting substances the rate of the process is limited by steric hindrances. 


The experimental data obtained in the present work throw more light on the role of iron ions during the sulfuric 


acid leaching of uranium from ores using pyrolusite. The assumption that iron is a direct oxidant and the process occurs 
according to the scheme 


+- 2F —> 2F efaa); (a) 
+ 4H faq) —> 2Fefaq) + Mniaq) 4. 2H,O; 1) (b) 


at first glance is in contradiction to the higher oxidation-reduction potential of Mn**/MnO, compared with Fe™ /Fe**, 


(V) 


Fejag te (E£°=0,77%); 


2Hy0¢1) + 4H faq) -+ 2e” = 4,22v). 


Further, if we consider the necessity for contact between the two solid phases for oxidation of UO, by means of 
MnO, in pure solution and the absence of necessity forsuch contact in the presence of iron ions in the solution, then 
it becomes understandable why the sulfuric acid leaching of primary uranium mineral from ores, which almost always 


contain iron, occurs according to scheme (V). 
SUMMARY 


1, A study has been made of the effect of the intensity of mixing, the ratio UOg:MnOy», the H,SO, concentration 
preliminary grinding of oxides and other factors on the reaction between UO, and MnO, in sulfuric acid solution. 


2. On the basis of the obtained experimental data, possible schemes have been discussed. 


3, It has been found that the occurrence of this reaction in which two solid phases take part is effected to a 
large extent by steric factors, In connection with this the problem was considered of the occurrence of the process 
with the presence in the solution of Fe* and Fe** ions, eliminating these steric hindrances. 


In conclusion we would like to thank V. G.Romanova for taking part in the work and L. V. Zverev for his useful 


comments. 
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A STUDY OF THE PROPERTIES OF URANIUM HEXAFLUORIDE 
IN ORGANIC SOLVENTS 


N. P. Galkin, B. N. Sudarikov, V. A. Zaitsev, 


D. A. Vilasov,and V. G. Kosarev 


Translated from Atomnaya Energiya, Vol. 10, No. 2, pp. 143-148, February, 1961 
Original article submitted April 14, 1960 


Studies have been made of the solubility and kinetics of solution of uranium hexafluoride in carbon 
tetrachloride, chloroform, dichloromethane, asymmetrical dichloroethane,symmetrical tetrachloro- 
ethane, pentachloroethane, trifluorotrichloroethane, symmetrical trichloropropane,and tetrachloro- 
propane. It has been shown that solutions of uranium hexafluoride in carbon tetrachloride, tetrachloro- 
ethane, pentachloroethane,and trifluorotrichloroethane are completely stable for two weeks at 20°C; 
at this temperature solutions of uranium hexafluoride in chloroform, dichloroethane ,and dichloro- 
methane are unstable. It has been shown that reactions of uranium hexafluoride with these studied 
organic solvents at 60-100°C occur in the following way: at first uranium pentafluoride is formed 
which is reduced at first to the intermediate uranium fluorides, containing a large amount of tetra- 
valent uranium,and then to uranium tetra fluoride. 


It is known that uranium hexafluoride is soluble in carbon tetrachloride, chloroform, symmetrical tetrachloro- 
ethane and other halogenated hydrocarbons [1,2]. In [3] the solubility of uranium hexafluoride was studied in carbon 
tetrachloride in the temperature range 5-35°C, 


Solutions of uranium hexafluoride and chlorinated and fluorinated hydrocarbons have varying stability. For 
example, uranium hexafluoride forms a stable solution in symmetrical tetrachloroethane. The yellow solution loses 
its color on boiling; on cooling, the yellow color is restored. This change in color is explained by the formation of 
complex compounds of uranium hexafluoride with tetrachloroethane which becomes colorless on boiling. Solutions 


of uranium hexafluoride in pentachloroethane have a stability which is still greater than in symmetrical tetrachloro- 
ethane. 


In addition to this uranium hexafluoride reacts with some organic solvents, For example, 1,2-difluoro-1,1,2,2- 
tetrachloroethane reacts fairly rapidly with uranium hexafluoride. Gases are liberated and as a result the whole of the 
uranium hexafluoride is reduced to uranium tetra fluoride [1]. Some organic solvents can be used to reduce uranium 
hexafluoride to the tetrafluoride at fairly low temperatures, The reduction of uranium hexafluoride by liquid trichloro- 
ethylene at ~80°C was described in [4]. 


On heating solutions of uranium hexafluoride in carbon tetrachloride up to temperatures of 150°C in an auto- 
clave a reaction takes place, accompanied by an increase in pressure and temperature. The reaction products are 


uranium tetrafluoride, chlorine and a mixture of freons CClsF and CFCls. The equation of the reaction will have the 
following form: 


UF, + 2CCl, UF,+Cl, + 2CCIF. 


The aim of the present work is to determine the solubility of uranium hexafluoride in halogenated hydrocarbons 
which are the most stable to uranium hexafluoride, and to study the stability of the obtained solutions. Most of the 
selected solvents are readily available and cheap. 


The solubility of uranium hexafluoride was determined in a quartz vessel with a seal. The sealing liquid was a 
mixture of completely fluorinated hydrocarbons with a boiling point of ~140°C, A fixed volume of the organic solvent 
was added to the vessel and solid uranium hexafluoride in an amount needed to form a saturated solution at the parti- 
cular temperature. When mixing was complete and the liquid had been allowed to stand, a fixed volume of solution 
was removed for analysis. The uranium content in the organic phase was determined volumetrically. All the organic 
solvents used in the work were subjected to careful chemical purification and repeated fractional distillation. 
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Solubility of UF,, g/ml 
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Fig. 1. The effect of time of mixing on the 
solubility of uranium hexafluoride at 25°C 

in halogenated hydrocarbons: 1) CHgCl,; 2) 
CCl gF3; 3) CHCls; 4) CCly. 


The kinetics of solution of uranium hexafluoride in all organic 
solvents was studied at 25°C; experiments with dichloroethane were 

carried out at 10°C. Figures 1 and 2 give the solubilities of uranium 
hexafluoride in halogenated hydrocarbons for different times of mix- 
ing. 


As can be seen from these data, equilibrium in the solution of 
uranium hexafluoride in most organic solvents is established in 1 hr. 
Only in chloroform and trichloropropane is the equilibrium established 
in 3 hr. The values are given below for the solubility of uranium hexa- 
fluoride in the investigated halogenated hydrocarbons at 25°C; 


Organic solvent Solubility of UF,, g/ml 


Solubility of UF,, g/ml 


CH,Cl, 0.983 
CHCl, 1.095 
CCl, 0.799 
C,H,Cl, 0.626 
0,520 
C,HCl, 0.490 
0.983 
C3H;Cly 0.415 


C3H,Cly 


0.600 


05} 
4 
041g 
08 

0,2 


80 60 90 120 150 180 210 240 
Time of contact, min 


Solutions of uranium hexafluoride in the studied halogenated hydrocarbons had different colors; green—in 
carbon tetrachloride, orange—inchloroform, yellow—in tetrachloroethane and pentachloroethane, cherry red~in di- 
chloroethane, dichloromethane, tetrachloropropane and trichloropropane. An exception is the colorless solution of 

uranium hexafluoride in trifluorotrichloroethane. The intensity of the color of the solutions increases with the con- 
centration of uranium hexafluoride. The formation of colored solutions is probably due to the formation of complex 
compounds of uranium hexafluoride with the organic solvents. 


When studying the dependence of uranium hexafluoride solu- 
bility on temperature the time of mixing was 1 hr for all organic sol- 
vents; only for chloroform and trichloropropane was the time of mix- 
ing 3 hr. Figures 3 and 4 give the solubilities of uranium hexafluoride 
in halogenated hydrocarbons at different temperatures. In the investi- 
gated temperature range the solubility of uranium hexafluoride in all 
halogenated hydrocarbons increased with temperature. Chlorine deri- 
vatives of methane had greater dissolving capacity than the chlorine 

derivatives of ethane and propane, 


A study was also made of uranium tetrafluoride solubility and 
that of uranyl fluoride in carbon tetrachloride and chloroform at 25°C, 
These studies were of interest because these compounds can be present 
in uranium hexafluoride in small amounts, A determination was also 
made of the solubility of the double salt 3NaF*UF, at 25°C. The results 


Fig. 2. Effect of time of mixing on the solubility of the experiments are given in Table 1. 
of uranium hexafluoride at 25°C in halogenated 
hydrocarbons; 1) CgH4Cl4; 2) CpzH2Cly; 3) CgHClg; 
4) CsHsCly; 5) (at 10°C). 


As can be seen from the data, the solubility of uranium tetra- 
fluoride and uranyl fluoride in carbon tetrachloride and chloroform 
is very low (it is presumably of the same order in the other investi- 
gated solvents), the solubility of the double salt 3NaF-UF, is some- 


what higher, which is possibly due to the presence in it of hexafluoride which was partially absorbed during the pre- 
paration of the salt. 
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TABLE 1. The Solubility of Uranium Tetrafluoride, 


Uranyl Fluoride and the Double Salt 3NaF-UF, in 
Carbon Tetrachloride and Chloroform at 25°C 
| Solubility mg/liter U 
Ga 44 | | incucis 
| | UF, <15 15 
= eae 3NaF- UF, 72 
q 0 +— 
48-—- | 
ig he 4 TABLE 2. The Stability of Solutions of Uranium Hexa- 
yze fluoride in Organic Solvents as a Function of the Stand- 
, ing Time at 20°C 
20 40 50 


1,6 


Temperature, °C Organic Initial con- 
i 
Fig. 3. The temperature dependence of the solubility solvent centration prem: aoe vad 
of uranium hexafluoride in halogenated hydrocarbons; of UF*, g/ml | 
1) CHClg; 2) 3) CH,Cl,; 4) CCly. ing, & g/m 
CCl {0,273 0,270 0,272 
(0, 185 0,186 0,184 


0,199 0,188 0,189 
(symmetri- | {97130 0,127 0,128 
1,2 cal) 
10 27 
0,304 0,300 0,302 
> 98 CoClyFs {0206 0,203 0,204 
‘a (= TABLE 3. Degree of Reduction of Uranium Hexafluoride 
during 30 Days at 20°C 
+- = or | Organic Initial con- Content of (Degree of ‘ 
tent of US 4 reduction 
solvent 
0,0 70 2 30 70 30 g uranium, ‘/o 
Temperature, °C 
Fig. 4. The temperature dependence of the solubility CCl, 0,920 0,020 pai 
of uranium hexafluoride in halogenated hydrocarbons: CoHQCly 0,933 0,020 > a | 
1) 2) 3) 4) 5) 
CaHsCls. 0,995 0,021 2,0 
0,630 0,019 3,2 


The method for studying the stability of solutions at room temperature (20°C) was as follows. Solutions of uran- 
ium hexafluoride of known concentration in aluminum containers were placed in a desiccator, After certain intervals 
of time the contents of hexa- and tetravalent uranium in the solutions were determined volumetrically. No tetravalent 
uranium was found in any of the cases. The results of the experiments are given in Table 2. 


The constancy of the uranium hexafluoride concentration, the absence of tetravalent uranium in the solution, 
the stability of the color of the solution and the absence of precipitates in them indicated that solutions of uranium 
hexafluoride in the studied organic solvents are completely stable at 20°C for two weeks. 


A study was made of the stability of solutions of uranium hexafluoride in carbon tetrachloride, symmetrical 
tetrachloroethane, pentachloroethane and trifluorotrichloroethane with a standing time of about 30 days. The method 
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in Table 4. 


A study was made of the stability of solutions of uranium hexafluoride in chloroform and asymmetrical dichloro- 
ethane at 20°C and a standing time of 7 days. In both cases precipitates were formed; the volume of the solution was 
5 ml. The studies were carried out according to the above-described method, The results of the experiments are given 


TABLE 4. Degree of Reduction of Uranium Hexafluoride 
during 30 Days at 20°C 


| 

Organic _ Initial con- Content ust | Degree of 
solvent tent Us*, g g reduction of 
uranium, % 

CHCl, 0,603 0,162 27 

CyH,Cly 0,500 0,150 30 

(asymmet- 
rical) 


TABLE 5, Temperature Dependence of the Degree of 
Reduction of Uranium Hexafluoride by Dichloromethane 


Temperature | Content of U*" in | Total Content 
solution, g of 
10 0,021 0,048 
20 0,012 0,059 
25 0,009 0) ,066 
| 


TABLE 6. Degree of Reduction of Uranium Hexafluo- 
ride by Organic Solvents (in Percentages) 


Organic Temperature, °C 

solvent 60 | 80 | 100 
CCl, 0 0 57,5 
CHCl, 1,1 6 56,0 
CHCl, 10,8 27,7 57,1 
CoHyCl, 20,5 36,8 50,7 
3,9 5,7 18,5 
CHCl, 2,4 2A 2,4 
C3H,Cls 2,9 29,9 57,7 
47,5 70,8 82,7 


Under other conditions white products were obtained with a greenish hue of varying intensity (Table 7). These 
products were analyzed in the following way: the precipitates were separated from the organic phase on filter paper, 
washed with carbon tetrachloride and the content of hexa- and tetravalent uranium in them was determined. 


of investigation was as before,but in connection with the formation of precipitates the content of tetra- and hexavalent 
uranium was determined in the whole contents of the containers. In all cases the volume of the solutions was 5 ml. The 
results of the experiments are given in Table 3. 


It can be seen from these data that regardless of the initial concentration of the uranium hexafluoride and the 
type of organic solvents, the tetravalent uranium content remains constant; this is probably due to the reaction of 
uranium hexafluoride with the material of the container. 


It can be seen from these data that solution of uran- 
ium hexafluoride in chloroform and asymmetrical dichloro- 
methane are less stable than in the previously studied 
organic solvents, Solutions of uranium hexafluoride in di- 
chloromethane were also unstable. Partial reaction was 
observed as soon as the solid uranium hexafluoride was 
mixed with the solvent. The conditions of the experiments 
were as follows: the volume of solvent was 3 ml, UF,— 

8 g; time of mixing was 30 min. The results of the experi- 
ments are given in Table 5. 


The stability of solutions of uranium hexafluoride in 
some chlorinated derivatives of methane, ethane and pro- 
pane at 60, 80 and 100°C were studied in the following 
way: solutions of uranium hexafluoride with the same con- 
centration were placed in a steel bomb which was kept in 
a thermostat for 1 hr at a given temperature. After the 
completion of the experiment the total content of hexa- 
and tetravalent uranium were determined volumetrically 
in the obtained products, 


In experiments in which the degree of reduction of 
uranium did not exceed 2.5% (for CHClg at 60°C, for 
C,HCls at 60, 80 and 100°C, for CsHsClg at 60°C), the 
reaction products were a homogeneous organic phase, In 
view of the very low solubility of uranium tetrafluoride in 
organic solvents, the presence in the solution of tetravelent 
uranium was probably due to the solubility of the intermedi- 
ate fluorides of uranium, formed during the reaction of uran- 
ium hexafluoride with organic solvents, In the other experi- 
ments precipitates were formed. Due to the reaction of uran- 
ium hexafluoride with chloroform and trichloropropane (at 
80°C), dichloromethane and dichloroethane (at 60° and 80°C), 
tetrachloroethane (at 60, 80 and 100°C) and tetrachloropro- 
pane at (60°C) white compounds were obtained which dis- 
proportionated in air in 15-30 min and were almost instant- 
aneously hydrolyzed in water with the formation of uranium 
tetrafluoride and uranyl fluoride. In the experiments the 
volumes of the solutions were 3 ml and the concentration 
of UF, was 0,48 g/ml. The results of the experiments are 
given in Table 6, 


TABLE 7. The Relationship of Hexa- and Tetravalent The investigations show that the reaction of uranium 

Uranium in the Obtained Products hexafluoride with various organic solvents occurs in the fol- 
lowing way: initially uranium pentafluoride is formed which 
is first reduced to the intermediate fluoride of uranium con- 


Organic Ratio of to Degree of re- 
solvent duction 


uranium, taining a large amount of tetravalent uranium,and then to 
CHCl, 56,0 uranium tetra fluoride. 


CH,CI, 36 57,1 The following conclusions can be drawn from the study 
CHCl ; 57,7 of the solubility of uranium hexafluoride in some halogen 
Csl,Cl, \7 82,7 derivatives of methane, ethane and propane and a study of the 


1, Chlorine derivatives of methane have a greater dissolving capacity than those of ethane and propane. 
2. The solubility of uranium tetrafluoride and uranyl fluoride in carbon tetrachloride and chloroform is very low. 


3. At 20°C solutions of uranium hexafluoride in carbon tetrachloride, tetrachloroethane, pentachloroethane and 
trifluorotrichloroethane are stable,and in chloroform, dichloromethane and dichloroethane they are unstable. 


4, In the reaction of uranium hexafluoride with the studied organic solvents uranium pentafluoride is formed and 
is first reduced to intermediate fluorides of uranium and then to uranium tetrafluoride. 
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METHODS OF REDUCING URANIUM HEXAFLUORIDE 


N. P. Galkin, B. N. Sudarikov,and V. A. Zaitsev 


Translated from Atomnaya Energiya, Vol. 10, No. 2, pp. 149-155, February, 1961 
Original article submitted February 8, 1960 


On the basis of published material the present article gives brief characteristics of different methods 
for reducing smelting. Values are given for the free energy of some investigated reactions of the reduc - 
tion of uranium hexafluoride. The effect of various factors is considered on the purity and the granular 
specific gravity of the obtained uranium tetrafluoride. The advantages and disadvantages of different 
schemes for reducing uranium hexafluoride are discussed. 


Uranium hexafluoride is the most important product in the nuclear engineering industry since it is essential in 
the separation of uranium isotopes by the gaseous diffusion method. 


Many papers have appeared recently on the fluoride method of treating irradiated nuclear fuel (the fluoride - 
vaporization process) with which the uranium in the form of uranium hexafluoride is separated from the active fission 
products [1-5]. Papers have also been published on the preparation of nuclear-pure uranium hexafluoride from uranium 
ore concentrates and rich uranium ores by the method of direct fluorination [6,7]. 


Uranium hexafluoride obtained by a given method should later be processed to metallic uranium or some of its 
compounds (uranium dioxide, uranium tetrafluoride, etc.), which have been used as nuclear fuel. In the conversion of 
uranium hexafluoride to the dioxide the most convenient methods will presumably be those based on the hydrolysis of 


uranium hexafluoride by water or by aqueous solutions with subsequent precipitation and roasting of the ammonium 
diuranate [8]. 


To prepare metallic uranium the chemical engineering conversion of uranium hexafluoride is carried out in two 
stages: at first the uranium hexafluoride is converted to uranium tetrafluoride, then the latter is subjected to reducing 
smelting. The reduction of uranium hexafluoride is therefore of considerable practical importance. 


Various authors have studied the reduction of uranium hexafluoride to the tetrafluoride by hydrogen, hydrogen 
chloride, carbon tetrachloride, ammonia, ethylene, propane, trichloroethylene, sulfur dioxide, silicon tetrachloride, 
thionyl chloride and by other reducing agents, Some of the work was developed to pilot-plant and industrial scale, The 
most promising reducing agents for uranium hexafluoride from the point of view of obtaining uranium tetra fluoride 
suitable for reducing smelting are presumably hydrogen, carbon tetrachloride, hydrogen chloride, tichloroethylene 
and ammonia. Table 1 give the values of free energy for some of the investigated reduction reactions. 


As can be seen from these data, in the series HCl, CCl4, HBr, Hg there is a reduction in the value of free energy 
for the reduction of uranium hexafluoride. In this direction the thermal effect of the reaction increases and the chemi- 
cal affinity of the reacting substances increases. In connection with the large thermal effect (83.6 + 2.5 kcal/mole 
UF, at a temperature of -30°C) it should be expected that the value of free energy in the reaction of uranium hexa- 
fluoride with ammonia will be lower than for the case of reduction of uranium hexafluoride by hydrogen. 


Reduction of Uranium Hexafluoride by Hydrogen Chloride 

In [9] there was a description of the reduction of uranium hexafluoride by hydrogen chloride; it was stated that 
at 200°C a mixture is obtained of brown, green and black products which become green on being degassed in a vacuum, 
The optimum temperature range of the reduction was 250-3800°C, 


An improved method for reducing uranium hexafluoride by hydrogen chloride was described in [10]. The reac- 
tion was carried out at 200—400°C in a vessel having the shape of an inverted cone. Vapors of uranium hexafluoride 
and hydrogen chloride were introduced tangentially into the chamber near its bottom where the stream of gases was 
relatively fast and had high turbulence. On reaching the upper part of the chamber where the area of cross section 
increased, the stream lost its speed and its turbulence was reduced, which helped the precipitation of uranium tetra- 


| 


TABLE 1. The Values of Free Energy for the Reduction 
of Uranium Hexa fluoride 


Reaction 


Free energy, 
kcal/mole 


UFas) UF 2HF @) 


—68,9 


UF es) + 2HBrg) —> 
UFas) + 2H Fig) + Brag) 


—43,5 


UFe(s) +-2CCl a1) 
—» UF ds) 2CC15F 
with a molar ratio 
CC), UFs=4: 1 


UF es) +- 2HClg) > 
— UF + 2HKg)+ Chg) 


fluoride particles. The particle sizes of the product could be 
controlled by the volumetric rate of feed of the reagents, A 
6-8-fold excess of hydrogen chloride in the mixture was used 
with different diluents, The granular weight of the uranium 
tetrafluoride varied within the limits of 1.6—2.0 g/cm’, The 
authors of [10] state that the uranium tetrafluoride obtained 
in this way is suitable for direct reducing smelting in order 
to produce uranium, 


Reduction of Uranium Hexafluoride by Hydro- 
gen 

Many papers have been devoted to the reduction of 
uranium hexafluoride by hydrogen. In [11] an account was 
given of the practically instantaneous reduction of uranium 
hexafluoride to the tetrafluoride by hydrogen at room tempera- 
ture, however the experiments carried out in 1943 by Spencer- 
Palmer [12] did not confirm these data. At temperatures of 
78 —100°C uranium hexafluoride is practically unreactive 
with hydrogen, It is probable that the materials used in [12] 
contained impurities which were catalysts for this reaction. 


The reduction of uranium hexafluoride by hydrogen 


was studied in greater detail in [12]. A study was made of the dependence of the degree of reduction of uranium hexa- 
fluoride by hydrogen on the temperature (Table 2) and the time of reaction (Fig. 1). 


was reduced by hydrogen according to the reaction UF¢g) + Hxg)> UF 4s) + 2HF, 


TABLE 2. Effect of Temperature on the Degree of Reduction of Uranium Hexa fluoride by Hydrogen 


Temperature of reaction, °C 


226 


250 


275 | 285 | 292 305 


Degree of reduction of uranium, % 


3,8 


22,3 


40,1 50,5 


42,5 | 44,3 


Note. Time of process 30 min, molar ratio Hg; UF, = 5:1. 


At a temperature of 121°C no reaction was observed between uranium hexafluoride and hydrogen; the reaction 
started in the temperature range 215—288°C, These data agree well with the results of the Spencer-Palmer work. 
Based on preliminary data, the authors observe that the reduction of uranium hexafluoride by hydrogen is a first order 
reaction and proceeds both in the gaseous phase and at the walls of the reactor. High activation energy is needed for 
it. This fact is mentioned in a number of other papers. It is mentioned that the reduction of uranium hexafluoride by 
hydrogen is accelerated catalytically in the presence of certain chlorides (hydrogen chloride, mercury monochloride, 
ferric chloride); however, in these cases also the reaction is comparatively slow [13]. 


At the Second Geneva Conference on the Peaceful Uses of Atomic Energy the US delegation presented a report 
[14] on the successful reduction of uranium hexafluoride by hydrogen on an industrial scale. Uranium hexafluoride 


+68.9 kcal by two methods, In the 


first of them (the method of reactor with “hot™ walls) the heat needed for the start of the reaction was communicated 
to the system through the wall of the reactor, in the second (the method of reactor with “cold” walls) the required 


heat was obtained by the reaction of hydrogen with elementary fluorine which was added together with the uranium 
hexafluoride. 


The layout of the reactor with hot walls is shown in Fig. 2. The reactor (150-mm, with a tube of monel metal) 
was split up into five temperature zones; the wall temperature of the reactor varied within the limits 426—537°C. The 
feed rate of the uranium hexafluoride was 45,36 kg/hr, a 100% excess of hydrogen was used above the stoichiometric 
amount, The mean granular specific gravity of the obtained uranium tetrafluoride was 3.6 g/cm®, The total of metal- 
lic impurities in the uranium tetrafluoride varied from 2.0°10~* to 1.3-10°%, averaging 4.0-10°%%, The faults in this 
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Fig. 1. Effect of Reaction time on the degree 
of reduction of uranium hexafluoride by hydro- 9 7 8 
gen (at temperature 249—250°C). Molar ratio 


Fig. 2. Layout of reactor with hot walls: 1) 
nozzle; 2) reactor housing; 3) pneumatic 
vibrator; 4) Cottrell filter; 5) heater; 6) 
mechanical vibrator; 7) unloading worm; 8) 
motor; 9) receiver for UF,. 


method for reducing uranium hexafluoride, as shown by the authors themselves [14], are the considerable melting of 
the product on the reactor walls, especially in the upper zone and the formation of hardened material on the walls. 
After every 8 hr of operation the apparatus was stopped and the layer of fused uranium tetrafluoride was removed 
mechanically. Vibrating the reactor walls did not remove the fused product. 


In the reduction of uranium hexafluoride in a reactor with cold walls an attempt was made to preheat the uran- 
ium hexafluoride and hydrogen to the required temperature; however, it was found that at 426°C the uranium hexa- 
fluoride reacts with the material of the tube (monel metal) forming solid crystalline products, presumably a mixture 
of intermediate uranium fluorides. The result was that the tubes choked up and the apparatus became unserviceable, 
Elementary fluorine was therefore added to the uranium hexafluoride to increase the temperature in the reaction zone. 
From the heat liberated during the reaction F,+ Hy > 2HF + 128 kcal, the starting materials were heated to the required 
temperature. The reduction of uranium hexafluoride by hydrogen in this case proceeded entirely within the flame with 


TABLE 3, Effect of the Gravimetric Ratio on the Granular Weight of UF, 


Gravimetric ratio F,:UF, 0,040 | 0,045 | 0,050} 0,055 0,060 | 0,00 0,070 0,07 0,075} 0,075 


Granular density of UF,, g/cm*® 


2,18 | 2,59 | 2,68 | 2,99 | 2,89 | 3,08 | 2,99 | 3,08 | 3,58 | 3,71 


the reactor walls at temperatures of 92-204°C. The layout of the apparatus is shown in Fig. 3. The gas rate at the 
nozzle varied from 15 to 40 m/sec with a rate of production of uranium hexafluoride of 4.53-11,34 kg/hr, The excess 
of hydrogen above the stoichiometric amount was 100-500%; the stream of fluorine was 0.05-0.08 kg F,/UFs. 


For low flow coefficients with regard to fluorine the reduction of uranium hexafluoride was incomplete and the 
final product was contaminated with intermediate fluoride, The granular weight of the obtained uranium tetra fluoride 


was also determined by the gravimetric ratio of fluorine and uranium hexafluoride (Table 3). The thickness of the 
uranium tetrafluoride layer on the walls did not exceed 0.8 mm. 


10 


16 
Fig. 3. Apparatus for reducing uranium hexafluoride Fig. 4. Apparatus for reducing uranium hexafluoride 
by hydrogen using a reactor with cold walls; 1) H, with trichloroethylene; 1) UF, cylinders; 2) receivers; 
cylinders; 2) receiver; 3) diaphragm; 4) control 3) mixers; 4) heaters; 5) Cottrell filter; 6) pneumatic 
valve; 5) reactor; 6) porous filters; 7) traps; 8) un- vibrator; 7) reactor; 8) water-cooled condenser; 9) 
. loading worm; 9) UF, receiver; 10) UF, cylinders; ice-cooled condenser; 10,11) receiver for freons; 12) 
: 11) container for UF. porous filter; 13) unloading worm; 14) receiver for 


UF,; 15) mechanical vibrator; 16) cylinders with CpHCls. 


With a certain modification the method of the reactor with cold walls was described in [15], where the energy 
needed for the start of the reaction was obtained by the recombination of atomic hydrogen: 


2H — H, +103 keal. 


Reduction of Uranium Hexafluoride by Ammonia 


It has been stated that uranium hexafluoride reacts with ammonia at a temperature of -72°C, forming a mix- 
ture of solid products. 


At 300°C gaseous uranium hexafluoride reacts with ammonia in a reactor made of nickel and forms a mixture 
of ammonium uranium pentafluoride and ammonium fluoride (98% NH,UF,+ 5% NH,F). X-ray structural analysis shows 
that there is no uranium tetrafluoride in this product. A product of this composition was obtained on fusing uranium 
tetrafluoride with ammonium bifluoride at 450°C [16]. The reaction of uranium hexafluoride with ammonia at 450°C 
also leads to the formation of ammonium uranium pentafluoride [15]. 


More detailed information on the reaction of uranium hexafluoride with ammonia is given in [17], where a 
study was made of the composition of products from the reaction of uranium hexafluoride with ammonia in the temper- 
ature range -50 to 200°C and the total equations are proposed for the reaction: for temperatures of -50 to -30°C 


6UF, + (8 +6n)NH, —> 6UF,-nNH, + 6NH,F + Ny, 


where n= 0,73; 
for temperatures of 0-25°C 


4UF,+-8NH, 2UF; +} 2NH,UF, 4NH,F + N,; 
for temperatures of 100-200°C 


3UF, + 8NH, 3NH,UF, 3NH,F 


The rate of reaction of solid uranium hexafluoride with gaseous ammonia was examined in the temperature 
range -20 to 20°C, The reaction rate decreased with decrease in temperature; however, even at -20°C the reaction 

was finished in 3-5 min, The reaction of uranium hexafluoride with ammonia occurred with a large evolution of heat 
even at low temperatures (Table 4). At 500°C ammonium uranium pentafluoride decomposed into uranium tetrafluo- 
ride and ammonium fluoride. The nitrogen content in the final product did not exceed 10° %o [18-20]. 
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TABLE 4. The Thermal Effect of the Reaction of Solid 
Ura nium Hexafluoride with Liquid Ammonia 


Temperature of reac- | Thermal effect of reaction 


tion, °C kcal/mole UF, 
—50 50,8+1,5 
—40 67,0-+2,0 
30 83,6+2,5 


S 


Solubility of UF, in CCly, g/ml 


5 100 1 20 25 30 39 
Temperature, °C 


Fig. 5. Effect of temperature on the solubility of 
uranium hexafluoride in carbon tetrachloride. 
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Fig. 6. Apparatus for reducing uranium hexafluoride by 
carbon tetrachloride: 1) CCl, cylinder; 2) rotameter; 3) 
vaporizer; 4) line of preliminary heating of CCl,; 5) line 
of preliminary heating of UF,; 6) rotary furnace; 7) 
manometer; 8) Kronberger flowmeter; 9) diaphragm; 

10) differential manometer; 11) UF, cylinder; 12) cycl- 
ones; 13) condenser; 14) condensate receiver; 15) tank 
containing NaOH solution; 16) scrubber; 17) receiver 

for solution. 
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Reduction of Uranium Hexafluoride 
by Trichloroethylene 

Trichloroethylene in the liquid and gas phases re- 
acts with uranium hexafluoride according to [14,21] 


UF, -+CHCI=CCl, —» UF,+HF + 


-+ mixture of chlorofluorohydrocarbons 


The mixture of chlorofluorohydrocarbons is formed 
due to the reactions of addition, substitution, and poly- 
merization and also to splitting of trichloroethylene 
molecules and subsequent formation of chlorofluoro deri- 
vatives of methane. 


The reduction was carried out at 82°C in a spray- 
ing column and the uranium hexafluoride was first dilu- 
ted by nitrogen. The resulting uranium tetrafluoride was 
collected as a slurry at the bottom of the column and 
separated from the excess trichloroethylene by a centri- 
fuge. In [14] it was mentioned that during operation the 
gas nozzles became blocked and the separation of the 
solid and liquid phases was seriously hindered. 


In the gas phase the reaction was carried out at 
232°C with a molar ratio CyHClg:UF, = 3:1 and the initial 
concentration of uranium hexaflouride was 0.5-3.5 
mole %, The rate of flow of the uranium hexaflouride 
varied within the limits 0.9-2.3 kg/hr. 


Under these conditions the uranium hexafluoride 
was completely reduced to uranium tetrafluoride with 
a granular weight of 0.22-1.6 g/cm’, The granular weigh 
of uranium tetrafluoride was increased to 2.96 g/cm® 
on roasting its powder in an atmosphere of hydrogen fluo- 
ride at 900°C for 2 hr (Fig. 4). 


The obtained uranium tetrafluoride was contamina- 
ted to a large extent with carbon,and to remove this the 
product was roasted in an atmosphere of nitrogen, hydro- 
gen and oxygen. Only treatment of uranium tetrafluoride 
with oxygen led to a considerable reduction in the carbon 
content, however in this case to prevent the formation 


of uranyl fluoride it was essential to add hydrogen fluo- 
ride, 


The described process was carried out on a pilot- 
plant scale, The gravimetric ratio of the trichloroethyl- 
ene and uranium hexafluoride varied within the limits 
0.7-3,25, The temperature of the gases at the input to 
the reactor varied within the range 82-304°C; the reac- 
tor temperature was 74-483°C. The density of the obtained 
uranium tetrafluoride did not exceed 1.6 g/cm’, As 
regards the carbon content the product was below stand- 
ard in all cases, 


A serious fault of this method reducing uranium 
hexafluoride (apart from the low granular weight of 


uranium tetrafluoride and the high content of carbon in it) were the complications connected with the blocking of the 


tubes with readily condensing organic reaction products, especially hexachlorobenzene. 


= 
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Reduction of Uranium Hexafluoride by Carbon Tetrachloride 


It is known that uranium hexafluoride is readily soluble in carbon tetrachloride; the solution is stable at room 


temperature for several weeks (see page 356 of [13]). Figure 5 shows the solubility curve for uranium hexafluoride in 
carbon tetrachloride as a function of the temperature, 


When a solution of uranium hexafluoride in carbon tetrachloride was heated in an autoclave to 150°C or when 
the vapors of uranium hexafluoride and carbon tetrachloride were reacted at 200°C and above,the uranium hexafluo- 
ride was reduced to the tetrafluoride according to the equation [22]. 


UF, 2CCl, —> UF, + Cl, + 2CCI,F. 


It has been mentioned that in both cases this reaction proceeded stepwise: as well as the uranium tetrafluoride and 


elementary chlorine there were mixed chlorofluoro derivatives of methane, which underwent further change during 
the reaction. 


Traces of uranium tetrafluoride had a catalytic effect on the reduction of uranium hexafluoride by carbon tetra- 
chloride. 


In the condensed phase the reaction was carried out in the following way; 1 liter of CCl, was added per 1 kg of 
UF, in an evacuated autoclave; the autoclave was then heated from outside to 150°C and further to 200°C from the 
heat of reaction. In 2 hr the pressure in the apparatus was raised to 25 atm, then the autoclave was cooled to 50°C, the 
gases were passed through alkaline scrubbers, where the chlorine was trapped, and then the gases were put into the 
atmosphere. There was 99,7% reduction of uranium hexafluoride to the tetrafluoride; the obtained uranium tetrafluo- 
ride was heated in a current of air at 120°C to remove chlorine and other adsorbed gases. When air was drawn in and 


when the vessel was poorly evacuated the reduction of uranium hexafluoride was accompanied by the formation of 
phosgene, 


Uranium hexafluoride was quantitatively reduced by carbon tetrachloride with a molar ratio of components of 
1:4 or with still greater excesses of carbon tetrachloride. A reduction in the excess of carbon tetrachloride led to an 
increase in the content of intermediate uranium fluorides in the reaction products, 


In the gaseous phase there was continuous reduction of uranium hexafluoride by carbon tetrachloride at atmo- 
spheric pressure and temperatures of 200-400°C in a horizontal furnace of diameter 108 mm and length 1676.4 mm. 
The layout of the apparatus is shown in Fig. 6. 


At 200°C for the quantitative reaction of uranium hexafluoride with carbon tetrachloride with the formation of 
solid nonvolatile products (intermediate uranium fluorides) the time needed is 2 min; the reduction of these products 
in a stream of carbon tetrachloride to uranium tetrafluoride is accomplished in 10 min. 


Figure 7 shows the dependence of the degree of reduction 
of uranium hexafluoride by carbon tetrachloride on the reac- 


tion temperature, the rate of production and the ratio of the 

reagents. 
: %, CON be By reducing uranium hexafluoride with carbon tetrachlo 
o 310} % o% Op ride uranium tetrafluoride is obtained with a high granular 
: 290. *x NT weight and a content of impurities within the limits of the 
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1 1 a W. Page et al.,Report BNL-174 (March 1952), 
eft 2. G. Cather and R. Lenze, Nuclear Engineering and Science 
Congress, Ohio (December 1955). 
Fig. 7. The effect of temperature, rate of produc - 3. R. Steunenberg et al,, Report No, 539 presented by the 


Molar ration CC1,/UF, 


tion and ratio of reagents on the degree of reduc- USA to the Second International Conference on the Peace- 
tion of uranium hexafluoride by carbon tetrachlo- ful Uses of Atomic Energy (Geneva, 1958), 
ride, Rate of production; 1) 1 kg/hr uranium; 2) 4, | Steunenberg and Vogel, Transactions of the Second Inter- 


0.6 kg/hr uranium. 
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LETTERS TO THE EDITOR 


THE MECHANISM OF REACTION OF FAST NUCLEONS 
WITH NUCLEI 


V. S. Barashenkov, V. M. Mal'’tsev,and K. Mikhul 


Translated from Atomnaya Energiya, Vol, 10, No, 2, pp. 156-158, February, 1961 
Original article submitted July 14, 1960 


Several papers have been published on the experimental investigation of the mechanism of reaction of protons 
with energy of 9 Bev with nuclei (see bibliography in [1-5]). However, there is no single opinion on this problem: 
some are of the opinion that at energies of 9-10 Bev the mechanism of an intranuclear cascade, realized at lower 
energies [6], is replaced by a mechanism of simultaneous reaction with the whole “tube” of the nuclear material 


[4]; others consider that at an energy of 9 Bev the reaction of fast nucleons with nuclei occurs by means of intra- 
nuclear cascades [1,2]. 


In order to find which of the two discussed mechanisms is realized in the reaction of a nucleon with 9 Bev energy 
with a nucleus, we carried out calculations of an intranuclearcascade and compared them with experimental data. 

The calculations were made by the Monte Carlo method for the case of relativistic three-dimensional kinematics, 
Allowance was made for the multiple production of particles due to collision between fast particles and nucleons of 
the nucleus. The angular and energy distributions of the particles after elastic and inelastic (NN)- and (™N)-collisions 
were taken from experiments. A Fermi gas nuclear model was considered with the same parameters as those given in 
[6]. 


The results of the calculations and the known experimental data on shower and cascade particles* are given 
in Tables 1 and 2 and in Figs. 1-3 (the dotted line shows the experimental data of [1]). 
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Fig. 1. Angular distribution of shower particles: a) group of light nuclei; b) group of 
heavy nuclei. 


* We apply the term shower particles to nucleons with energies of 0.5-9 Bev and ™ -mesons with energies of 0,08-8 
Bev. We use the term cascade particles for nucleons with energies of 0,03-0.05 Bev and ™ -mesons with energies of 
0.15-0,08 Bev (for more details see [1]). All particles with lower energies will be referred to as evaporation particles. 
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TABLE 1. Mean Number of Particles Produced in One 
Act of Inelastic Nucleon-Nuclear Collision 


Relative number of particles, % 
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Relative number of particles, % 


Group of light Group of heavy 

Particles 
| calc. lexp{1.2) | calc. lexp{1,2] 
Shower 2,9 |3,0+0,2] 4,1 |3,5+0,3 
Cascade 1,5 |1,44+0,1] 4,0 |4,1+0,5 
Evaporation 3,5 |3,3+0,1] 5--6 |6,1+0,6 
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TABLE 2. Mean Kinetic Energy of Particles (in 


Bev), Produced in One Act of Inelastic Collision 
of a Proton with the Mean Nucleus of a Photo- 


emulsion 
Particles Calcula- Experiment 
tions 1,2 
Shower protons 2,9 3,0+40,5 
Shower pions 0,52 0,63+0,1 [5] 
Cascade protons 0,15 0, 120-40 ,012 
Cascade pions 0,048 0,040-+0 ,003 
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Fig. 2. Angular distribution of cascade particles: a) group of light nuclei; b) group of 
heavy nuclei. 
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Fig. 3. Energy distribution of cascade particles (T is the kinetic energy): a) group of light 
nuclei; b) group of heavy nuclei. 


As can be seen, the experimental and calculated data for the shower and cascade particles agree well; the 
experimental data for the evaporation particles with black traces can also agree with the calculated data; however, 


40 
30 | 
20 
| 
= 10 q 
0,6 | 
a 
a 
151 


the results of the calculations depend to a large extent on actual assumptions as to the model of the nucleus. 


Together with arguments against the “mechanism of a tube” given in [1], the results of our calculations indicate 
that at an energy of 9 Bev and at lower energies a mechanism of intranuclear cascade is realized. 


LITERATURE CITED 


V. Barashenkov et al,Nucl. Phys. 14, 522 (1959/60). 

V. S. Barashenkov, et al., Atomnaya Energiya 1, 4, 376 (1959). 
N. P. Bogachev et al., Zhur, Eksp. i Teoret. Fiz. 38, 432 (1960). 
G. B. Zhdanov et al. , Zhur, Eksp. i Teoret. Fiz. 37, 7, 620 (1959). 

Yu. T. Lukin et al,, zhur, Eksp. i Teoret. Fiz. 38, 1074 (1960). 

N. Metropolis et al., Phys. Rev. 110, 185 204 (1958). 


i 
1. 
3. 
4, 
5. 
6. 
: 
: 
152 
iz 


MEASURING THE RADIATION CAPTURE CROSS SECTIONS 
OF FAST NEUTRONS OF 1??? 


Yu. Ya. Stavisskii, V. A. Tolstikov,and V. N. Kononov 


Translated from Atomnaya Energiya, Vol. 10, No. 2, pp. 158-160, February, 1961 
Original article submitted July 14, 1960 


The isotope of iodine 127 can be used as a standard [1,2] when carrying out activation measurements, It has a 
convenient half life (25 min), a fairly high fast neutron radiation capture cross section and a known thermal neutron 
capture cross section, 


However, until recently there were no data on fast neutron capture cross sections of ”' over a wide range of 
energies, Paper [3], the data of which were often used as a basis, is in sharp contradiction to the recently published 
data [4-6] and in all probability its results are erroneous, However, in [4-6] also,the range of energies of 0,01 -2.5 
Mev which is important in practical applications is not covered by measurements made with one method. The results 
of these papers also differ from one another, 


In the present work, carried out in 1958-1959, the dependence of the radiation capture cross section of neutrons 
with energies of 0.02-2.5 Mev of I!" were determined by the activation method [7]. Specimens of iodine and a fission 
chamber with a layer of U™® were irradiated together by a stream of fast neutrons. The applied 8 -activity was meas- 
ured by end-window beta-counters, The course of the cross section of radiation capture of 1’ neutrons was calculated 
from the dependence of the fission cross section of u™** on the neutron energy. Mean values of the fission cross sections 
of U5 were used (see the continuous curve in [8]). 


Beam of protons Beam of protons __ | 


7 2 3 3 2 1 
a b 


Fig. 1. Geometrical conditions under which the specimens and fission chamber were irradia - 
ted, 


The source of fast neutrons was the reaction T(p.n)He®, taking place in a Van de Graaf accelerator. A tritium 
target was used of 20 kev thickness. 


The irradiation of neutrons with energies > 300 kev was carried out at an angle of 0° with respect to the direc- 
tion of the beam of accelerated protons (Fig. 1a). The specimens were in the form of a disc of 20 mm diameter; they 
were placed at a distance of 45 mm from the center of the target close to the layer of u*®, trradiation by neutrons 
with energies of <300 kev were carried out at an angle of 100° with respect to the direction of the beam of accelera - 
ted protons (Fig. 1b). The specimens were in the form of rings of height 0.7 mm, external diameter 40 mm and internal 
diameter 28.5. The fission chamber was also in the form of a ring. 


The measurements at angles of 0 and 100° with respect to the direction of the proton beam were carried out 
with a certain “overlapping” with regard to the neutron energy. The coincidence of cross sections in the overlapping 
regions of energies indicated the reliability of the measurements, 
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The effect from neutrons scattered by the walls of the container did not exceed 0.3%, It was determined from 
the deviation of the law of neutron flux decay with distance from the law 1/R’, 


10 


0,01 2 3 5678901 2 8 45678940 
En, Mev . 

Fig. 2. The dependence of the neutron radiation capture cross section of 17 on the 

neutron energies, Data of particles: @—the present; ®-[1]; O—[4]; O-[15]; x —[6); 

A-[10}; O-[11]; +—-[12]; @-[13]; A-[14]. 


The contribution of neutrons scattered in the target backing was found by calculation. The correction to the ratio 
(Aji7/ Atjoas)” due to the scattered neutrons did not exceed several percent, which was confirmed by the coincidence 
of values of cross sections obtained in measurements at angles of 0-and 100°, Calibration measurements with thermal 
neutrons were carried out in the thermal neutrons was taken equal to 5.6 + 0.3b [8]. The fission cross section of a 
by thermal neutrons, equal to 582+ 4b, was taken from [9]. 


Errors in the obtained values of radiation capture cross sections of 1'7 neutrons were mainly determined by errors 
in the fission cross sections of U™*, which were ~25% for En < 60 kev, ~12% for E, = 60-150 kev. ~8% for neutrons of 
higher energies [5]. 


The results of the measurements together with the data of other authors are given in Fig. 2. 


In the investigated range of neutron energies the radiation capture cross section of “08d decays monotonously; 
within the limits of error in the measurement the experimental points lie on a curve of E~°*’, The smooth dependence 
of the radiation capture cross section of I' on the neutron energy is a very valuable factor when using this cross section 
as a standard in activation measurements. In the range of energies of 0.02-1 Mev our measurements agree well with 
the data of [5] within the limits. In the range of energies 0.15-1 Mev the values of the cross sections obtained in [4] 
are about 20% lower than our values, whereas in the region of energies 1-2.5 Mev they agree well. 


The data of [10, 11] (the data of the last paper were obtained by the transmission method in spherical geometry) 
are also in good agreement with our work, For neutrons with energies of 0.01-0,.2 Mev the cross sections given in [6, 
12,13] are about 15-20% lower than the values obtained in the present work, These difference may be due to errors in 
the fission cross sections of U%* which in this region of energies are known with a relatively poor accuracy. 
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A BETA-SOURCE BASED ON Au!*® FOR THE INVESTIGATION 
OF PHYSICAL PROPERTIES OF SUBSTANCES DURING IRRADIATION 


M. A. Mokul'skii and Yu. S. Lazurkin 


Translated from Atomnaya Energiya, Vol. 10, No. 2, pp. 160-162, February, 1961 
Original article submitted May 7, 1960 


For investigations into the action of radiation on a substance in many cases it is sufficient to have an apparatus 


giving in a small volume (about 100 mm’) an intensity of irradiation of the order of several hundred rads per sec- 
ond, 


A working instrument of this type is described below (Fig. 1). The source of radiation is a gold tip of a bronze 


needle, in the form of a tube of gold foil with walls of 0.2 mm thickness, external diameter 0.85 mm and 10 mm 
length, 


The needle with the gold tip 1 is fastened in an aluminum holder 2 and kept in a protective container of lead 
3 weighing 14 kg. 


The holder is fastened on the aluminum bar and is lowered into one of the channels of the reactor [1], where 
the thermal neutrons flux is 1.6-10'% cm~-sec-! (when operating at nominal power), After irradiation in this stream 
for 10° sec, 3-10" atoms of Au! appear in the gold tip, decaying according to the reaction Au! +g ~ + Hg'® with 
a half life of 64.6 hr and maximum energy of the 8 -spectrum of 0.97 Mev. The decay is accompanied by y -radia- 


tion with an energy of 0.411 Mev. The amount of Au’ obtained due to this irradiation in the reactor corresponds to 
an activity of ~2,5C, 


Irradiation of Specimens in a Working Container, The specimens are in the form of a tube 4 (see Fig. 1) of dia - 
meter 1-2/5-8 mm of the investigated material into which the gold needle is inserted. For low-temperature irradia - 
tion the specimens are placed at the axis of a Dewar vessel of foam plastic 5, placed in a steel working container 6 
with walls of thickness 8-100 mm, The protective container with the needle on a steel rod is placed on the working 
container and the needle is put into the specimen, Screw clamps 8 are used to hold the rod 7. 


With an activity of the source of 2.5 C, in the gold itself and in the specimen due to retardation of the elec- 
trons the energy liberated is 5-10* erg-sec™*. 


Electron Paramagnetic Resonance (EPR) in Irradiated Specimens, By means of a spectrometer similar to that 
described in [2] with high frequency modulation of the magnetic field measurements were made of the EPR spectra 
in specimens irradiated by the method described at room temperature and at 77°K [3]. The EPR spectra of polymethyl 
methacrylate, teflon and polyethylene, irradiated at a temperature of 20°C with different doses have the usual form. 


Using the EPR, measurements were made of the dose distribution along the radii of the specimens, A strip of 
teflon of thickness 0,1 mm was coiled in a tube similar in dimensions to the usual specimens, After irradiation the 
strip was cut up into pieces corresponding to different layers and using EPR the number of peroxide radicals in each 
piece was measured, The number of radicals is proportional to the dose obtained by the material. This holds for doses 
(of the order of several megarad), which were used in our experiment, The dose distribution along the radius of the 
specimen is shown in Fig. 2. The path of the fastest 8 -electrons of Au! in teflon does not exceed 2 mm [2]. Conse- 
quently, the layers of specimen at distances greater than 2 mm from the surface obtain a dose only from the y -rays. 
The absolute value of the dose from y -rays was measured by dosimetric instruments at distances (from the radiator) 
of 250 mm and more. Knowing the geometry of the source it is possible to calculate from these data the dose obtained 
from y -rays by any layers of the specimen (see Fig. 2, broken curve) and also to obtain the absolute value of the dose, 


corresponding to the points of Fig. 2 lying between 2 and 3 mm, Consequently, it is possible to calculate the absolute 
value of the total dose in each layer (the continuous curve of Fig. 2). 


It has been found that outside the source ~10% of the energy of 6 -electrons is liberated, i.e. a layer of speci- 
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Fig. 1. of "gold needle" ap- 


paratus, 
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Fig. 2. Dose distribution along the radius of the teflon 
specimen irradiated by the gold needle, 1) distance from 
the inside surface of specimen (tube); —total dose; 
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Fig. 3. EPR spectrum of polymethyl methacrylate 
during the irradiation process; A) spectrum; B) 
dependence of the value of the signal on time, 
The arrow indicates the instant of removal of the 
needle from the specimen, 


men of thickness 1 mm obtains a mean dose of ~200-250 
tad/sec. The contribution of the y -rays to the total dose 
is ~16%, for the whole teflon specimen of the given 
thickness. 


* Similar experiments were first conducted in the Institute of Chemical Physics, Academy of Sciences, USSR [4] using 


a linear electron accelerator at 1.6 Mev. 
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ferent materials* . 


The specimen in the form of a tube was placed in a Dewar vessel in the resonator of the EPR-spectrometer. The 
needle on a long bar was placed inside the specimen, This led to a slight detuning of the resonator (no greater than on 
the introduction of ordinary specimens), Figure 3 shows the EPR spectrum in irradiated polymethyl methacrylate and 
the change in value of the signal during irradiation and after its completion. 


Creep in Teflon Under Irradiation, In [3] there was a description of the reversible acceleration of relaxation 


Observing EPR in the Process of Irradiation, The small dimensions of the source (gold needle) and the compara- 
tively high intensity of irradiation make it possible to use this source to obtain EPR spectra in the irradiation of dif- 


processes in polymers under the action of irradiation in a reactor, It was of interest to find out whether this effect 
would be observed under the action of irradiation by light particles when the number of produced nuclear recoils is 
negligibly small compared with the case of irradiation in a reactor with high fast-neutron fluxes. 


During irradiation by electrons the needle was placed parallel to the plane of the specimen at a distance of 1 
mm from its center, The dose obtained by the specimen was distributed unevenly throughout its thickness and length. 
The mean dose rate throughout the volume of the specimen was 10-20 rad/sec, The needle was brought to the speci- 
men periodically (six times in all), The tensile deformation of the specimen was recorded continuously* by the épp- 
09 recording apparatus, The obtained curve is given in Fig. 4. Despite the certain number of random inflections on the 
curve (connected with friction in the instrument) the effect of reversible acceleration of creep can be clearly seen. 


The described experiment was not intended to give quantitative results, However, it is fundamentally important 


and indicates that the reversible acceleration of relaxation processes in polymers is not connected with the specific 
action of heavy particles, 
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Fig. 4. Creep in teflon. Stress 0.35 kg/mm?, tempera - 
ture 18°C, The instants of switching on and switching 
off the irradiation are shown by arrows (top-switching 


on). 
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* The measurement was made by M. E. Fiveiskii on an instrument which he designed for the remote recording of creep 


3. A. G. Kiselov, M. A. Mokul'skii,and Yu. S. Lazurkin, 
4, Yu. N. Molin et al., Doklady Akad. Nauk SSSR 123, 


5. M. A. Mokul’skii, Yu. S. Lazurkin and M. B. Fiveiskii, 


The small dimensions of the sources based on mdio- 
active gold, the ease of protection from soft y -radiation, 
the comparatively long half life and the large activation 
cross section make these sources very convenient for much 
work on the irradiation of materials. Using these sources it 
would presumably be possible to irradiate materials at the 
temperature of liquid helium simultaneously with the re- 
cording of the EPR spectra or with other measurements, The 
configuration of the radiators for many experiments should 

be different (gold disc, two parallel gold discs, thin-walled 
tube, etc.), The choice of the shape of the source is especial- 
ly important in those cases where it is essential to have ap- 


proximately the same dose rate throughout the whole volume 
of the specimen, 
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A GENERATOR PRODUCING A HIGH FLUX OF 14 
OR 2.5 Mev NEUTRONS 


V. 


I. 


Petrov 


Translated from Atomnaya Energiya, Vol. 10, No. 2, pp. 163-166, February, 1961 
Original article submitted September 1, 1960 


High-voltage accelerators of 150-200 kev are generally used for the production of 14 and 2.5 Mev neutrons, In 
the majority of cases, the maximum D-T neutron flux is slightly more than 10” n/sec. An increase in flux to 10" 
n/sec and more (in the case of 2.5 Mev neutrons, to more than 10° n/sec) introduces certain technical difficulties 
which are associated both with the production of a sufficiently intense beam of deuterons in the accelerator tube and 
with the limiting thermal conditions for tritiated and deuterated targets. 


In this letter a neutron generator is described which produces a D-T neutron flux greater than 10"! n/sec and a 
D-D neutron flux greater than 10° n/sec. 


The construction of the generator is shown in the figure. A high-frequency-discharge ion source is used in the 
generator, The basic characteristics of the source are as follows; deuterium ion current, 9-11 ma; deuteron content, 
~89o; deuterium consumption, ~53 cm¥hr; length of ion extraction system channel, 10 mm; channel diameter, 3.5 
mm; extraction voltage, 7—8 kv; high-frequency-generator power consumption, 500-600 w; generator frequency, 25 
Mc. The power supply circuit for the source is assembled in a separate unit, Transmission line voltage to the unit is 
fed through four transformers which maintain the voltage at 250 kv (transformer ratio is 1:1). Operational control of 
the source is carried out from a control panel by adjustment of the voltages which are fed into the primary windings 
of the transformers. During operation, the ion source is continuously cooled by a stream of air from a fan. 


The construction of the neutron generator acceleration chamber provides a minimum distance from ion source 
to target, That type of construction simplifies the problems of focussing and accelerating an intense ion beam, The 
chamber is made up of two sections in which the electrodes for the focussing and accelerating gaps are located. A 
porcelain cylinder serves as insulator for the focussing gap. The accelerating gap electrodes are located within an 
insulator made of a plastic similar to lucite, The shape and location of the electrodes keep scattered charged particles 
from striking the internal insulator surfaces. 


High vacuum is produced in the accelerating chamber by a diffusion pump (pumping speed 1000 liters/sec) and 
a forepump. The diffusion pump is connected to the chamber through an angle valve on which an oil vapor trap cooled 
by liquid nitrogen is built. The 300-liter forevacuum tank is also used as a base for the acceleration chamber. The 
working vacuum in the chamber with an ion beam current of several milliamperes amounts to about 1 x 10-5 mm Hg. 


An electromagnet which produces a transverse magnetic field across the path of the ion beam is placed at the 
exit of the acceleration chamber, The electromagnet winding is made of one layer of water-cooled copper tubing 
with a 5x3 mm bore, The need for the introduction of a magnet became apparent during the first experiments on 
ion beam acceleration. With an ion beam current of several milliamperes, the absence of a magnetic field led to the 
formation of an opposing electron current in the acceleration chamber comparable in magnitude to the ion current, 
The electron current upset the stable operation of the ion source; in a number of instances, it melted source walls, 
caused electrode heating, loss of vacuum, and high-voltage breakdown. The introduction of a transverse magnetic 
field very definitely removed the opposing electron current and associated phenomena, In addition, target operating 
conditions have been improved by the introduction of an electromagnet since, in the presence of a transverse field (of 
a particular strength), the target is freed from thermal loading by molecular ions in the beam, and is significantly less 
contaminated by deposits of organic vapors. 


The accelerator high-voltage supply is obtained from two high-voltage devices: for the accelerating gap, from 
a 200 kv constant voltage generator at 20 ma load, and for the focussing gap, from a 30 kv, 20 ma rectifier. The high 
voltage generator is based on a one-stage voltage doubling circuit. The basic elements of the generator are a type 
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Construction of Neutron Generator (dimensions in milli- 
meters): 1) Ion source; 2) focussing gap; 3) insulator; 4) 
accelerating gap; 5) electromagnet; 6) target. 


IM-100-0,1 condenser, a KR-220 kenotron, and a type 
IOM-100/25 high voltage transformer. For the high-volt- 
age supply of the focussing gap, a high-voltage rectifier 
based on a full-wave bridge circuit employing a KRMM- 
110 kenotron is used, Voltage to the rectifier is fed through 
a special 200 kv transformer. Under operating conditions, 
the voltage to the focussing gap electrodes is 23-25 kv, 
that to the accelerating gap electrodes is 170-180 kv, The 
ion source extraction voltage of 7-8 kv must be added to 
that, 


Zirconium targets, saturated with tritium or deuter- 
ium on silver or copper backings 25 mm in diameter and 
1,5-2 mm thick are used in the neutron generator, During 
operation, the targets are cooled by a stream of running 
water, Despite a satisfactorily high rate of cooling, the 
maximum permissiable ion current at the target at full 
accelerating voltage does not exceed 2,5-3 ma. An 
increase in the permissible current has been achieved by 
a circular displacement of the beam on the target through 
the use of a 4 kv, 5 kc, a.c. voltage, The voltage is fed, 
with phase shift, to two deflection plates located ina 
special exit structure, In this way, the permissible current 
at the target has been increased to 5 ma, 


The neutron generator was placed in operation in 
1959. Assembly and adjustment of the apparatus were 
carried out with the assistance of V. P. Zyuzin, R. N. 
Morozov, and V. P. Suslov. The following data, character- 
istic of machine operation, were obtained: maximum flux 
of 14 Mev neutrons, 5.3x 10" n/sec; maximum flux of 
2.5 Mev neutrons, 5.0 x 10° n/sec; operating deuteron cur- 
rent at the target, up to 5 ma; full accelerating voltage, 
200-210 kv; maximum current of focussed deuterium ions 
at full accelerating voltage, 10 ma. 
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THE EFFECT OF RADIATION ON THE ELECTROCHEMICAL 
BEHAVIOR OF 1Kh18N9T STEEL 


V. V. Gerasimov and V. N. Aleksandrova 


Translated from Atomnaya Energiya, Vol. 10, No. 2, pp. 164-166, February, 1961 
Original article submitted September 1, 1960 


The aim of this work was to study the electrochemical behavior of IKh18N9T steel during irradiation by thermal 
neutrons, 


In the tests a special electrolytic glass cell was used (Fig. 1) of diameter 20 mm and height 400 mm. The thick- 
ness of the cell wall was ~1 mm, In the upper part of the cell there was a calomel electrode and an electrolytic switch, 
Tests were carried out on specimens of 1Kh18N9T steel. The specimen dimensions were 1 x 10 x 235 mm. The chem- 
ical composition of 1Kh18N9T steel is as follows; 0.07 weight% C; 1.23 weight.% Mn; 19.1 weight.% Cr; 10.5 weight. 

%o Ni; 0.53 weight.% Ti. 
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Fig. 1. Electrical cell; 1) Specimen; 2) auxilliary Fig. 2, Anodic and cathodic polarization of 
electrode; 3) calomel electrode; 4) electrolytic 1Kh18N9T steel in a 0.01 N solution of Na,SO,: 
switch; 5) tube; 6) ampoule; 7) case; 8) to circuit; ; —@-—under conditions of irradiation at 80-90°C 
9) to potentiometer. (—-@-—— curves taken under the same conditions 
in the reverse order; —x——without irradiation at 
85°C (—-x———curves taken under the same condi- 
tions in the reverse order). 
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The auxiliary specimen was a 1-mm diameter wire of 1Kh18N9T steel insulated from the investigated specimen 
by a 5-mm diameter glass tube, The investigations were carried out in a 0,01 N solution of sodium sulfate and a 0.01 
N solution of sodium chloride. The volume of liquid poured into the cell was ~20 ml. 


The current intensity was measured by the M-82 and M-91 milliammeters, the potential was measured by the 
LP-5 potentiometer. The current was provided by dry BAS-80 anode batteries. The electrolytic cell was protected from 
mechanical damage by a cover made of foil (1Kh18N9T steel) with wall thickness 0.1 mm. Between the cover and 
the wall of the cell near the middle part of the specimen there were the junctions of a chromel-copel thermocouple, 
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Fig. 3. Anodic and cathodic polarization of 1Kh18N9T 
steel in a 0.01 N solution of NaCl; -®——under condi- 
tions of irradiation at 80-90°C; ——x——without irradia - 


tion at 85°C, (For the explanation of the broken curves 
see the caption to Fig. 2). 


potential and temperature were measured every 10 min, The current was assumed to be steady if, for a given value of 
the potential, it did not change in 1 hr. At first the cathode polarization curve was taken and then the anode polariza- 
tion curve, Three to four hours after the end of polarization the stationary potential was again measured. The obtained 


polarization curves are shown in Figs. 2 and 3. 


The irradiation did not change the kinetics of the anode process of IKh18N9T steel in a 0.01 N solution of sod- 
ium sulfate. The rate of the anode process in the passive region and in the region of superpassivation, the values of 
the breakdown potential during irradiation are the same as without irradiation. 


Before the test the surface of the specimen was carefully cleaned with No, 180 polishing paper. The potential 
of the specimen was compared with that of a saturated calomel electrode. To avoid contamination of the investigated 
solution with potassium chloride the electrolytic switch (after drawing in potassium chloride) was washed with the 
working solution on the side of the specimen, The distance from the end of the electrolytic switch to the specimen in 
a given experiment does not affect the measured potential; with maximum removal of the end of the switch from the 
specimen (by ~5 mm) the potential decreases only by 0,015 v compared with the potential measured with minimum 
distance, The potential is practically unchanged along the length of the specimen; with a current of up to 100 wa the 
potential at the ends of the specimen is 0,02 v greater than the potential of the middle part. 


The Stationary Potential of 1Kh18N9T Steel at 80-90°C 
in Irradiated and Nonirradiated Solutions (Volts with 
Respect to a Standard Hydrogen Electrode (SHE)) 


Without irradiation} During irradiation 
before j3 hr before 3 hr 
Medium polariza-| after lariza -| after 
tion larization lariza- 

ion Fion 


0,01 NV 
solution 0,133 0,318 0,403 0,673 
NagSU, 


0,01 NV 
solution 0,083 0,243 0,503 0,583 
NaCl 


In order to study the electrochemical behavior of 
1Kh18N9T steel during irradiation the cell was placed in 
the active zone of the reactor. The experiments were car- 
ried out with a flux of 10” neutrons/cm?-sec and a cell 
temperature of 80-90°C, Each experiment was performed 
twice. For comparison the experiments were repeated with- 
out irradiation in a thermostat at 85°C, The results of the 
measurements are given in the table. 


The kinetics of the electrode processes were studied 
by the potentiostatic method. One hour after immersing the 
cell in the channel of the reactor the stationary potential 
was measured and the external current applied. The current 


: 


In a 0,01 N solution of sodium chloride during irradiation the character of the anodic polarization curve of 
1Kh18N9T steel is preserved. The rate of the anode process in the passive region is practically unchanged. However, 
the breakdown potential and the section of the anodic polarization curve corresponding to the superpassivation region 
are displaced toward the positive side by 100-150 mv. 


In both tested solutions during irradiation the sections of the cathodic polarization curves corresponding to the 
ionization of oxygen were displaced toward the positive side. The latter fact points to an increase in the rate of the 
cathode process, which was caused by the presence in the solution of radiolysis products—short-lived radicals and hydro- 
gen peroxide. 


An increase in the rate of the cathode process leads to displacement of the stationary potential of steel toward 
the positive side (see table) due to which the region of the passive state decreases. In solutions where the region of 
the passive state is small, for example in chlorides and especially in the presence of tensile stresses in the metal, the 
displacement of the stationary potential toward the positive side under the action of irradiation can bring steel into 
the superpassivation region, The latter fact can lead to an increase in the effectiveness of the corrosion process. During 
irradiation in solutions of sodium sulfate a stationary potential is set up corresponding to the region of the passive state, 
in a solution of chlorides a stationary potential is set up close to the breakdown potential or corresponding to the region 
of superpassivation (see table and Figs. 2 and 3). 
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A METHOD OF INVESTIGATING PROCESSES OF RETARDATION 
OF FISSION FRAGMENTS IN METALS AND ALLOYS 


N. A. Protopov, Yu. B. Shishkin, V. M. Kul*gavchuk, 
I. 


and V. Sobolev 


Translated from Atomnaya fnergiya, Vol. 10, No. 2, pp. 166-168, February, 1961 
Original article submitted November 21, 1959 


The growing interest in the problem of making fuel elements, especially for the operation of nuclear reactors 
in energy-stressed systems and with deep combustion of the fissile isotope, requires a more detailed study of the ele- 
mentary processes of retardation of fission fragments. Up to the present these processes have been insufficiently studied 
and it has only been established that there are regions of a thermal peak and their volume [1]. Some work was carried 
out in 1946 [2] with a noncollimated beam of fragments; in 1957 work was published [3] in which, using the mass 
spectrograph method for the separation of fragments, paths and retardation curves were determined for a number of 
metals such as aluminum, nickel and gold. To measure the energy of the fragments a CsI(T1) scintillation counter 

was used, the sensitivity of which was determined by means of a-particles of polonium, However, the complexity of 
this method unfortunately limits its wide application when studying the retardation of fragments in a number of metals 
and alloys which are of the most interest for the technology of fuel element preparation. 


In application to the metallurgy of nuclear fuel the problem of developing a similar method is considerably 
simplified if it is borne in mind that in this case considerable importance attaches to the duration of the regions of 
the thermal peak and the density of energy liberation, depending on specific losses in energy of the fragment, which 
in the final analysis characterize the scale of local damage in a solid and the rate of change in its physicomechan- 
ical properties, The difficulties belonging to problems of nuclear physics and connected with the identification of 
fragments by the masses and the nuclear charge are therefore eliminated since the specific losses in energy and paths 
of the fragments can be determined directly experimentally. Therefore, for the problems under consideration of con- 
siderable importance is the curve of specific losses in energy of the fragment as a function of its path. 


AAAAAAA 


Fig. 1b. Arrangement of experiment for study - 
F 7 Mt ing specific energy losses of fission fragments 
in metals by the method of substitution. 


Figure la shows a scheme of experiments which demonstrate the 
Fig. 1a. The main arrangement of the essence of the method developed by the authors, The main part of the 
apparatus for studying the passage of scheme is the ionization fission chamber. The impulses formed in both 
U fission fragments in metals and channels of the chamber are fed through amplifiers I and III to a circuit 
alloys. of double coincidences II with a resolving time of 0.5 usec, The impulses 

are counted by the circuit IV with a mechanical counter. The fission chamber 


is made of brass and has a leakproof valve (not shown on the diagram) to 
evacuate the system and fill it with argon. 


The chamber consists of the following elements (see Figs. la): 1—-UgOg preparation (~ 0.1 mg/cm?, enriched 
up to 90% with U™*), applied to aluminum foil of thickness 0.4 mg/cm? and placed in the middle of the chamber; 
2~—collimators; 3—foil stack with variable thickness of the investigated metal; 4—collecting electrodes. 
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Fig. 2. Dependence of the specific losses of 
a light group of fission fragments on the path 


in gold;——curve taken from [3]; data 
of authors, 
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Layer of material, mg/cm* 
Fig. 3. Dependence of yield of fission 
fragments on the path in gold: a) curve 
taken from [2]; b) data of the authors. 
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Layer of material, mg/cm? 
Fig. 4. Dependence of yield of fission fra- 
gments on the path in aluminum; a) curve 
taken from [2], plotted with an:ionization 
fission chamber; b) curve taken from [2], 
plotted by radiochemical measurements; 
c) data of authors, 


The thickness of the collimators and the pressure of the argon 
(~ 100 mm Hg) in the chamber are chosen so that the energy losses 
of the fragment on passing through the collimator can be neglected 
(~3-5 Mev). The use of a fission chamber and coincidence circuit 
eliminates errors connected with counting false impulses arising 
from the y -quanta and the recoil nuclei. This means that all the 
measurements can be made in the nuclear reactor with a thermal 
neutron flux, changing according to the conditions of the experiments 
in the limits 10°- 10° neutrons/cm?.sec, 


We will illustrate the operation of the circuit with the example 
of liberation of fragments with the same path, we will determine 
their path and the specific energy losses. We will suppose that in 
each channel of the chamber behind the collimator there are foil 
stacks of the investigated material of the same thickness, We will 
designate the initial number of coinciding impulses by No for zero 
thickness of the material. If d is increased the same in each channel 
then the number of coincidences will decrease. For a certain value 
d=d,, the number of coincidences will be equal to zero or a very 
small part a of No, depending on the sensitivity of the circuit, In 
this case d,, will presumably characterize the total path of the con- 
sidered fragments. 


To measure the specific energy losses in this case we proceed 
in the following way. In both channels the same foil stacks are placed 
(with thickness d=d,) of a material for which the curve of the depen- 
dence of specific energy losses on the path of the particles is known, 
Then in one of the channels all foils of the material with known 
losses, starting with the first, are successively replaced by foils of 
the investigated material (Fig. 1b); by selecting the thickness, con- 
stancy of the initial number of coincidences is obtained. Therefore, 
by the substitution method we determine the specific losses in the 
investigated material for the whole path of fragments, Throughout 
all measurements the second channel is kept unchanged. 


To determine the path and specific losses of energy of the 
heavy and light groups of fragments separately, as in the previous 
case, we first determine Ny. Then in one of the channels the foil 
stack of the investigated material is dismantled until the number 
of coincidences becomes small compared with Np. In this case the 
fragments with the maximum path will presumably be liverated. 
Depending on the selected value of the number of coincidences a 
the investigations are carried out for the same spectrum width of 
the fragment paths, Then, in the second channel of the chamber 
the foil stack of the investigated material is dismantled until the 
number of coincidences is reduced to zero. In this case the number 
of foils in the second channel characterized the total path of the 
corresponding heavy fragments. A similar substitution method is used 
to determine the specific energy losses of heavy and light groups of 
fragments separately, using known data on losses in any metal, for 
example,in aluminum, 


Figure 2 shows the curve for specific energy losses during the 
retardation of fragments of a light group in gold, which we plotted 
with the described method, In these measurements a= 0.05 and fra- 
gments were liberated with a maximum path, lying near the value 
16.5 mg/cm*, For comparison we give a curve of losses plotted by 
means of the graphic calculation of dE/dx from the energy-path 
curve taken from [3]. As can be seen, the results agree well. 


_ 


Figure 3 shows a curve for the passage of fragments b, which we obtained for retardation in gold. For comparison 
a similar curve a is given, taken from [2]. As can be seen, , the curves differ considerably. In the first place the maxi- 
mum thickness of gold needed to absorb 100% of the fission fragments of U5 is not 9 mg/cm? (as in [2]) but 16.5 
mg/cm?®, We notice that our measurements compare best with data of a later work [3], according to which the path 
for gold was equal to 13 mg/cm?. The low value of the path in gold obtained in [2] is presumably due to the use of 
an insufficiently developed apparatus and method, On the other hand, our results agree with those of recently published 
work [4], from which it follows that the number of very light fragments with maximum paths is much greater than 
would be expected on the basis of earlier radiochemical measurements. In the second place, curve a is displaced some 
distance to the left. This is presumably due to the fact that in [2] an uncollimated beam of fragments was used. 


Figure 4 gives a similar curve c for passage of fragments in aluminum and for comparison curves a and b are 
given, taken from [2]. The above remarks refer mainly to this case, 


We notice that the curves for passage of a collimated beam of fragments take on a new importance, In fact, if 
the curve for passage in the material is determined by a certain function N = f(x), then the expression 
'max 
+ \ xf' (x) dx 


min 


(where lyin and 1,4, are the minimum and maximum paths of fragments in the substance respectively; No is the 
total yield of fragments) is the total path taken by all fragments formed in unit time in unit volume of the solid. For 
a known value of the effective cross section S of the physical process accompanying the passage of a fragment, for 


example the process of homogenization of an alloy, as in [5], the product LS =dv/dt is the rate of homogenization of 
the alloy, 
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THE MELTING POINT AND OTHER PROPERTIES 
OF THE LOWER OXIDES OF NIOBIUM 


O. P. Kolchin and N. V. Sumarokova 


Translated from Atomnaya Energiya, Vol. 10, No. 2, pp. 168-170, February, 1961 
Original article submitted October 4, 1960 


Several lower oxides of niobium have been reported in the literature (Nb, NbO, Nb,O,, Nb,O3, NbO», Nb,O, 
and NbgO;), although as a result of the investigations [1,2] the existence of only two of them has been established; the 
oxide NbO and the dioxide NbO,. A knowledge of the properties of these oxides is essential both to determine their 
own possible applications and also to obtain an accurate conception of the processes taking place during the reduction 
of Nb,Os, in which the lower oxides occur as intermediate compounds, and of the processes for removing oxygen from 
niobium by vaporation in the form of lower oxides by heating the metal in vacuum, Yet contradictory and even in- 
accurate information about the properties of the lower oxides is frequently met with in the literature. For example, a 


“trioxide” of niobium, NbOg, with a melting point of 1772°C has been described [3], although in actuality this is not 
a separate oxide but a two-phase product. 


Published data on the true lower oxides of niobium are very meager [2,4-8]. 


We have attempted to determine a number of properties of these lower oxides, 
Initial Materials 


Niobium dioxide was prepared by reducing the pentoxide with lamp black in vaccum at 1300 and 1700°C, and 
the oxide was prepared by the same method from the dioxide. The initial pentoxide contained (wt.-%); <0.5 TagOs, 
0.06 SiOg, 0.1 TiO,, and 0.04 FeOg; during reduction in vacuum the greater part of the iron, silicon, and titanium 
was removed, The three specimens of dioxide that were investigated contained (wt.-%); 0.02-0.04 C and 0.006 N; they 
had the compositions NbO, 942, NbO; 9g, and NbO, 9g, which all lie within the homogeneous range of the dioxide [5]. 
X-ray analysis showed that the product was single-phase and had the parameters a_= 4,824 0,02 kX and c= 2,994 0.02 
kX, which agree with the values published in [2]. 


The specimens of niobium oxide contained (wt.-%): 0.04-0.06 C and 0.03 N and had the compositions NbOp,95, 
NbO, 1, and NbO, 99, which all lie within the homogeneous range of the oxide [5]. 


Melting Point 


The dioxide and oxide were ground to powder of particle size <0.15 mm and compressed under a pressure of 
1.5-2,.0 tons/cm? into compacts 10 mm in diam and 10 mm high. Two compacts of the oxide or dioxide were placed, 
one on top of the other, in a container, 17 mm in diam and 100 mm high, made of niobium sheet. Then the container 
was placed inside the graphite heater of a vacuum (0.5-1.0" Hg) electric resistance furnace in such a way that the 
molten top portion of the upper compact was most strongly heated. By arranging the change in this way, it was possible 
to avoid interaction between the molten oxide and the metal and carbon heater, It also ensured suitable conditions for 
measuring the temperature of the molten oxide by means of a disappearing-filament optical pyrometer without the 
error [3] due to the unknown emissivity of the oxide, since the distance from the surface of the oxide to the top of 
the container was almost five times the diameter. To confirm that the true temperature was being measured, compacts 
of the oxide, lamp black, and graphite were placed in the same container and in separate experiments theirtemperature 
was measured in the region close to the melting point of the oxide. In spite of the considerable difference in the 
emissivities of carbon and niobium oxide, the difference between the measured temperatures did not exceed 10°C, 
i.e. it lay within the limits of experimental error of measurements made with the pyrometer. A correction was made 
in the pyrometer readings only for the light absorbed by the observation window in the vacuum furnace. 


Under the conditions described, the melting point of the specimens of niobium oxide investigated was 1935°C, 
and that of the dioxide 2080°C. The approximate error in the temperature determination was 15°C. 
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Figure 1 shows a photomicrograph of the melted portion of the oxide compact, and from this it can be seen that 
there was no second phase present in the oxide whose melting point had been determined. The single-phase character 
was also confirmed by x-ray analysis. This showed that the product has a parameter a = 4,201 kX. It is to be noted that 
even with a small excess of oxygen in the product it was possible to observe the presence of a second phase, niobium 
dioxide, which forms a eutectic with the oxide (Fig. 2); in the x-ray phase analysis, the sensitivity of which is ~1%, 
this phase was not observed. 


Fig. 1, Structure of pure melted niobium oxide. 
Etched in HF + HNOg (3:1) for 3 sec, x200, Fig. 2, Structure of melted niobium oxide containing 

niobium dioxide. Etched in HF + HNOg (3:1) for 3 sec. 
x 200, 


The discrepancy between the values found for the melting point of the dioxide and those given in the literature 
[9] is no doubt due to the fact that in [9] the specimens of dioxide used contained niobium oxide as impurity, and 
hence the data obtained from them do not relate to the true melting point of the dioxide. 


It is of interest to note that the melting point we have determined for niobium dioxide is close to the erroneous 
value determined in 1907 for metallic niobium (1950°C) [7] and for long regarded as reliable (see, for example, [10]). 
It is evident that the product of the reaction of niobium with oxygen was used to determine the melting point of the 
metal; this is quite possible since the melted oxide is very difficult to distinguish externally from the metal; it has a 
metallic luster and in external appearance is similar to an ingot of niobium, though very brittle. 


Volatility. Composition of the Gaseous Phase 
Weighing the upper compact after heating in vacuum (0.5-1 ¢ Hg) at various temperatures (accuracy of weighing 
0.0002 g) showed that the oxide and dioxide of niobium began to evaporate at a significant rate at 1700°C, At 1850°C 


the whole of the oxide and 45 wt.-% of the dioxide had evaporated after 4 hr, and the whole of the dioxide compact 
had evaporated after 8 hr. 


A mass-spectroscopic investigation of the dioxide specimens we obtained showed [11] that the gaseous phase 
over the dioxide contained only NbO, molecules, i.e. there was neither dissociation nor association of the molecules. 
Microhardness 

The microhardness of various specimens of melted niobium oxide, the compositions of which lay within the 
limits of homogeneity of the oxide, was 1930 kg/mm? and the mean value of the microhardness of the eutectic (see 
Fig. 2) was 794 kg/mm? (load = 50 g), The microhardness of the dioxide was 1720 kg/mm’, although this value may 
be inaccurate since small cracks appeared in the dioxide crystals during the microhardness determinations. 


Character of the Electrical Conductivity 
Niobium oxide is characterized by the electrical conductivity of metals, while the dioxide is a semiconductor, 
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Oxidation in Air 

Studies were made of two specimens of dioxide, one of which was obtained at 1300°C and had a composition 
corresponding to the formula NbO2,99, while the other was obtained at 1700°C and had the composition NbOj,99. After 
grinding, both specimens formed black powders with a particle size 0.15 mm. A specimen of niobium oxide with the 
composition NbOpo 9:4 was obtained at 1700°C and ground to a gray powder haveing a metallic luster and a particle 
size of 0.15 mm. Samples of the oxide were heated in air at a temperature of 100-300°C, After this oxygen content 
was determined from the increase in weight on heating in air at 800°C; the accuracy of this method of determination 
is ~4 0.1%, The oxygen content of niobium oxide did not change on heating for up to 6 hr at 100, 150, 200, 225, 250, 
and 275°C; the oxygen content of the dioxide increased slightly at 275°C. However, at a temperature as low as 150°C 
the surface of the particles of oxide and dioxide took on a yellowish tinge and at 200°C a bronze color, Both the oxide 
and the dioxide were completely oxidized to the pentoxide after 6 hr at 300°C, 


The author wishes to thank L. V. Mel'nikova for carrying out the metallographic analysis of niobium oxide. 
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THE HARDNESS OF SOME NIOBIUM-BASE ALLOYS 
AT HIGH TEMPERATURES 


I. I. Kornilov and R. S. Polyakova 


Translated from Atomnaya Energiya, Vol, 10, No. 2, pp. 170-172, February, 1961 
Original article submitted June 18, 1960 


One of the principal factors in raising the high-temperature strength of alloys is the formation of multicom- 
ponent solid solutions. On account of the complexity of the chemical composition of the solid solution, and increase 
takes place in the strength of the chemical bond between the atoms of the various elements, 


The object of the present work was to investigate the change in hardness of niobium and its alloys with tempera- 
ture. The "hot" hardness method described in [1] was used in the experiments. 


As a result of these comparatively short-time hardness tests at high temperatures, estimates of high-tempera- 
ture strength of the alloys were obtained which were in the same order as those obtained by long-time creep tests 
(2). 


The capacity of niobium to form solid solutions is used in practice to strengthen it at high temperatures and 
to produce commercial alloys. For the investigation 2-to-6- component alloys were synthesized, the compositions 
being calculated to fall within the single-phase and supersaturated -solid-solution regions. The composition of these 
alloys (in wt.-%) are given in the table. 

The alloys were made from the following materials: 99.2% 
metallic niobium, 99.9% molybdenum, 99.8% iodide zirconium, 
commercial silicon, 99.9% aluminum, and carbon. The alloys 
were melted in an arc furnace with a nonconsumable tungsten 
electrode and an argon atmosphere. To avoid metal losses the 
alloys were pressed into 20-g charges before melting. To obtain 
uniform ingots the specimens were turned over and remelted 
several times, 


Mo 
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Cylinders 14 mm in diam and 5 mm high were machined 
from the specimens obtained, All the specimens were given a homogenizing anneal in a TVV-2 vacuum furnace at 
1600°C for 10 hr. These conditions were sufficient to bring the specimens into an equilibrium state. 


Microstructures of annealed alloys are given in Fig. 1. 


Hardness measurements were carried out on a VIM-1 machine in vacuum at temperatures of 20, 100, 200, 300, 
400, 500, 600, 700, 800, 900, and 1000°C. The diamond pyramid was under load for 1 min, Three impressions were 


made on each specimen at each temperature, The diagonal of the impression was measured on a PMT-3 microhard- 
ness machine, 


Figure 2 gives curves es the temperature dependence of the hardness of pure niobium, with a room-temp- 
erature hardness of 13.8 kg/mm?, The hardness curve for niobium decreases linearly with rise of temperature to 400- 
600°C and then falls more steeply at 600-1000°C, At the latter temperature the hardness of niobium was 10 kg/mm’, 
With greater complexity of chemical composition the hardness of the alloys increases, as can be clearly seen from 
curves 2-6, The course of curve 2shows that the introduction of a second element raises the hardness of the alloy, 
particularly at high temperatures. A similar change in hardness was observed in the three-component alloy. As the 
number of components in the alloys increased, the absolute value of the hardness rose markedly, However, the soften- 
ing process proceeded somewhat more rapidly than in the case of binary and ternary alloys. 


In order to demonstrate more clearly the part played by alloying in strengthening the alloys, a characteristic 
was introduced which indicates the degree of strengthening in relation to the composition of the alloy [3]. As a 
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Fig. 1. Microstructures of annealed niobium-base 
alloys: a—90% niobium, 5% molybdenum, 5% 
zirconium (X 315);b—89%% niobium, 5% molyb- 
denum, 5% zirconium, 1% silicon (X 315). 
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Fig. 2. Temperature dependence of the hard- 
ness of niobium alloys: 1) niobium; 2) two- 
component alloy; 3) three-component alloy; 
4) four-component alloy; 5) five-component 
alloy; 6) six-component alloy. 
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Fig. 3. Values of the coefficient of strengthen- 
ing at 20, 800, and 1000°C, Alloys: 1) two- 
components; 2) three-component; 3) four- 
component; 4) five-component; 5) six-com- 
ponent, 


measure of the degree of strengthening a coefficient was adopted, 
in the calculation of which the hardness of niobium was taken 
as unit at any temperature in the range from 20 to 1000°C, 


Figure 3 shows the values of the coefficient of strengthen- 
ing at 20, 800, and 1000°C. The two-component alloy had a 
coefficient 1.5 at room temperature (20°C), 5.0 at 800°C, and 
10.8 at 1000°C, The coefficient increased with rise in the number 
of components in the alloy. Thus, the coefficient of a six-com- 
ponent alloy was 2.96 at 20°C and 19,3 at 1000°C, i.e, at 
1000°C the six-component alloy was 19.3 times stronger than 
niobium, 


From the foregoing the following conclusions may be 
drawn: 


1) Niobium, like nickel, iron, and cobalt, can be strength- 
ened by alloying with elements which enter into solid solution 
with it or form supersaturated solutions with separation of excess 
phase. 


2. The “hot” hardness methodenables the high-tempera - 
ture strength of the alloys to be determined to a first approxi- 
mation. 


3. Multicomponent alloys having great hardness at 1000°C 
require further investigation and testing by the standard method. 
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THE CHARACTERISTICS OF IRRADIATED GLASSES* 


Zdenek Spurny 


Institute of Nuclear Studies, Czechoslovakia Academy of Sciences, Prague 
Translated from Atomnaya Energiya, Vol. 10, No. 2, pp. 172-173, February, 1961 
Original article submitted September 22, 1960 


In a number of papers [1-5] resuits have been given of investigations into the dependence of darkening of glass 
after irradiation with a certain dose on the composition of the glass, the effect of impurities in the glass on the forma- 
tion of color centers,and descriptions have been given of the dosimetric properties of glasses of different compositions. 
To compare the action of ionizing radiation on glass of different types it is essential to determine certain indices. The 
use of Czech optical glasses for dosimetric purposes [6] has shown that the number of color centers forming in the glass 
during irradiation depends on the following factors: the dose of irradiation, the energy of the radiation, the chemical 
composition of the glass. If these factors are constant then the number of color centers forming in the glass also depends 
on the temperature, the dimensions of the glass and the time after irradiation, 


In [7] four indices are suggested for the quantitative determination of the radiation effect in glasses; the index 
of stability of glass to the action of ionizing radiation K, (the logarithm of the dose at which the light transparency 
is reduced by 5%), the index of saturation K,,, (the logarithm of the dose at which the glass is saturated with color 
centers), the minimum value of light transmission T,,j, (transmission in the region of saturation by color centers) and 


the coefficient of intensity of darkening Q (the ratio of the coefficient of light transmission before irradiation to the 
coefficient of transmission after irradiation at a dose of 10° eB, 


To measure these indices the following conditions must be satisfied: 


1, The source of irradiation should be Co® of about 100 r/min intensity in order to always observe the same 
conditions in measuring the energy and intensity of the source. 


2. The temperature of the glass should be kept constant (20°C) in order to eliminate its effect. 


3. The transmission of light should be measured not later than 10 min after irradiation, in order to eliminate 
the effect of fading (the dependence of “breakdown” of the color centers on the time after irradiation). 


4. Glass specimens of the same thickness should always be used (preferably more than 1 mm) in order to elimin- 
ate the effect of scatter. 


5. The light transmission should always be measured with a source of monochromatic light of a certain wave- 
length (for example A = 3500 A), 


The index K, should be called the "threshold of noticeable changes” since it refers to the minimum measured 
number of color centers, This term is used in the literature and it should not be confused with the concept of stability 


of color centers, The index K,,, should be determined as the point of intersection of two tangents to the curve for 
dependence of light transmission on the dose of irradiation. 


For the total quantitative characteristics of the radiation effect in glasses it is essential to introduce another 
very important constant, characterizing the stability of the color centers which form. It is well known that after forma- 
tion color centers “breakdown,” Curves for the breakdown of color centers have been obtained experimentally for a 


number of optical glasses of different compositions. It has been found that the fading of glasses is best described by 
the empirical formula 


1 
log 


* This article was received from Czechoslovakia. 
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Indices of Czech Glasses After Irradiation 


Indices 


T of glass 
80 ath Ky | Q | K 


nu |Tmin| Ty 


SF-10 3,3 4.143 
F-8 3,1/90 14,3) 4 |2,4 6 
BK-12 3.4119 14,4] 4 | 3,4] 25 

Ordinary 
40 sodium 3,6141 14,6| 4 14,4] 35 
ZK-5 3,719 4,7 1 4,8 42 
FK-5 4,4| 3,615.44 4 |5,5 120 
0 20 40 60 80 


Curve for breakdown of color centers. 


where yp is the coefficient of the transmission of light, measured 1 hr after irradiation ur is the same coefficient 
measured 7 hours after irradiation (T<T; Ke is a constant). 


The figure shows the curve for breakdown of color centers for ordinary sodium glass, obtained from formula (1), 
and experimental points. The constant K, in formula (1) can sufficiently well characterize the stability of color 
centers (darkening of glass) which are more stable the greater the values of K,; this stability can also be characterized 
by the time (in hours) needed to break down 25% of the formed centers (the constant T,). Then, 100% of the number 
of color centers should be taken as the light transmission of glass after irradiation with a dose of 10°r, 


Observing the mentioned conditions of measurement, using a Co™ source we irradiated several types of Czech 
optical glasses of the standard Jena composition, we measured the indices of the author of [7] K,, Ksaz, Q, Tryin and 
the author of the present work (K,, T,). The results of these measurements are given in the table. The measurements 
show that the indices suggested in [7] for glasses sensitive to radiation are apparently difficult to measure. A more ac- 
curate determination of these indices should therefore be given. 
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THE BUILD-UP FACTORS FOR HETEROGENEOUS SHIELDING 


L. R. Kimel’ 


Translated from Atomnaya Energiya, Vol. 10, No. 2, pp. 173-175, February, 1961 
Original article submitted June 15, 1960 


The energy build-up factors of y -irradiation of Co® in heterogeneous (double layer) shielding were determined 
for a plane-parallel beam normally incident on the shielding, the beam being provided by the method described in 
[1]. The detector was a halogen counter with filters [2]. The intensity was determined graphically by integrating its 
distribution before and after passing through the shielding barrier. The shielding materials were plates of lead, iron 
and aluminum measuring 75x 75 cm. The build-up factors were determined for the combinations Pb + Al; Al+ Pb; 

Pb + Fe; Fe + Pb; Fe+ Al; Al+ Fe. The material occurring first in the combination was nearer to the radiation source. 


Tables 1, 2 and 3 give experimental values of the build-up factors (energy) for the given combination of mater- 
ials. 


TABLE 1, Build-up Factors for the TABLE 2. Build-up Factors for the TABLE 3, Build-up Factors for the 


Combinations Pb + Al (top values) Combinations Pb + Fe (top values) Combinations Fe + Al (top values) 
and Al+ Pb (bottom values) for and Fe+ Pb (bottom values) for and Al+Fe (bottom values) for 
hv)= 1.25 Mev hv 9= 1.25 Mev hv 9= 1.25 Mev 
Thick- | Thickness of material Thick- Thickness of material Thick- ickness of material 
oa : placed nearer to the ness of : placed nearer to the ness of : Rist nearer to the 
teri materia material! detector, Ux 
detector, ux placed detector, plac 
nearer nearer nearer 
to the 4 5 to the 5 to the 415 
source, source, source, 
x x 
|6,6 2,3] 2,8]3,5] 4,4]5,1 2, 


3,4 


3,8 


4,1 


From an analysis of the tables it can be seen that the build-up factor for Pb + Al combinations can be determined 
as the product of the build-up factors for lead and aluminum of the appropriate thicknesses (in ux). This determination 
of build-up factors for the cases where the light material follows the heavy material was recommended in [3]. For the 
combinations of Al+ Pb the build-up factor is somewhat greater than for lead of a thickness (in x) equal to the thick- 
ness of Al+ Pb, For the combination Fe + Al and Al+ Fe at small thickness the build-up factors almost coincide, 


It can readily be shown that the build-up factor for double-layer shielding can be determined from the express- 


hv 
Jr (ey) e [hy, (hv) dhv 
By, Byy + 


: 
| 
| 5,4) 6,9 18,4 9 5,2 |6,3 9 | 3,3/4,1]5,4| 6,6 17,0 
2,7|3,113,3] 3,6 |3,8 2,0] 2,3] 2,6] 3,0] 3,3 3,3/4,115,2| 6,016,9 
3 |3,7/5,0]6,7| 8,5 |10,3 3 3.3] 4,2|5,2| 6,3]7,4 3 7,7 19,3 
3,0/3,4|3,7| 3,8] 4,0 2,5] 2,8]3,2| 3,5]3,8 4,215,3|6,3| 7,2|8,0 
4,4|6,0| 8,0] 10,2/12,5 ‘ 3,815,016,1| 7,5|8,8 5,316,5|7,81 9,4 |10,7 
3,2|3,613,9] 4,0) 4,3 3,113,4]3,7] 4,114,4 5,316,4|7,4] 8,2] 9,0 
5,3] 7,0|9,4 | 12,0]14,6 4,3|5,8|7,1] 8,7 |410,2 5 | 8.7]8,0] 9,3] 10,6/12,2 
° 4,3) 4,6 13,6] 4,014,3| 4,6 | 5,0 6.417.518.5| 9.4110.3 
| 
ion 
(1) 
114 


TABLE 4, A Comparison of the Experimental (Upper) The expression (1) makes it possible (for cases where 
and Calculated (Lower) Values of the Energy Build-Up the spectrum of the scattered radiation is known after pass- 
Factor for the Combination Pb+ Al ing through the first barrier) to calculate the build-up 
factors. This calculation was performed for the combina - 
tion Pb+ Al, the spectra of radiation scattered after the lead 
Pb, ux { | 2 | 4 being taken for an infinite medium [4]. With this selection 


Thickness of Thickness of Al, ux. 


the error is small, since the build-up factor for a barrier 
of thickness 1 ux is only 5% less than the build-up factor 
: nae ot for an infinite medium at a depth 1 ux (for hv = 1 Mev) 
— [5]. In the calculations no allowance was made for the 
3,0 4,4 6,9 angular distribution of radiation leaving the first layer, 
- 2,74 3,92 6,46 which is presumably not very large. In Table 4 the results 
_" of the calculation are compared with experimental data. 


The differences between the calculated and experimental 
data, mainly caused by the “rigid operation” of the detector 

in the region of soft energy, do net exceed +15%, The good coincidence of experimental and calculated data for the 
combination 5y xFe +1.85 j:xPb is shownin[6]. From this work it can also be seen that the ratio of the energy build-up 
factors of radiation from a point isotropic source for the combination 1.85, xPb +5y xFe and 5uxFe+ 1.85" xPb is 1.5, 
which agrees well with our results (6,3 +4.0 = 1.58). The coincidence of the calculated and experimental data means 
that this method can be used to calculate build-up factors for arbitrary energies and material if the spectrum of scat- 
tered radiation is known, 


The author would like to thank O, I. Leipunskii for his constant interest in the work. 
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SOLUTION OF THE KINETIC EQUATION FOR A MEDIUM 
WITH A POINT MONODIRECTIONAL SOURCE 


E. B. Breshenkova and V. V. Orlov 


Translated from Atomnaya Energiya, Vol. 10, No, 2, pp. 175-177, February, 1961 
Original article submitted August, 1960 


The problem of the transmission of radiation in an infinite medium with various sources has been studied in 
many articles [1-3], Point isotropic sources and infinite plane sources with various angular distribution of y -radiation 
and neutron emission were studied. In such cases, due to the symmetry of the problem ,the radiation flux depends 
only on three variables, which fact significantly simplifies the solution of the kinetic equation. The solution of this 
equation for a point monodirectional source is associated with five variables. The steady problem of the transmission 
of radiation through a medium with a point monodirectional source was examined in general form in [4], but the 
direct solution of this problem entails considerable mathematical difficulty due to the large number of variables. 
Satisfactory quantitative results have not been obtained. 


In solving the single-velocity problem for a medium with point monodirectional source, I. I. Bondarenko used 


the reversibility of single-velocity diffusion and the problem was reduced to the solution of the kinetic equation for 
a medium with isotropic source. 


The diffusion of radiation in a slowing medium (due to elastic and nonelastic collisions of neutrons with nuclei 
and Compton scattering of y -rays by electrons) is an irreversible process, However in this case the problem may be 
solved by using a more general property of the kinetic equation, namely the reciprocity theorem [5,6]. 


In the present paper the solution of the kinetic equation for a medium with point monodirectional source is 
reduced with the aid of the reciprocity theorem to the solution of the conjugate kinetic equation for a medium with 
an isotropic source. This equation contains three independent variables and may be solved by the method of moments. 


Suppose G(r, £, Q; ro, M,Qq) is Green's function for the kinetic equation describing the flow of radiation with 
energy E and direction® at the points (with the condition that at the point fo there is a unit point source with energy 


Ep and direction® ). If the function G‘(r,//,L; to,9,Qo) is Green"s function for the conjugate kinetic equation, then 
by the reciprocity theorem 


G(r, E, To, Ep, Q,)=G* (To, Eo, «o, I, E, Q). 


Integrating this equation with respect to dQ, we obtain 


Go (r, E; To, Ey, (To, Eo, Qo; E), (1) 


where Go—\ GaQ. 


Thus if we must know the distribution of total flux Gp from a monodirectional point source, then it is sufficient 
to find the solution of the conjugate equation (including the angular distribution) for a medium with isotropic source. 


In the case of y -radiation the kinetic equation of a medium with a point monoenergetic monodirectional source 
has the following form: 


Qgrad F(r, Q, A)= —p(A) F(r, Q, A) + \ A) an’ 8(1—Q-Q' 


a—2 
F (r, A’)4-8 9) 8 (@— Qo) 8 (A— Ay), 


where F (r, Q, 2.) is the fediation flux at the pointr with wavelength A and direction of velocity 9; K(A', A)= 
(VX)? E (A’—A) +-(A’—A)? ] is the probability density of the Compton scattering with change 


‘ 
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in wavelength from \' to A; u (A) is the coefficient of absorption, Together with this equation we may, corresponding 
to formula (1), solve the conjugate kinetic equation for a medium with isotropic source: 
@, 2) = —w (A) 9,9) + (a, O°, x 
1 


(2) 


In this the distribution of total y -radiation flux in a medium with a point monodirectional source may be found with 
the formula 


Fo(r’, ’)=F* (to, Qo, Ao): 


Performing the substitution r -'=R and considering the symmetry of the source, from which it follows that F+(R, Q, A4)= 
=F*(R, &,A), where € = QR)/R, we obtain together with formula (2), the following equation: 


a 
3 
xF*(R,§", 2" 4-2) 4 2 (3) 


We expand the conjugate flux in a series of Legendre polynomials 


substitute formula (4) into equation (3), multiply the relation obtained by #(€) and integrate through by dQ. Asa 
result we obtain an equation for Fj* (R, 4) 


at2 (5) 


Multiplying equation (5) by R" and integrating it over all of space, we obtain an equation for the moments bj,, ,(A) 


a+2 


\ R°F}(R, 4) 4nR@R, 
0 


It is easy to see that system (6) determines only the moments by (A), for which l<n and the indices? and n have 
the same parity. We seek Ftz(R, ) and F$7,4,(R, A) in the following form: 


N 
uR 


n=-0 


(6) 
where 
177 


N 
—wR 
Dy Cater (A) HR), 
n=0 


where the coefficients ag; , , (A) and coi 41, » (A) are solutions of the systems of N+ 1 equations; 


N 


n=0 


N 
1)! 


n=0 


form=f 41,2 


The solution of system (7) allows the expression of the coefficients ay .,,(A) and Cg + 4,9(A) in terms of the 
even-even and odd-odd moments found from equation (6), It is easy to see that the number N is limited by the num- 
ber of known moments, The more moments that are known, the more accurate may the quantity FT (R, A) be calcula - 
ted. 

To solve the problem considered, it is useful to exclude the unscattered flux of radiation from the equation for 
moments (6), which may be done in the usual manner. 

In conclusion the author expresses thanks to V. F. Turchin for many useful remarks in the preparation of this 
article, and to G. I. Marchuk and Sh. S. Nikolaishvili for fruitful discussions of the results. 
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THE EFFECT OF INELASTIC SCATTER OF NEUTRONS 
IN URANIUM ON THE MODERATION LENGTH IN WATER 


B. A. Levin, E. V. Marchenko, and D. V. Timoshuk 


Translated from Atomnaya Energiya, Vol. 10, No. 2, pp. 177-179, February, 1961 
Original article submitted July 7, 1960 


In the calculation of uranium-water reactors an important part is played by the moderation length of the neutrons 
in the lattice. Measurements and calculations [1-3] of the moderation length, performed for different concentrations 

of uranium and different diameters of blocks in the lattices, show that the moderation length is considerably affected 
by the inelastic scatter of neutrons in uranium. The values of the neutron ages [T = 1/6)#”] determined in assemblies 
using a thin source (convector) and in exponential experiments in subcritical systems, differ greatly. Besides certain 
reasons,this difference can be explained by the fact that in experiments of the first type there is practically no in- 
elastic scatter in the source,whereas in experiments of the second type the inelastic scatter in the “parent” block can 
noticeably change the primary energy of the neutrons. 


In the present work direct measurements were made of the effect of inelastic neutron scatter in uranium on the 
moderation length in pure water, For the most complete exclusion of the effect of elastic scatter the measurements 
were made in spherical-symmetrical geometry, The source of diameter 16 mm, imitating the spectrum of fission 


neutrons, was surrounded concentrically by spherical layers (of thickness 2 cm) of metallic uranium depleted with 
isotope. 


The mean square of the distance traveled by the neutrons from the point isotropic source on moderation to an 
energy of, for example, 1.46 ev is determined by the expression 


A(x) rtar 
(1) 
A(r) r°dr 


where A(r) is the activity of the indicator (indium foil in a cadmium can) as a function of the distance r from the 
source, The application of expression (1) to the distribution of moderation density in this particular case gives a value 
of r? which is determined not only by the energy spectrum of the fast neutrons, but also by the dimensions of the cavity 
formed by the sphere in the moderator. In order to exclude the effect of dimensions of the cavity the measurements 

of r? were performed with uranium layers of different radius a at a constant thickness of the layer. By extrapolating 
the relationship r’(a) to a= Othe value of r? was determined corresponding to a point surce of neutrons. To study the 
effect of dimensions of the cavity on the value of r? with change in a over wider limits and to obtain r § for neutrons 
ofthe initial spectrum,measurements were made with thin-walled hollow spheres. 


The measurements were made in a water-filled tank of diameter 100 cm and height 110 cm, The spheres were 
suspended on two steel wires of diameter 0.5 mm along which was lowered a light plexiglas holder until it was arres- 
ted by the sphere; the holder had grooves for locating the casettes with indium foils (at a given distance from the 
surface of the sphere), The casettes were sealed with wax. The foils of thickness 70 mg/cm? had a working surface 
diameter of 1.7 and 2.5 cm. Suspending the foils in the neutron field showed that the activities of indicators of the 
same diameter differed from the mean value by not more than 0.5%, As a rule, the targets were exposed in two pieces 
with a distance between them of not less than 8 cm. The activity of the foils was measured simultaneously from both 
sides by two scintillation counters (FEU -29 with a crystal of stilbene of diameter 32 mm and thickness 1 mm) in lead 
shielding. The sensitivity of the apparatus was checked by a preparation of Sr” and was kept at one level with an ac- 
curacy of +0.5%, The measurements at all distances were repeated from two to six times. At distances up to 15 cm 
from the surface of the uranium sphere the measurements were made by small targets, at distances of 10-21 cm by 
large targets in cadmium (thickness 0.5 mm) casettes and at distances of 15-31 cm by large targets in aluminum 
cassettes, In order to obtain the density distribution of moderation of the neutrons to an energy of 1.46 ev over the 
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whole range of distances, separate sections of the curve were placed together at the points where the measurements 
overlapped. In measurements with thin-walled spheres small targets in cadmium casettes were used. In addition, for 


a cavity with a= 3,3 cm, measurements were made with targets of different diameters in cadmium and aluminum 
cassettes. In both cases within the limits of error in the measurements (1.2%) the same values of T? were obtained. 


The values of r? obtained from relationship (1) from the space distributions A(r) are shown in the figure and in 
the table. Integration of the curves of A(r) r” and A(s)r‘ at distances greater than r-a=15 cm was performed analytical- 
ly. In accordance with the results of the measurements it was assumed that the decrease in Apr’, starting with r-a = 
15 cm, becomes exponential with a relaxation length of 7.114 0.08 cm for hollow spheres and 6,52 4 0.06 cm for 
uranium, The relaxation lengths were determined by the method of least squares, 


The Values of r? Measured at Different Radii a of the Cavity in Water 


Cavity With uranium With lead 


3,3 5,0 7,0 3,59 5,80 6,80 3,56 6,81 


244,8+ | 272,64 | 311,24 206, 3+ | 234,14 | 255,14 | 242,64 | 305,24 
+1,9 | +5,3° | 7 +1,8 | | 42,1] | 43,8 


0,588-+-0 ,006 0,480+0,027 0,582-0,039 


r2, 187,1+1,1 161,544,0 182,8+5,7 


*k— coefficient for a in equation (2). 


The neutron density was assumed to be constant within the 
spheres. As can be seen from the figure, {¥? is well approxi- 
mated by a linear dependence on a 


Vr +ha, (2) 


which permits reliable extrapolation of the value 1 # and 
a=0, Fromthe curve 1 for neutrons of the imitating source 
which we used, r? is equal to 187,14 1.1 cm’, 


The results of measurements with layers of uranium are 
shown by curve 2. The difference in the slopes of curves 1 and 
2 is presumably connected with the softening of the spectrum 
Lo with inelastic reaction of neutrons with uranium nuclei. The 
correction for capture of the resonance neutrons in uranium 
was determined by comparing the course of density of neutrons 
with an energy of 1.46 ev near the uranium and hollow spheres 
01 2 82 4 § &§ 7 8 Y 1 4a,cm on the assumption that the distributions of densities in the 
absence of resonance absorption are similar. The correction 
is 0.5% of the value of r*, The linear extrapolation of curve 
2 to a=0 gives 161.54 4.0 


Dependence of JF? on the radius a of a spherical 
cavity in water; 1) (O)-without scatterer; 2) (@)— 
with uranium; 3) (Li)—with lead, 


Measurements with lead (see curve 3 of the figure) were 
made to find the effect of elastic scatter in a spherical layer on the moderation length, As shown by the measurements 
made by another method, inelastic scatter of fast neutrons is very small for the source which we used in lead. A com- 
parison of curves 1 and 3 shows that the difference in the values of 7? obtained with spherical layers of lead and hollow 
spheres is very small and is within the limits of error of the experiment. Consequently, the observed reduction inf 
for uranium is not connected with elastic scatter of neutrons in the spherical layer. 
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In measurements with uranium there was a 9% increase in the number of neutrons connected with the fission 
of uranium mainly by fast neutrons (the uranium and lead spheres were in cadmium cans). Since the mean paths of 
fission neutrons and neutrons of the surce in uranium are close, in view of the closeness of the mean energies the 
effect of inelastic scatter will be approximately the same and the addition of fission neutrons is equivalent to an 
increase in source strength. A decrease in f” due to fission neutrons which have not undergone inelastic scatter in 
uranium due to the difference in the ages of the fission neutrons (28 cm?, see below) and neutrons of the source which 
we used (31 cm?) is 1%, Therefore, allowing for this correction the inelastic scatter of neutrons in a layer of metallic 
uranium of thickness 2.0 cm leads to an increase in F? of (12.7 4 2.2)%. 


On the basis of known cross sections, estimating the fraction “of inelastically scattered neutrons in a 2.0 cm 
thick layer of uranium as being equal to 0.21 and using the relationship f?= ar*jnej, +(1—a)rg, we find that for in- 
elastically scattered neutrons the second spatial moment inel is equal to 64 cm?, From the relationship T (E) [4] the 


mean energy of inelastically scattered neutrons in uranium is ~0.5 Mev. 


Curve 1 (see figure) gives a quantitative allowance for the effect of dimensions of the source on T , obtained in 
[5], for fission neutrons. 


If we assume the end of the cylindrical tube with converter (U™) to be a sphere of the same radius, then for T 
we obtain the value 28+ 1.5 cm?, which is closer to the theoretical value T die” 26 + 0.5 cm? [4] for neutrons of the 
fission spectrum in water. : 


The results of this work show that in uranium-water systems with blocks of sufficiently large dimensions a de- 
crease in the moderation length of fission neutrons due to inelastic scatter in uranium is considerable and should be 
allowed for in reactor calculations, 


We would like to thank G. A. Bat" for discussing certain theoretical problems, L. E. Morozova, G. S. Stolyarova, 
and L. A. Serdyukova for helping with the measurements and calculations. 
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An international conference on radioisotope applications in the physical sciences and industry, organized by the 
International Atomic Energy Agency (IAEA) in collaboration with UNESCO, met in session from September 6 to Sept- 
ember 16, 1960, in Copenhagen, Participants at the conference included over 550 persons from 40 nations and 11 inter- 


national bodies; the number of reports heard was 150 (including two evening lectures), 20 of which were delivered by 
Soviet scientists, 


The reports submitted to the conference can be grouped according to subject matter under 12 different headings. 
Below , we endeavor to make a brief summary of the most interesting papers, from our viewpoint, and the scientific 
and engineering data contained therein, grouping them under several broad headings. 


Radioisotopes in Geophysics and Meteorology 


Fifteen reports presented in these fields took up problems related to meteorite studies, propagation and migration 


of radioisotopes in the atmosphere and in the sea, carbon-14 dating, and applications for nuclear-geophysical techniques 
in mineral exploration and research on underground waters, 


At the first session of the conference, papers presented by Soviet scientists were "Carbon-14 Dating” by A. P. 
Vinogradov, and “Liquid Scintillators for Radiocarbon Dating in Archeology” by I. E. Starik and colleagues, These papers 
made available data on a carbon-14 measurement procedure worked out in the Soviet Union for ascertaining the a- 


mount of C™ present in archeological and geological specimens, and information on age determination of numerous 
geological and archeological objects in the USSR. 


A paper presented by O. A. Schaeffer and associates (USA), “Meteorites as interplanetary space samples, for 
studying cosmic rays," proved that the solar component of the cosmic radiation exerts no discernible effect on the 
composition of meteoritic matter, on the basis of a study of radioactive isotopes of different half-lives found in metero- 
ites. K. Govel and P. Schmidlin of CERN presented a wealth of empirical data on determinations of the “cosmic® age 
of meteorites, i.e. the duration of their irradiation in cosmic space, The values derived ranged from 2 to 200 million 
years, W. Gehr et al. (West Germany) presented data on an application of the Re/Os-technique to age determinations 
of iron meteorites, Although the data reported by the authors on the age of iron meteorites, 5.6 to 8.7 10° years, is 
acknowledged to be of preliminary and tentative character, they are of great interest, since determination of the age 
of iron meteorites by other methods has been hampered by great difficulties. G. Wilson and A. Mac Near (Britain) 
presented findings on a determination of the half-life of Rb®™, reported as (5.254 0.1)-10® years, is the first value de- 
rived by the direct method; it is of great importance for calculations of the absolute age by the Rb/Sr method. 


It is fitting to mention (bearing the foregoing discussion in mind) that papers on studies of the content of radio- 
isotopes in meteorites are viewed as highly significant from the standpoint of age studies and the elaboration of theories 
bearing on the formation of the matter comprising the planetary system. This question, aside from its value as added 
knowledge, is of direct practical significance in that the conquest of outer space necessitates a knowledge of the 


changes produced by cosmic radiation in matter found in outer space; and some idea of the intensity of the cosmic 
radiation proper is also highly desirable. 


Interesting and highly useful data on the irregular latitude distribution of C™ over the earth's surface obtained 
on the basis of four years of measurements (1956-1959) by G, Tauber (Denmark) and generalizations of the findings 
reported by other workers were made available in Tauber's paper entitled "The effect of geographic latitude on the 
migration of radioactive carbon from the stratosphere to the troposphere.” The author in particular demonstrated that 
the pathways of migration and the mechanism allowing access from the stratosphere to the troposphere for fallout iso- 
topes resulting from nuclear bomb explosions may be traced by determining the c™ content in contemporary flora 
specimens collected in various sites, This is of great interest both for prognoses of radiation hazards applicable to dif- 
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ferent regions and in relation to utilization of the radiocarbon method for determining absolute age. The author's 
recommendations on the feasibility of measurements of c™ content in atmospheric carbon dioxide at a large num- 
ber of sampling stations in order to determine c" transport from stratosphere to troposphere also merits attention. 
This report and two review papers presented by the World Meterological Organization showed the feasibility of utiliz- 
ing radioisotopes to estimate displacement and the character of the transport of air and water masses over the earth's 
surface. This avenue of research unquestionably holds great promise not only for the solution of purely meterological 


aerological, and hydrological problems, but, as has already been stated, is also useful for prognoses regarding radio- 
active fallout levels. 


Problems related to the use of tritium as a research tool in underground water studies were discussed in papers by 
Soviet (F. A. Alekseev et al.) and Netherlands (E. Kupper and H. Hout) authors. 


A report by Yu. P. Bulashevich et al. (USSR) “Nuclear geophysics in exploration of ore and coal deposits," in 
which practical techniques for radioactive logging which have found many applications in ore and coal fields through- 
out the Soviet Union, was met with deep interest. 


Radioisotopes in metallurgy and solid state physics 

These areas were the subject of 15 papers presented. A paper by L. Lacombe et al. (France) discussed the ques- 
tion, associated with severe technical difficulties, of the separation and separate investigation of self-diffusion with- 
in grain boundaries and between grain boundaries, The paper demonstrated the fact that the formula derived by the 
Soviet scientist P. L. Gruzin can be extended and that, by measuring the total activity in iron specimens, the frac- 
tions of activity associated with self-diffusion within and between grain boundaries can be successfully distinguished. 


A paper by C. Besnar and J. Talbot (France) discussed the use of isotopes in studying sulfur penetration in to iron 
when the iron has absorbed hydrogen from a cathode. The authors found an intensive penetration of the sulfur into 


the bulk phase of hydrogen-saturated iron, with the sulfur following a path along rain boundaries. This paper was of 
corsiderable practical interest, 


A report by P. Albert and J. Betti (France) was devoted to the elaboration of techniques for determining trace 
impurities in high-purity materials ( metallic and nonmetallic alike) for those instances where determination of 
such impurities by chemical and spectral analyses is not feasible. The possibility of determining traces of 45 elements 


(the number may be extended to 60 as more experience is acquired) was demonstrated by the example of neutron- 
irradiated aluminum and iron, 


J. Crawford and J. Cleland (USA) explained the results of detailed investigations of thermal neutron effects on 
germanium and its properties (magnitude and mode of germanium conductivity), It was shown that of five german- 
ium isotopes, three isotopes are transformed to chemical impurities (Ga"!, As”, se™) which, upon forming in the 
germanium lattice, effect a change in the electrical properties of germanium, Data on the number of atoms of the 
isotopes formed and their half-lives are presented. A similar investigation was also carried out with respect to silicon 
and indium, The use of this method holds great promise for the study of numerous complex semiconductors (inter - 


metallic and chemical compounds) which present great hindrances to the introduction of certain interesting im- 
purities by the chemical pathway. 


A paper by E. G. Miselyuk et al. (USSR) communicated findings of a study using radioactive isotopes Agi? 
Fe®, ca!®, Te!” and Ta'® of diffusion and solubility of the corresponding isotope species in germanium. The interest 
placed in diffusion and solubility studies of these elements is particularly due to the fact that their impurity centers 
in germanium are electrically active and capable of affecting both conductivity and charge carrier recombination, 
by greatly increasing the rate of recombination. Soviet scientists also submitted papers dealing with diffusion, re- 
crystallization, and phase transformations (S. Z. Bokshtein et al.), measurement of self-diffusion constants of cations 
insilicate melts with Na”4, K®, and Ca® employed as tracers (V. I. Malkin and B. M. Morgunov), determination of 
the thermodynamic activity of antimony in a-iron (I. A. Tomilin), Tomilin advanced a new method based on a 
study of the distribution of antimony between two phases, liquid lead and solid iron. The gist of this method is that 
a specimen of thin iron foil 10-30 microns thick and weighing about 0.1 g is contacted with molten lead containing 
a predetermined concentration of Sb! tagged antimony. On being left to stand at a preset temperature, the anti- 
mony diffuses into the iron, ending up evenly distributed throughout the entire thickness of the foil. Use of the radio- 
active antimony isotope has made it possible to accurately measure the concentration of that element in specimens 


weighing less than 0.1 g. Quantitative determination of the composition of such tiny weights by straightforward 
chemical analysis is extremely difficult. 


The study of lattice imperfections is one of the most burning problems in solid state physics. Methods widely 
used in such studies are the method of chemical additives and radiation effects on the lattice structure. A paper by 

G. Lambe (USA) discussed several possible uses of radioactive disintegration for bringing about imperfections in solids. 
The possibility of using this method for the above purpose was demonstrated by the author with tritium and krypton as 
examples, The method advanced may find a wide range of applications in obtaining and investigating free radicals in 
organic substances which exhibit reasonably high stability to radiation effects, 


Radiation Detectors, Monitoring and Measuring Instruments, Activation Analysis Tech- 
niques 

Under this heading, one paper of interest with respect to technique used is a communication by P. Keukelsbergs 
and associates (Belgium) in which data are given on specially selected counter gas mixtures allowing the possibility 
of measuring tritium in methyl chloride in a Geiger counter with an excellent plateau (length 1300 v, slope 2%). In 
a paper presented by D. Sprecklen and G, Stevens (Britain) on the use of autoradiography for quantitative analysis of 
radioactive specimens, a novel technique is advanced for preparing sample preparations with the use of post-scoring 
techniques, as well as a technique for recording many specimens on a single film. L. Wiesner and associates (West 
Germany) described an ingenious modification of the procedure of vapor chromatographic analysis using pure argon 
flow counters operating at reduced pressure in the region between ionization chamber and counter modes of operation. 


Papers on activation analysis, especially several contributions devoted to the use of a new technique for carry- 
ing on activation analysis, are of some interest, 


A paper presented by K. Lungren (Sweden) described a method of activation analysis which allows identifica - 
tion of isotopes in a complex mixture of radioactive substances without prior chemical separation, an accomplishment 
of enormous value for rapid execution of a large number of analyses. The paper presents a theory to explain the method 
and cites examples of the determination of several elements in mixtures. An arrangement consisting of a dual-channel 
facility with Nal-crystal scintillation counters. The channels are included in a triple-coincidence unit, at the output 
of which is an electronic circuit which determines the amplitude of the incoming signal. 


A report by B, Carr (USA) is devoted to the use of pulsed neutron sources (based on the T (d, n) He reaction) in 
activation analysis. Miniature accelerators developed by the authors produce neutrons at 14 Mev energy and at an 
intensity adequate for the formation of the required numbers of short-lived radioisotopes in the test samples. A report 
by W. Meinke (USA) wad devoted to the use of a thermal neutron generator, producing a flux of 10° to 10° neutrons/ 
/cm*-sec, for the same purpose, and to the use of pneumatic devices on reactors with a flux of 10” neutrons/cm?*-sec. 


It is worth noting that sources of that type may be found applicable in stationary facilities for neutron produc- 
tion and in portable instruments, particularly when the latter are used for well logging. The use of generators such as 
those proposed by Meinke broadens the field of applications for activation analysis, bringing it to the level of a routine 
industrial method, and the use of a pneumatic rabbit makes it possible to carry out an analysis for short-lived isotopes 


and in combination with high-speed spectrometers, i.e. makes it possible to shorten the time required and to enhance 
the sensitivity of the method. 


The report submitted by W. Kuykendall and R. Wainerdi (USA) is highly interesting in this light. They developed 
system of automated activation analysis incorporating an anticoincidence gamma-ray spectrometer with a 256- 
channel analyzer and an IBM-1704 electronic computer, This automatic counting arrangement takes into account the 
character of the gamma-ray spectrum and the decay time of the spectral components, and allows high-speed process- 
ing of the data for specimens containing a considerable number of different isotopes. At the present time, it takes 15 
min to analyze a specimen containing 15 components. Kuykendall and Wainerdi worked out routines for qualitative 
and quantitative automatic analysis of specimens containing as many as 25 unknown constituents. 


Radioisotopes andChemistry 


Almost one third of the reports read at the conference (49 reports) were devoted to the use of radioisotopes in 


chemistry. Solutions of analysis problems facilitated by the use of various radioisotope techniques furnished the sub- 
ject matter for a large grouping of reports. 


K. Evans and J. Harrington (Britain) described an interesting procedure followed in determining hydrogen content 


in metals by isotopic dilution with tritium. Experimental data on hydrogen content in beryllium, aluminum, and 
uranium specimens were reported. 


A paper submitted by D. Ulen and S. Ziegler (USA) dealt with a determination of the composition of matter in 
the form of a gas or a solution containing compounds with reducing and oxidative properties, by binary chemical ex- 
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change. The method consists in evolving the radioactive gas Kr® from Kr® quinolate, and the content of the unknown 
compound is calculated from the amount of gas evolved. The method is interesting, novel, and may be found useful 
in industry. The authors employed this method for a meterological analysis of ozone in the upper layers of the atmo- 
sphere and for an analysis of air pollution by sulfur dioxide gas. 


A. Schrodt (USA) proposed a variant of the radiometric method of chemical analysis based on measurements of 
the radioactivity of an excess of titrating reagent. This method results in fairly high accuracy, is simple in its experi- 


mental procedure, lends itself readily to automating adaptations, and is expected to meet with a wide variety of practi- 
cal applications, 


Modifications of the isotope dilution method devised to facilitate identification and quantitative estimating of 
small amounts of chemical substances were also interesting as described in several reports. In particular, a paper presen- 
ted by R. Bailey (Britain), devoted to the use of isotope dilution in protein analysis, offers some new approaches in 
methods for determining racemization of amino acids. H. Phillips and W. Criddle (Britain), with the aid of a combina - 
tion of paper chromatography and isotope dilution techniques, identified and evaluated the amount of decomposition 
products (labeled with C' isotope) of gamma -irradiated saccharides. In a report by J. Carroll et al. (Britain), a highly 


sensitive microanalysis method using F’* in fluorine analysis was described; the method is based on the effect of ab- 
sorption of the fluorine ion by a glass surface. 


Several interesting papers were devoted to applications of radioisotopes in organic chemistry. 


In a paper presented by T. Pratt and R. Wolfgang (USA), a theoretical investigation of the process of tritium- 
labeling of methane based on exchange between gaseous tritium and methane, induced by radiation, is presented, and 
a treatment of the reaction mechanism is given. Together with the paper by K. Wilzbach (USA) which takes up gamma 
radiation effects on tritium-labeled toluene, this paper may be viewed as a serious attempt to study the elementary 


events brought about in the molecules of various organic compounds by exchange between hydrogen and gaseous tri- 
tium. 


In papers presented by H. Simon andD. Palme (West Germany), E. Hondett and R, Rowton (USA), G. Ropp (USA), 
D. Marinsky and D. Krasner (USA), G. I. Avdulov, N. K. Semenova, and I. F. Tupitsyn (USSR), the use of isotope effects 
of hydrogen and carbon in studying the mechanisms of organic reactions is described, The study of thermodynamic and 

kinetic isotope effects of hydrogen in isotope exchange of hydrocarbons dissolved in liquid hydrogen bromide and 

liquid ammonia reported by the authors of the last-mentioned paper aids in promoting a correct choice of methods 


for producing tagged (deuterated and tritiated) aromatic hydrocarbons and for bringing about isotope exchange condi- 
tions, 


A paper presented by L. de Boekner and associates (Belgium) gives findings of measurements of the diffusion rate 
of a polymer solution through a semipermeable membrane and of adsorption of the polymer on the surface of the 
membrane. The paper shows that the possibilities inherent in tracer studies of the properties of polymeric systems 
have been greatly underestimated, and that tracer methods enjoy unquestionable advantages over the conventional 
and much less sensitive (e.g.,gravimetric) methods, 


A large group of papers were devoted to the various applications of isotopes in physical chemistry (investiga- 
tions of kinetics and reaction mechanisms, of the structure of complexes, ion and isotope exchange, adsorptive pheno- 


mena, radiolysis, solubility, corrosion, etc.). Many of these papers hold considerable theoretical and methodical inter- 
est for specialists working in those fields. 


For example, a paper by H. Lawrence and D, Stranks (Britain) expounds a new method for studying reactions of 
halide atoms, using the nuclear reaction (n,y ) to produce radioactive atoms, which thereupon enter into a reaction 
with the substrate, Values were obtained for the rate constants and energy of activation for exchange between iodine 
and methyl iodide ions. A paper by E, Kereszk et al. (Hungary) also took up exchange between iodinated organic 
molecules and the element iodine, particularly the effect of solvents on isotope exchange rate. Studies of this kind 
have paved the way for evaluating the strength of iodine atom bonds in molecules, and may harbor practical signifi- 
cance for production of tagged organic (iodinated) compounds through the isotope exchange method. 


An investigation of the structure of complex compounds through exchange of radioisotopes was discussed in a 
paper by B. Jazowska-Zebiatowska and J. Czolkowski (Poland). In this work, isotope exchange methods were success~ 
fully combined with other physical-chemical techniques. Data on exchange rates for carbon isotopes in the KgMn(CNg) 
+KC™ N and KsMn(CN)sNO + KCN systems allowed the authors to present justification for their view on the elec- 
tronic structure of complexes, and to prove the associative character of the exchange process. The fact of the effect 
of visible light on exchange rate will be found interesting. 
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A paper by R. Dietz and associates (USA) presented data characterizing the process of extraction of In’ and Ga’ 
ions from aqueous HC] solutions by dichloroethyl ether. The extraction coefficients were derived radiometrically. In 
addition to new results referable to the ionization constants of HCl and HBr in an organic phase and the instability 


constants of InCl4 and InBr4, the paper present such practical value as, for example, elaboration of ore extraction 
technology for trivalent metals. 


Radioisotopes in Industry 


Over 300 papers were presented on the subject of industrial isotope applications. These papers took up both tracer 
applications and various arrangements for utilizing the penetrating power, absorption, and scattering of radiation. 


Such problems as friction, wear, and lubrication were also treated. Reports submitted by the USSR, Great Britain, 
USA, Czechoslovakia, Japan, Canada, France shed quite a bit of light on new methods using radioisotopes to probe 
processes responsible for wear of materials, and to reveal the laws governing those processes. A large share of the 
attention was given to problems related to the operating mechanism of lubricants, rate at which lubricants are used up 
in practice, etc, Under this heading, the only papers by foreign scientists which presented any fundamentally new 
data were those by R. Campbell and L. Gruenberg (Britain). Two papers, by G. Cartledge (USA) and A. Findeis (Britain) 
were review papers on the subject, while the remaining papers offered either concrete cases of applications for already 
familiar radioisotope techniques to the solution of particular problems, or failed to present adequate proof that the 
methods suggested in the papers are superior to methods now in use. In particular, one of the papers by Campbell and 
Gruenberg presented experimental data on engine wear obtained by the method of radioactive inserts, These data, a- 
long with findings from numerous other investigations resulting from use of the method of radioactive inserts practiced 
in the USSR, place in doubt the feasibility of employing the method of activating large-size diesel engine cylinder 
sleeves, a technique evolved in the USA, especially in the light of a comparison of the simplicity and ease of activa - 
tion of inserts as against the complexity of shielded irradiation arrangements, 


Much interest was evoked by a paper submitted by Yu. S. Zaslavskii and G. I. Shor (USSR) on the use of tracer 
methods to evaluate the performance properties of lubricating oils. Data were presented on a radiometric method for 
investigating electrokinetic processes and sedimentation in oils with additives to which tagged carbon black had been 
introduced, The technique is based on the use of beta counters performing the double function of detectors and elec- 
trodes establishing the electric field, an arrangement which makes it possible to automatically record the rate of dis- 
placement of tagged carbon black in oils with additives in response to electric and gravitational fields. The method 
is suitable for investigating various disperse systems having a labeled disperse phase, The authors also render an ac- 
count of a complex research method which they worked out for probing multicomponent antiseize compounds with 
highly efficient action against scoring, in the absence of corrosive attack, 


An interesting point appearing in papers by Campbell and Gruenberg is the use of pulsating current to achieve 
local heating of surfaces. The experimental data obtained are of interest to the oil refining industry. Kol, Kopecky, 
and Frynda (Czechoslovakia) explain an interesting and simple procedure for investigating wear on large turbogenera - 
tor parts, in the paper they presented. 


W. Grimes (USA) and colleagues presented experimental data on the corrosive action of fluorine -containing 
compounds on metals. A report by G. Cartledge (USA) devoted to a study of corrosive attack on metals and ways of 


coping with it attracted the interest of specialists in several branches of industry, even though the paper was intended 
as a survey, 


J. Cameron and P, Berry (Britain) explained their method for checking the leaktightness of vessels by using Kr® 
which, in the view of the authors, is eminently suitable for leak detection, superior in that respect to other radioac- 
tive gases, There are two variants to the method, In the first variant, the gas is admitted inside the vessel and out- 
leakage is measured, In the second variant, the gas is allowed to surround the component to be measured, The techni- 
que may be used in industry for automatic monitoring of small vacuum volumes, although contamination difficulties 
may turn out to be a problem (e.g., occlusion of the gas by the test component), 


Several papers were devoted to uses of beta-excited x-rays in the solution of various problems. Among these, 
we may note the report by R. Pegg and D. Pollack (Britain) dealing with use of a tritium source for continuous analysis 
of sulfur in hydrocarbon media and gases. The authors developed a device in which absorption of x-radiation emena- 
ting from tritium is put to use, The instrument may be incorporated into a process control loop. The sensitivity of 
the instrument is 0.1% of the sulfur present to an accuracy of 0.03% in a gas and 0.5% of the sulfur present to an ac- 
curacy of 0,04% in kerosene. The instrument will bring in sizable economic returns according to the authors’ estimates. 
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Several reports took up the application of radiofstopes in thickness gaging and weighing of the outside skin 
or jacketing of various industrial products. A paper presented by Soviet authors B. N. Vasil'ev and V. K. Latyshev et al. 
on radioactive thickness gaging of hot-rolled sheet product (5-50 mm thickness) was heard with great interest. 


Let us mention, finally, the paper by J. Gregory devoted to the latest achievements in isotope applications in 
Australia. This paper presents a description of research efforts completed during the past two years. The isotope 1'*! 
was employed to label a small quantity of hot water at the intake of the cooling pond of an electric power generating 
station, and the movement of this water was traced till it was completely mixed with the water in the pond. Scintilla- 
tion counters placed below water level recorded the radiation. A recording scintillation detector mounted in the pump 
house of the power station measured the rate of water transport shown by the labeled water returning to the station pond. 
The paper also contains a description of an investigation employing Br® isotope to probe the flow paths in a raw-sugar 
tank. An application of P*?-tagged lubricant to uncover the reasons for the troublesome appearance of oil spots on 


valuable railway cargo was described, as well as the performance of a high-precision liquid density gage measuring 
fluid flow in concrete piping. 


Production of Isotopes and Radiation Sources and Synthesis of Labeled Compounds 


These questions were discussed in 27 papers embracing various aspects of the problem of radioactive prepara - 
tion fabrication. 


Several interesting reports dealt with methods and production technology relating to pile production of radio- 
isotopes and the use of accelerators in isotope production. The report by J. Roy (Canada) devoted to production of the 
isotope Be’ was found highly valuable and stimulating. This paper traces the development of the concept of producing 
radioactive isotopes by using double reactions and successive reactions of the (T, n) and (T, p) types, where tritons are 
first produced by irradiation of lithium using pile neutrons in the reaction Li’(n, p)T. At present, this method is the 
one used to produce F'® and other isotopes. The authors studied four different nuclear reactions and made an evalua- 
tion of the efficiency of alternate pathways in Be’ production. Their suggestion is to produce Be’ via a double reac- 
tion on protons and deuterons obtained through recoil when targets are irradiated by fast neutrons. The reactions are 
B'°(p,)Be”, Lil(a n)Be", and Be°(n,3n)Be”. The authors for instance demonstrated the fact that the reaction 
Li'(p,n)Be" has practical significance in isotope production when LiOH is used as target: 1 gram of LiOH subjected to 
20 days bombardment in a neutron flux of 2:10" neutrons/cm’ sec intensity, yields several microcuries of Be’. T his 
paper is particularly important from a fundamental standpoint, since it demonstrates the feasibility of seeking out ways 
of producing isotopes through double reactions, 


The short-lived isotope Fr’, which may be produced from actinium via a paper chromatograph technique, may 
hold great promise for medical applications (paper by J. Fouarge and V. Meinke, Belgium). The same applies to iso- 
tope P**, The development of a method proposed earlier for producing P® via neutron bombardment of enriched s® 
isotope was discussed ina report by S. Forberg and T. Vestermark (Sweden). 


Under this heading, we also take note ofa Belgian paper (J. Goste and A. Demildt) on cyclotron production of 
carrier-free Ta’® through deuteron irradiation of tungsten, and on the possibility of producing short-lived isotopes 
Ta'*™ (deuteron energy 12 Mev and higher) and Ta’® (deuteron energy 22 Mev) at satisfactory yields. H. Morinaga 
told of a new technique elaborated in Japan for producing radioisotopes of high specific activity on an internal betatron 
target by means of 25 Mev bremsstrahlung and at a dose rate of ~10* r/min. Eight new radioactive isotopes were 
obtained in this way: Co®, Ga", As®!, Tm", 


Considerable interest was stimulated by a paper presented by T. Butler, E. Lamb, and A. Rapp (USA) on new 
technological radioisotope production processes at Oak Ridge. These new processes are similar to techniques developed 
in the USSR. A shift to solvent extraction for separating out fission products is characteristic and significant in that 
solvent extraction was not used for those purposes previously in the USA. It is noteworthy that the processes recently 
developed in the USA have already been put into practice and are being utilized to produce isotopes having activities 
ranging high into the kilocurie range (cs™, cel, pm!®, sr™, To), It is particularly fitting to note the scaled-up 
pilot production of Pm’ and Tc” and the high purity and specific activity reported for Pm’ (425 curies per gram 
of oxide, which is one half the highest theoretical computed activity). 


Reports by I. Filosofo et al, (USA) and J. Cameron and J. Rhodes (Britain) on the production and investigation 
of beta-excited soft x-ray sources are of considerable interest. In these papers, data on the spectra and photon yields 
are determined experimentally and computed theoretically as functions of energy of beta radiation, thickness and 
atomic number of the target material; the design of four different types of sources is shown (target in front, reflecting 
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target, target formed of layers, and a source where the radioactive substance and target material are mixed). The 
papers present data for H®, Kr®, sr+¥®, 7174, and pm’ targets, and describe a high level Pm! source for industr- 
ial radiographic work. 


J. Puig and J. Sandier (France) presented a paper on the preparation of solid Kr® sources. A team of Soviet 
scientists (N. I. Bogdanov and others) submitted a report on a thermal method for preparing Sr™-base sources. 


Several reports read by scientists from various countries (USA, Italy, UAR, Belgium, West Germany) on the sub- 
ject of synthesis of labeled compounds by various physical-chemical techniques (in an electric-arc discharge, under 
exposure to radiation and ultraviolet light) were heard at the conference. All ofthese papers offer only the underlying 
principles and possibilities inherent in the methods, with no concrete synthesizing techniques, since most of the workers 
failed to carry out adequate separation and purification of the substances containing radioisotopes, and therefore failed 
to adduce any data on contaminations of the basic product, It should also be mentioned here that these efforts constitute 
a novel and promising trend in the production of labeled molecules, 


Soviet scientists presented several reports on research efforst using the energy of recoil nuclei to obtain radio- 
active preparations, 


In a report by L. N. Kurchatova and B. V. Kurchatov, the production of concentrated cl® and Br® preparations 
(with increased specific activity) via the method of recoil nuclei under irradiation in high-density neutron fields was 
discussed. The authors suggested and tried out an ingenious method for producing concentrated preparations of radio- 
halogens based on a modification of the Szilard-Chalmers effect. As a result, Br® preparations were produced with 
specific activity reaching 15 mC/mg, and Cl® preparations with specific activity to 0.23 mC/g. 


A. N. Nesmeyanovy et al, discussed the results of an investigation of interaction between tritium recoil atoms 
(obtained from irradiation of lithium) and some organic compounds. The authors studied the effect of tritium recoil 
atoms on cyclohexane, methyl cyclohexane, cyclohexanol, cyclohexylamine, benzene, cyclohexadienes, and demon- 
strated the possibility of producing tritiated cyclic hydrocarbons via “hot synthesis” with rather high yields. A report 
by B. G. Dzantiev and N. M. Barkalov entitled "Synthesis of radiosulfur-labeled compounds as a result of reactions 
between hot S*° atoms and cyclic hydrocarbons" summarized the results of many years’ work by the authors on synthe - 
sizing sulfur-tagged open-chain and cyclic mercaptans and sulfides. 


The role of vapor chromatography in the synthesis of labeled organic compounds was stressed in a paper sub- 
mitted by Pichat et al, (France). A paper by R. Henri et al, (France) described methods for synthesizing some organic 
iodine-containing compounds by the isotope exchange approach. 


A paper by Soviet authors O. I. Andreeva and G. T. Kostikova discussed carbon isotope exchange reactions 
between potassium cyanide and lithium carbonate, sodium carbonate, potassium carbonate, and barium carbonate. * 
The authors worked out a new method for synthesizing C'*-tagged potassium cyanide, which is the basic building 
block in the synthesis of many organic C compounds via isotope exchange BaC4O3->KCN with subsequent separation 
of the potassium cyanide from the carbonate by ammonia liquid extraction. Labeled KCN is obtained in up to 90% 
yield and containing up to 96% basic material. Finally, we take note of the fact that work is being pursued in other 
lands on development of biosynthesis approaches to the production of labeled compounds, Interesting data on biosynthe- 
sis of three important s*-jabeled amino acids are to be found in a report by P. Chapeville et al. (France), S. A. 

Barker and associates (Britain) described the synthesis of two complex compounds (cyclic derivatives of C™-tagged 
saccharides) based on biosynthesis by microorganisms. 


Other papers described methods of synthesis involving the use of various labeled organic compounds, 


The proceedings of the conference will be made available in printed form by the International Atomic Energy 
Agency early in 1961, 


* Translators note “boron carbonate" in text. BaCOg is probably intended, and is mentioned in next sentence. Boron 
does not form a stable carbonate, as far as I know. 
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BRIEF COMMUNICATIONS 


Translated from Atomnaya Energiya, Vol. 10, No. 2, pp. 190-191, February, 1961 


USSR. The first steam turbine for the Novo-Voronezh nuclear power station has come off the line at the Kirov 
factory in Kharkov. The installation has a power rating of 70,000 kw, The turbine, designed for 3000 rpm, will be 
operated on saturated steam at 29 atmos pressure. The turbine will be automatically controlled. 


USSR. The IV All-Union Conference on Physico-Chemical Analysis met in Moscow, December 6-10, 1960, 
Reports dealing with physico-chemical analysis of metals used in nuclear engineering, and with new problems in 
physico-chemical analysis related to nuclear radiation effects on materials were discussed at the plenary sessions. 
Panel sessions took up uranium metals and rare (refractory and rare-earth) metals. 


189 


BIBLIOGRAPHY 


NEW LITERATURE 


Translated from Atomnaya fnergiya, Vol. 10, No. 2, pp. 192-200 


Books and Symposia 


I. E. Irodov. Compendium of Problems in Atomic Physics. 2nd edition. Moscow, Atomizdat, 1960, 252 pages. 
50 kopeks.* 


This compendium contains about 850 problems and fairly adequate hints for the solution of the most complex 
and tough ones. 


Each chapter begins with a brief review of the basic concepts and relationships needed to solve the problems 
dealt with in the chapter; the book terminates with a reference appendix of physical constants and useful tabulated 
data, 


The number of problems has been increased 50% over the number in the first edition of the book, and illustra- 
tions have tripled in number, Three new chapters have been introduced: the Schroedinger equation, Characteristics of 
the atomic nucleus, and Elementary particles, The appendix has been supplemented with new tables. 


This problem and answer book may be used as a textbook for a general course in atomic physics, for students 
majoring in physics. 
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Edited by Academician B. N. Petrov. Moscow, Atomizdat, 1960, 313 pages, 1 ruble, 41 kopeks.* 


This book makes a presentation of the basic information on technical means of automatic control, and describes 
up-to-date techniques in design work and calculations for automatic checking, monitoring, and control systems applic- 
able to technological processes in the hydrometallurgy of rate and radioactive metals. 


The fundamental tenets of linear automatic control theory are presented, and some of the techniques are ap- 
proached on the basis of concrete examples in setting up control loops and studying their stability. Some of the methods 
employed in studying nonlinear automatic control systems are entered into briefly. Process control variables of interest 


in hydrometallurgy are discussed, and recommendations are offered on the choice and adjustment of automatic control- 
lers and on the design of control components. 


Particular aspects of information on the design of systems subjected to random distrubances and adaptive control 
systems are outlined in concise form. 


Recommended reference literature is listed at the end of the book. 


The book is written to reach engineering and technical workers engaged in hydrometallurgical production, but 


may also be found useful by students in advanced institutes specializing in the field of hydrometallurgical process 
control. 


M. Curie. Radioactivité, Translated from the French, 2nd edition, Moscow, Fizmatgiz, 1960, 516+ xvi pages. 
2 rubles, 2 kopeks. 


This book is a classical monograph written on the basis of a lecture course which the author taught for many 
years at the Sorbonne. The first section of the monograph offers brief information on gaseous ions, electrons, and 


* Prices are in terms of the new “heavy” ruble. 
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rays formed in rarefied gases as electric current is passed through them. The second and principal section of the mono- 
graph is devoted to naturally radioactive elements and their emitted radiations. 


The book is written both for practicing scientific workers (physicists and chemists) and for students. 


A. H. Vapstra, G. I. Hiech, R, Van Lishut. Nuclear Spectroscopy Tables.Translated from the English. Moscow, 
Atomizdat, 1960, 178 pages. 85 kopeks. 


The book contains reference material dealing for the most part with alpha, beta and gamma ray spectroscopy. 


Most of the data is presented in the form of graphs and tables, and only a small portion appears in textual form 
as concise explanations included in each section. The experimental data appearing in the book are accompanied by 
information of a theoretical nature, Particularly interesting from this standpoint are sections VIII and IX, which are 
devoted to angular distribution problems, correlations, and nuclear models, 


The book may be used as a reference handbook in nuclear physics labs, 


Plutonium Recycling in Thermal Reactors, Translated from the English. Moscow, Foreign Literature Press, 1960, 
112 pages, 40 kopeks, 


Problems involving plutonium recycling in nuclear reactor systems had to crop up sooner or later in nuclear 
engineering discussions. These problems form the subject of one of the 13-volume USAEC publications, 


This volume embraces a rather wide range of questions, starting with physical problems related to plutonium 
utilization, and ending in teactor design for plutonium-burning reactors, It is generally known that one such reactor 
(PRTR) is under construction in the USA, and is expected to go into operation in the not too distant future. The problems 
of plutonium reactor physics are considered in this book in their relation to the PRTR reactor. Data of possible use in 
reactor neutron physics calculations are presented. However, the lion's share of the attention is given to plutonium 

fuel cycling, which is approached from the vantage point of maximizing savings in fuel costs and consequently optim- 
izing the over-all economicsof the reactor plant as a whole. 


In discussing the problem of determining plutonium costs, we are led to the conclusion that this cost can be ar- 
rived at by assuming the plutonium to be a source ofheat required for the generation of electric power, in other words 
through a comparison between plutonium and uranium as neutron sources, 


The second large grouping of questions relates to plutonium technology with plutonium approached as a fuel. 
In this section, the design of fuel elements, plutonium alloys and their properties, plutonium ceramics, the fabrica- 
tion of plutonium fuel and plutonium fuel elements, and data on pile testing of plutonium fuel are discussed. Methods 
used in chemical reprocessing of plutonium fuel elements are described, 


The third section is devoted to the design of facilities for handing plutonium or using plutonium, In particular, 
a pilot plant for fabricating plutonium fuel and a plutonium recycle test reactor are described. A report on this reac- 
tor was presented to the second Geneva conference in 1958 and appeared in Russian translation.* 


On the whole, the book may prove highly useful to specialists in power reactor design, since it sheds light on 
problems directly related to the perspectives of nuclear power development. 
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and the effect of radioactive radiations onthe productivity of spring wheat crops. 
Priborostroenie No. 10 (1960) 
A. V. Klimushev, V. S. Merkulov, 30-31, Criterion in the choice of protective membranes in measurements 
using the method of attenuation of beta radiation, 
Trudy inst. biologii (Ural'sk. filial akad. nauk SSSR) No. 13 (1960) 
N. V. Timofeev-Resovskii, 73-96, Development and present state of the art of radiation genetics, 
Trudy Sverdlov. sel'sko-khoz. inst. vol. 7 (1960) 
P. R. Borodin et al., 239-43, Some results of the effect on ionizing radiations on corn. 
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Atomkernenergie 5, No. 9 (1960) 
W. Burkhardt, D. Herrmann, 324-32, Calibration of beta dosimeters, 

R. May, H. Schneider, 333-35, Artificial radioactivity of spring water and river water. 

W. Gerlach, K. Stierstadt, 335-40. Studies of radioactive fallout and radioactivity of air (VI). 


Atomkernenergie 5, No, 10 (1960) 
K. Lindackers, 379-83, Reactor radiation shielding design, II. 


Atompraxis 6, No. 9 (1960) 
R. Plesch, 348-52, Shielding design calculations observing existing radiation regulations, 
Atompraxis 6, No, 10-11 (1960) 


N. Hemmer, H. Wilimzig, 392-96, Scram monitoring system for nuclear reactors, 
T. Musialowicz, F. Wachsmann, 404-407, Note on the problem of beta shielding. 


H. Egelhaaft et al., 407-11. Use of the film blackening method for checking observation of radiation shielding 
measures in laboratories where radium and radioisotopes are handled, 
F, Wachsmann, 411-413. Note on the preceding article. 


F. Habashi, T. Shonfeld, 414-15, Gamma-ray spectrometry techniques for detecting settling of “fresh” fallout 
fission products into the biosphere. 


Atomwirtschaft V, No. 9 (1960) 
J. Pfaffelhuber, 384-90. A new law on radiation shielding in West Germany. 
Atomwirtschaft V, No. 10 (1960) 


L. Merz, 437-41, "Uncontrolled runaway” and “scramming” as shielding problems in nuclear electric power 
stations, 

N. Hemmer, H. Wilimzig, 449-57, Scram monitoring system for nuclear reactors. 

G. Boucke, 463-65, Counters in radiation measurements techniques. 

P. Rinn, 465-72, Trends inscintillation counter development. 

J. Schulz, 473-74, Use of transistors in dosimeters. 

S. Miiller, 475-78, Modular assembling of dosimeters. 


Kerntechnik 2, No. 9 (1960) 
M. Oberhofer, 276-81, Graduation and calibration of dosimetric instruments. 
Kerntechnik 2, No. 10 (1960) 


T. Friese, 322-23, Shielding devices in the power measurements hole of the Berlin reactor. 


H. Fessler et al., 324-326. Simple anticoincidence circuits with two proportional flow counters. 
E. Leiter, 327-29. A portable fast-neutron dosimeter. 


Nucl, Energy No. 149 (1960) 
A. Gould, 464-67, Radioactive contamination problems. Required safety measures, 


Nucl, Energy No. 150 (1960) 
J. Coltman, 530-532. Structural materials for thermionic generators, 


Nucl, Engng. 5, No. 52 (1960) 
K. Orton, 401-402. Review of papers read at the conference on radiation dosimetry, Vienna, June, 1960. 
Nucl, Power 5, No. 55 (1960) 
D. Taylor, 104-108, International conference and exhibit on instrumentation and measurements. 


Nucl, Sci. and Engng. 8, No, 3 (1960) 
D. Anderson, K. Shure, 260-69, Thermal neutron flux distribution in metal-hydrogenous shields. 
S. Wexler, 270-73, Standard boron solutions for neutron absorption measurements, 
J. Ayer et al., 274-76, Engineering aspects of the water vapor permeability of glove materials. 
Nukleonika V, No. 9 (1960) 
Z. Kachliki, 569-73. A linear pulse amplifier. 
J. Andruszkiewicz et al., 576-82, Plastic scintillators. 
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Vv. Radioactive and Stabie Isotopes 


Labeled -atoms technique. Uses of radioactive radiations. Direct conversion of nuclear energy to electric energy 


Avtomobil. prom. No, 9 (1960) 
A. Kh, Eliava, 30-33, Procedure for utilizing tracers instudies on wear on carburetor engine parts, 


Azerb. khim, zhur, No. 2 (1960) 
M. M. Melik-Zade et al,, 127-34, Evaluation of film formation of oil additives on rubbing surfaces, using tracer 
methods. 


Vestnik akad, nauk SSSR XXX, No. 10 (1960) 
V. D. Goroshko, 64-66, Radiometric dry method for coal dressing. 


Vestnik vsesoyuz. nauchno-issled. inst, zhelez.-dorozhn, transporta No. 6 (1960) 


V. S. Shadikyan, A. N. Toropchinov, 13-17, Procedure for estimating the wear-resistance properties of diesel oils, 
using tracer techniques, 


Vestnik Mosk. univ., seriya 3, Fizika, astronomiya. No. 4 (1960) 
B. A. Nelepo, 64-70, Determination of the turbulent diffusion coefficient in sea water, using tracers. 


Voprosy rats. vedeniya vinogradarstva No, 1 (1960) 
P. I. Litvinov, 89-93, Radiophosphorus used in studying regeneration of grapevine roots and interrelations between 
separate parts of a grapevine. 


Doklady akad,. nauk SSSR _ 134, No. 2 (1960) 
£. E. Vainshtein, Yu. I. Belyaev, 322-25, Tracer studies of spatial distribution of atoms in a constant-current 
arc plasma in different atmospheres. 
I, I. Pokrovskii, M. M. Pavlyuchenko, 391-393, s®§-tracer investigation of the oxidation mechanism in oxida- 
tion of copper by molten sulfur. 


Doklady akad. nauk SSSR 134, No. 3 (1960) 
V. F. Oreshko et al., 636-38, Effect of ionizing gamma rays on the structural and mechanical properties of 
starch gels. 


Doklady akad. nauk SSSR 134, No. 4 (1960) 


V. E. Kazarinov, N. A. Balashova, 864-67, Tracer investigation of adsorption and desorption of iodine on smooth 
platinum. 


Izvestiya akad, nauk Armyan, SSR. Biol. nauki 13, No. 5 (1960) 
V. L. Ananyan, 87-90, On the radiometric method for determining total potassium content of a soil. 


Izvestiya akad. nauk Latv. SSR No, 7 (1960) 
G. Gaile, 162-71. On experience acquired in the assimilation of radioactive isotopes and nuclear radiations 
into the operations of enterprises in the national economy of the Latvian SSR, 
Izvestiya vyssh. ucheb. zaved, Fizika No. 4 (1960) 
—-—-, 152-59. Effect of gamma radiation on the dielectric properties of some electrical insulator materials. 
3. K, I. Vodop"yanov et al, Lacquers and varnishes. 4. K. 1. Vodop"yanov et al. Polyethylene. 
Izvestiya vyssh, ucheb, zaved. Chernaya metallurgiya No. 8 (1960) 
N. A. Chelyshev et al,, 48-58, Tracer studies of metal flow during rolling on a 750 rolling mill. 


Liteinoe proizvodstvo No. 10 (1960) 
V. I. Dobatkin, B. I. Simakovskii, 34-36, Tracer investigation of aluminum crystallization in continuous cast- 
ing. 
Mineral’noe syr*e (Vsesoyuz, nauchno-issled. inst. mineral’nogo syr'ya) No. 1 (1960) 
E. I. Zheleznova, D. V. Tokareva, 277-82. Use of a scintillation beta counter in determining zircon content. 
Prom, énergetika No. 9 (1960) 
A. D, Grober, 33-34, Applications for radioactive level gage annunciators in the cotton processing industry. 
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Trudy vsesoyuz. nauchno-issled. kinofotoinst. No. 37 (1960) 
A. V. Borin et al., 78-94, Tracer studies of adsorption of thiosulfate by microcrystals of silver bromide in photo- 
graphic emulsions, using radiosulfur tracer atoms, 


Trudy Mosk. inst, neftekhim, i. gaz. prom. No, 31 (1960) 
K. I. Yakubson, Sh. A. Guberman, 81-99, Experimental investigation of the space and energy distribution of 
neutrons, using a plastic model. 
Trudy. Kharkov. avtomobil.-dorozhn. inst. No. 23, (1960) 
A. I. Danilenko, 125-28. On the use of radiometric techniques in highway construction, 


Uspekhi fiz. nauk LXXII, No. 2 (1960) 
N.S. Gorbunov, V. I. Izvekov, 273-306, Tracer study of diffusion in metal oxides. 


Khim. i tekhnol. topliv i masel No. 9 (1960) 
A. A. Kuznetsov, V. S. Luneva, 61-64, Quantitative tracer determination of anticorrosion properties of high- 
consistency lubricants. 
Atomkernenergie 5, No, 10 (1960) 
—-— , 384, Possibilities of utilizing underground explosions for scientific and industrial purposes. 


Atompraxis 6, No, 10-11 (1960) 


H. Langel, 337-379, A summing instrument for comparing values obtained in tracer measurements of the thick - 
ness of materials. 


P, Gerke, 380-85, Applications for radiation-measuring instruments:in the petroleum and petrochemical industry, 
K. Waechter, 385-91. Choice of a proper detector for gaging the thickness of relative thick specimens. 


Atomwirtschaft V, No. 9 (1960) 
H. Krauch, 411-13, Radiochemistry in engineering and industry, II. 


Kerntechnik 2, No. 9 (1960) 
W. Kihn, F. Herrmann, 268-70. Possibilities of utilizing tritium bremsstrahlung in industrial applications. 


——-—, 285-86, A new cobalt facility for weighing commercial freight cars. 


Kerntechnik 2, No. 10 (1960) 
——-, 331. Use of radioactive tracers in West Germany. 


Nucl, Energy No, 149 (1960) 
—-—-, 490, ff. Copenhagen international conference on the uses of radioactive isotopes inscience and industry. 
Review of papers on industrial applications. * 

W. Arnold, V. Stonehocker, 492-93, Tracer applications for determining wear on reciprocating piston engines, 


Engng. 5, No. 52 (1960) 
—-—-—, 517-22. Radiography applications for accelerators, 


*See article, this issue JAE, page 182, 
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3 titles of immediate interest! 


LECTURES ON NUCLEAR THEORY 
by L. Landau and Ya. Smorodinsky 


Translated from Russian 


A concise presentation by these world-famous Soviet physicists of 
some of the basic concepts of nuclear theory. Originally published 
at $15.00 per copy. 


. a real jewel of an elementary introduction into the main con- 
cepts of nuclear theory . . “,-NUCLEAR PHYSICS 

“. .. @ decidedly worth-while addition to any experimental nuclear 
physicist’s library.”—E. M. Henley, PHYSICS TODAY 


cloth 108 pages $5.25 
Order from: PLENUM PRESS 227 W. 17th Street, New York 11, N.Y. 


— of major interest to all researchers in low-temperature physics. 


A SUPPLEMENT TO “HELIUM” 
By E. M. Lifshits and E. L. Andronikashvili 


Translated from Russian 


This notable volume consists of two supplementary chapters, by 
these outstanding Soviet physicists, which were added to the Rus- 
sian translation of W. H. Keesom’s classic book “Helium.” 

The first chapter, by Lifshits, presents a concise resume of the 
Landau theory of superfluidity. The second chapter reports in con- 


siderable detail the experimental work conducted by Peter Kapitza 
and E. L. Andronikashvili. 


cloth 167 pages illustrated $7.50 


Order from: CONSULTANTS BUREAU = 227 W. 17th St. 3 New York 11, N. Y. 


A NEW METHOD IN THE THEORY OF 
SUPERCONDUCTIVITY 
By N. N. Bogoliubov, V. V. Tolmachev and D. V. Shirkov 


Translated from Russian 


The Soviet authors put forth a systematic presentation of a new 
method in the theory of superconductivity, developed as a result of 
the research of N. N. Bogoliubov and V. V. Tolmachev—based on 
a physical and mathematical analogy with superfluidity. This new 
method is an immediate generalization of the method developed by 
Bogoliubov in formulating a microscopic theory of superfluidity. 
The authors give calculations for the energy of the superconducting ground state using 
Frohlich’s Hamiltonian, as well as of the one-fermion and collective elementary excited 
states. A detailed analysis of the role of the Coulomb interaction between the electrons in 
the theory of superconductivity is included. The authors also demonstrate how a system 


of fermions is treated with a fourth-order interaction Hamiltonian and establish the cri- 
terion for its superconductivity. 


cloth 121 pages illustrated $5.75 


Order from: CONSULTANTS BUREAU 227 W. 17th St. @ New York 11, N. Y. 


Descriptive folders upon request. 
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Research by Soviet Experts 
Translated by Western Scientists | 


Soviet Research on the LANTHANIDE 
AND ACTINIDE ELEMENTS, 1949-1957 | 


An important contribution to the literature of nuclear chemistry, this collection of papers is a 
comprehensive presentation of Soviet research on the chemistry of lanthanides and actinides. The 106 reports 


included in this collection appeared in the major Soviet chemical journals translated by Consultants Bureau, 
as well as in the Soviet Journal of Atomic Energy, 1949-1957. 


The five sections, totalling 657 pages, provide broad representation of contemporary Soviet _ 
research in this important aspect of nuclear science. This’ collection should be accessible to all nuclear 
researchers, whether theoretical or applied. 


Each part may be purchased as follows: 
Basic Chemistry (25 papers) $15.00 
Analytical and Separation Chemistry (30 papers)............ $20.00 
Nuclear Chemistry (and Nuclear Properties) (32 papers) ...$22.50 
Nuclear Fuel Technology (9 papers) ..............ee-eeeeeee $7.50 


RADIATION CHEMISTRY, 
PROCEEDINGS OF THE FIRST 
ALL-UNION CONFERENCE MOSCOW, 1957 


More than 700 of the Soviet Union’s outstanding research scientists participated in this 
conference sponsored by the Academy of Sciences and the Ministry of the Chemical Industry. Each of the 
56 reports read in the various sessions covers either the theoretical or practical aspects of radiation chemistry, 
and special attention is given to radiation sources used in radiation-chemical investigations. The general 


discussions which followed each report and reflected various points of view on the problem under analysis 
are also included. 


Primary Acts in Radiation Chemical Processes 


heavy paper covers 5 reports, plus discussion ...... illustrated .....-... $25.00 
Radiation Chemistry of Aqueous Solutions (Inorganic and Organic Systems) 
heavy paper covers 15 reports, plus discussion ...... illustrated ......... $50.00 
Radiation Electrochemical Processes 

heavy paper covers 9 reports, plus discussion ...... illustrated .....-... . $15.00 
The Effect of Radiation on Materials Involved in Biochemical Processes 
heavy paper covers 6 reports, plus discussion ...... illustrated ......... $12.00 
Radiation Chemistry of Simple Organic Systems 

heavy paper covers 9 reports, plus discussion ...... illustrated ......... $30.00 
The Effect of Radiation on Polymers 

heavy paper covers 9 reports, plus discussion ...... illustrated ......... $25.00 
Radiation Sources 

heavy paper covers .3 reports ...............eceeged illustrated ......... $10.00 


Individual volumes may be purchased separately. 


NOTE: Individual reports from each volume are available 
at $12.50 each. Tables of contents sent upon request. 


special price for the 7-volume set $125.00 
Payment in sterling may be made to Barclay’s Bank in London, England. 


CONSULTANTS BUREAU 
227 West 17th Street - New York, N.Y., U. ch 
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